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ABSTRACT 


An  experimental  and  theoretical  investigation  was  made  of  transverse 
pumping  saturation  effects  in  ZngY  ferrite  at  frequencies  far  below 
ferromagnetic  resonance.  Chief  interest  was  in  saturation  effects  at 
frequencies  below  the  bottom  of  the  spin  wave  manifold;  where  the  lowest 
order  phonon  instabilities  were  expected  to  compete  favorably  with  the 
higher  order  spin  wave  instabilities . 

Measurements  of  saturation  effects  in  single  crystal  Zn^Y  ,  with 
the  easy  plane  lying  in  the  plane  of  the  disk,  were  made  at  1.32  kMc. 

The  dc  and  rf  magnetic  fields  were  applied  orthogonally  in  the  plane  of 
the  disk.  The  large  rf  magnetic  fields  required  were  obtained  by  employ¬ 
ing  a  strontium  titanate  dielectric  resonator.  In  these  experiments  the 
order  of  the  first  resonant  spin  wave  instability  permitted  was  varied 
from  the  first  to  the  fifth  by  adjusting  the  position  of  the  spin  wave 
manifold  relative  to  the  pump  frequency.  The  nature  of  the  saturation 
was  investigated  by  observing  the  waveform  of  the  reflected  pump  pulse. 

Theoretical  on— resonance  spin  wave  thresholds  up  to  fourth  order 
and  off— resonance  spin  wave  thresholds  up  to  second  order  were  derived. 
The  first  order  theoretical  spin  wave  threshold  was  the  lowest  threshold 
in  the  range  of  low  dc  magnetic  fields.  At  higher  dc  fields  the  second 
order  theoretical  on-  and  off-resonance  spin  wave  thresholds  dominated 
the  higher  order  thresholds.  Taking  the  magnetoelastic  effect  as  the 
dominant  mechanism,  first  order  phonon  thresholds  were  derived  for 
transverse  and  longitudinal  elastic  waves  propagating  in  the  easy  plane 
and  along  the  c-axis.  The  lowest  theoretical  phonon  threshold  was  for 
transverse  elastic  waves  propagating  and  polarized  in  the  easy  plane, 
and  at  frequencies  below  the  bottom  of  the  spin  wave  manifold,  this 
threshold  dominated  the  theoretical  second  order  off  resonance  spin 
wave  threshold.  Theory  was  developed  for  the  growth  of  phonon  insta¬ 
bilities  from  thermal  level  in  order  to  relate  pulae  saturation  effects 
to  the  phonon  absorption. 
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Transverse  pumping  thresholds  observed  experimentally  vere  comparable 
with  the  above  theoretical  predictions  and  strongly  indicate  that  first 
order  phonon  instabilities  dominate  the  higher  order  spin  wave  insta*- 
bilities  at  frequencies  below  the  bottom  of  the  spin  wave  manifold. 

These  results  together  with  the  phonon  growth  theory  made  it  possible 
to  obtain  rough  estimates  for  the  phonon  Q  ,  one  magnetoelastic 
constant,  and  the  exchange  field  of  Zn2Y  in  the  easy  plane.  The 
observed  pulse  saturation  effects  and  the  results  of  two  experiments 
for  two  different  phonon  absorption  levels  indicated  the  presence  of 
nonlinear  elastic  wave  damping. 
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CHAPTER  I 


INTRODUCTION 

The  term  saturation  is  used  in  two  ways  in  discussing  the  properties 
of  ferrites.  The  most  common  use  applies  to  all  varieties  of  ferromagnetic 
materials.  It  refers  to  the  static  state  of  the  average  magnetization 
vector  when  the  atomic  magnetic  moments  are  aligned  in  one  direction  due 
to  the  application  of  a  sufficiently  large  dc  magnetic  field.  Thus  a 
material  is  said  to  be  saturated  vhen  its  magnetization  approaches  a 
constant  value  with  increasing  dc  magnetic  field  intensity. 

The  type  of  saturation  of  special  concern  here  is  the  saturation  of 
the  uniform  precessional  motion  of  the  magnetization.  The  simplest  motion 
of  the  magnetization  obtains  vhen  the  magnetization  itself  is  “Statically 
saturated  by  the  application  of  a  sufficiently  large  dc  magnetic  field. 

A  uniform  precessional  motion  of  the  magnetization  may  be  established 
by  the  application  of  a  transverse  rf  magnetic  field.  Saturation  of  the 
uniform  precessional  motion  is  usually  defined  in  terms  of  a  critical 
rf  magnetic  field  corresponding  to  the  point  at  which  greatly  Increased 
losses  obtain  if  the  rf  magnetic  field  is  further  increased.  These 
losses  are  due  to  Instabilities  arising  from  nonlinearities  in  the 
equation  of  motion  for  the  magnetization. 

The  saturation  of  the  uniform  motion  of  the  magnetization  at  far* 

1  2 

romognetic  resonance,  observed  by  Damon,  and  Bloembergen  and  Vang, 

3 

was  first  explained  by  H.  Suhl.  He  shoved  that  such  saturation  cpuld 
be  explained  on  the  basis  of  excitation  of  unstable  spin  waves  at  the 
pump  frequency,  and  further  shoved  that  a  subsidiary  absorption  having 
a  lover  saturation  threshold  would  occur  vhen  half  of  the  pump  frequency 
was  within  the  spin  wave  manifold.  This  is  due  to  unstable  spin  waves 
at  half  the  pump  frequency. 
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This  early  work  was  concerned  with  spin  wave  instabilities  excited 

by  a  pump  frequency  which  coincided  with  ferromagnetic  resonance  and 

hence  was  above  the  bottom  of  the  spin  wave  manifold.  Recently,  it  has 

been  shown  theoretically  that  higher  order  spin  wave  instabilities  may 

also  cause  saturation  of  the  uniform  motion  when  the  pump  frequency  is 

below  the  bottom  of  the  spin  wave  manifold- -where  the  first  and  second 

4  5 

order  instabilities  are  far  off  resonance.  *  Even  under  these  conditions 
there  still  exist  resonant  acoustic  waves  at  the  pump  frequency  and  at 
one -half  the  pump  frequency,^  and  there  exist  several  mechanisms  by  which 
phonon  instabilities  at  these  frequencies  may  be  excited. ^  It  is  unlikely 
that  such  phonon  instabilities  could  compete  favorably  with  spin  wave 
instabilities  when  the  pump  frequency  is  within  the  spin  wave  manifold 
because  weaker  coupling  mechanisms  are  involved.  Operation  with  the  pump 
below  the  manifold,  however,  allows  only  higher  order  spin  wave  instabilities 
to  exist  at  resonance;  and  the  lowest  order  phonon  instabilities  might 
be  expected  to  compete  favorably  with  higher-order  spin  wave  instabilities, 
since  the  thresholds  generally  increase  with  the  order  of  the  instability. 

The  purpose  of  this  study  is  to  report  a  theoretical  and  experimental 
investigation  of  the  transverse  pumping  saturation  effects  in  fer¬ 

rite  (2n2  Ba2  Fe^  C>22)  at  dc  magnetic  fields  corresponding  to  frequencies 
below  the  bottom  of  the  spin  wave  manifold.  The  object  was  to  investigate 
the  instability  mechanisms  and  to  determine  whether  saturation  is  caused 
by  spin  waves  or  phonons.  Experimentally,  the  pump  frequency  was  fixed 
at  I.32  kMc,  and  the  dc  magnetic  field  was  varied  to  adjust  the  position 
of  the  spin  wave  manifold  with  respect  to  the  pump  frequency. 

The  hexagonal  material  Zn2Y  was  chosen  because  it  has  a  magnetic 
"easy"  plane.  The  large  magnetic  anisotropy  of  7i makes  it  possible 
to  operate  below  ferromagnetic  resonance  with  relatively  low  applied  dc 
magnetic  fields.  In  the  easy  plane,  then,  the  stiffness  of  the  magneti¬ 
zation  is  low  and  it  is  possible  to  achieve  large  uniform  precession  1 
angles  with  a  much  lower  rf  magnetic  field  than  would  otherwise  be 
required. 

The  theory  of  instabilities  pumped  at  frequencies  below  the  spin 
wave  manifold  is  developed  in  Chapters  II  and  III.  The  derivation  of 


thresholds  for  spin  wave  instabilities  up  to  fourth  order  is  presented 
in  Chapter  II.  The  optimum  spin  wave  thresholds  were  obtained  for  each 
order.  Chapter  III  contains  a  derivation  of  the  first  order  thresholds 
for  phonon  instabilities.  Various  possible  instability  mechanisms  are 
discussed.  By  taking  the  magnetoelastic  effect  as  the  dominant  mechanism, 
thresholds  were  determined  for  transverse  and  longitudinal  acoustic  waves 
propagating  in  the  easy  plane  and  along  the  c-axis.  A  comparison  was 
made  between  the  lowest  theoretical  phonon  thresholds  and  the  spin  wave 
thresholds  of  various  orders  for  both  the  resonant  and  nonresonant  cases. 

Chapter  IV  is  concerned  with  the  growth  of  phonon  instabilities  from 
thermal  level.  The  phonon  power  absorption  was  examined  for  both  an  ideal 
pulse  and  the  actual  power  pulse  used  in  the  experiments.  The  threshold 
field,  at  which  the  effect  of  the  instabilities  can  be  observed  by  satuj- 
ration  of  the  reflected  pump  pulse,  is  expressed  in  terms  of  measurable 
microwave  parameters. 

The  microwave  apparatus  and  techniques  which  were  employed  In 
measuring  the  pumping  thresholds  are  described  in  Chapter  V.  Large  rf 
magnetic  fields  were  obtained  by  using  a  strontium  tltanate  dielectric 
resonator.  The  basic  resonator  theory  as  well  as  resonator  and  ferrite 
heating  effects  are  analyzed.  The  methods  used  to  measure  the  unloaded 
Qq  are  discussed. 

Chapter  VI  is  devoted  to  results  of  actual  instability  threshold 
measurements.  Transverse  thresholds  were  measured  in  a  single  crystal 
Zn^Y  disk  with  the  easy  plane  lying  in  the  plane  of  the  disk.  The  dc 
and  rf  magnetic  fields  were  applied  in  the  plane  of  the  disjc.  Results 
from  two  experiments  are  correlated  with  the  phonon  and  spin  wave  thres¬ 
hold  theory  developed  in  Chapters  II  and  III  and  with  the  phonon  growth 
theory  of  Chapter  IV.  In  correlating  theory  with  experimental  results, 
values  were  obtained  for  one  msgnetoelastic  constant,  the  phonon  Q  , 
and  the  exchange  field  of  Z^Y 

At  the  end  of  each  chapter  a  summary  is  given  which  describes  the 
important  features  and  conclusions  for  that  chapter. 
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CHAPTER  II 


SPIN  HAVE  INSTABILITY  THRESHOLDS 

This  chapter  is  concerned  with  the  calculation  of  transverse  pump 
thresholds  of  spin  wave  instabilities  up  to  fourth  order  for  operation 
far  below  ferromagnetic  resonance.  The  order  of  an  instability  n  is 
defined  as 


n 


2m 

“p 


(2.1) 


where  o>  is  the  frequency  of  the  spin  wave  whose  instability  is  being 
considered,  and  o>p  is  the  pump  frequency  of  the  rf  magnetic  field. 

In  the  experiment  described  in  Chapter  VI,  the  position  of  the  bottom  of 
the  spin  wave  manifold  (Fig.  2.2a)  is  varied  with  respect  to  the  pump 
frequency  by  changing  the  strength  of  the  applied  dc  magnetic  field  in 
the  ferrite  sample.  From  (2.1)  and  Fig.  2.2a  instabilities  of  resonant 
spin  waves  are  permitted  only  for  orders  n  such  that  noip/2  lies 
above  the  bottom  of  the  manifold.  Thus  the  order  of  the  first  allowable 
spin  wave  instability  varies  with  the  dc  magnetic  field,  and  proper 
interpretation  of  the  experiments  requires  the  study  of  several  orders 
of  spin  wave  instabilities. 

A.  GENERAL  EQUATION  OF  MOTION 

The  motion  of  the  magnetization  54  is  ultimately  based  upon  the 
motion  of  each  atomic  magnetic  moment.  The  classical  lossless  equation 
of  motion  of  the  average  atomic  magnetic  moment  on  a  macroscopic  scale 
is 

dM(r,t) 

— — -  *  -  7  »<£,t)  X  B(E,t)  ,  (2.2) 

dt 


-  k  - 


where  g  ,  the  total  magnetic  field  at  point  £  ,  results  from  internal 

fields  due  to  the  motion  of  as  well  as  from  externally  applied  magnetic 
fields.  The  gyromagnetic  ratio  y  is  given  by 


7 


(2.3) 


With  g  «  2.00  for  moments  due  only  to  the  spin  of  the  electron,  the 
value  for  yf 2*  becomes  2.8  Mc/Oe.  Losses  will  later  be  introduced  . 
into  (2.2)  through  phenomenological  damping. 

Hie  assumption  is  now  made  that  the  magnetization  is  saturated  by 
the  application  of  a  sufficiently  large  dc  magnetic  field.  It  will 
further  be  assumed  that  the  dc  magnetic  field,  as  well  as  the  value  of 
the  saturated  magnetization,  is  uniform  throughout  the  sample  volume. 
With  these  assumptions  it  is  convenient  to  introduce  the  unit  magneti¬ 
zation  vector  a(£,t)  ,  whose  components  are  the  direction  cosines ,  of 

the  magnetization,  as  follows: 


I Ste/t) 


(2.4) 


M 
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where  M  is  the  saturated  magnetization.  Equation  (2.2)  then  becomes 
s 


&£>*)  *  -  r  ste/*)  x  » 


(a-5) 


where  the  dot  hereafter  symbolizes  for  any  quantity  the  partial  derivative 
with  respect  to  time.  Since  the  unit  vector  preserves  its  magnitude,  the 
x -component  of  the  unit  magnetization  vector  is  given  by 


a  *  (1  -  (o?  +  a?)]1/2 
x  v  y  z 


-5  - 


(2.6a) 


The  x  component  has  been  singled  out  because  the  dc  magnetic  field  is 
applied  in  the  x-direction,  and  deviations  of  the  saturated  magnetization 
from  the  x-axis  will  be  small.  The  y  and  z  components  of  the  equation 
of  motion  may  be  found  from  an  expansion  of  (2.5)  to  be 


* 

a 

y 


a  h  ) 

X  z 


(2.6b) 


a 

z 


-  7(a  h  -  a  h  ) 

x  y  y  x 


(2.6c) 


Thus  (2.6)  constitute  the  equations  of  motion  of  the  magnetization  for  a 
dc  magnetic  field  applied  in  the  x-direction. 

In  obtaining  components  of  the  total  magnetic  field  g  ,  it  is  often 
useful  to  utilize  the  following  relationship  between  the  magnetic  field 
and  the  free  magnetic  energy  E  : 

S  -  -rvaE  >  <2-7> 

s 


where  represents  the  gradient  with  respect  to  the  direction  cosines 

of  the  magnetization.  The  total  magnetic  field  g  comprises  the  applied 
dc  and  microwave  frequency  (rf)  field  terms  as  well  as  terms  dependent 
upon  oc  .  Because  of  these  terms  of  H  dependent  upon  Of  ,  the 
equation  of  motion  [(2.5)  or  (2.6)3  will  be  nonlinear  in  Of  .  However, 
with  the  application  of  orthogonal  dc  and  rf  magnetic  fields  and  small 
deviations  of  Of  about  static  equilibrium,  the  equations  may  be  linearized. 
Under  the  quasi -static  approximation  to  Maxwell's  equations 
(V  •  B  *  0  ,  V  X  H  *  0)  there  is  a  uniform  solution  (uniform  preces¬ 

sion)  to  (2.6)  and  the  appropriate  magnetostatic  boundary  conditions, 

Q 

as  well  as  an  infinity  of  spatially -varying  solutions.  The  spatially - 
dependent  solution  modes  are  quite  complicated  if  the  exact  boundaiy 
conditions  are  matched,  but  usually  it  is  permissible  to  make  the 
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periodic  boundary  condition  assumption  and  to  regard  the  modes  as  plane 
waves  if 


2x/k  «  sample  dimensions  . 


The  uniform  precession  is  then  often  called  the  k  ■  0  mode. 

In  the  plane  wave  approximation  the  general  motion  of  the  magneti¬ 
zation  is  given  in  the  form  of  a  Fourier  series  expansion  of  Of  , 
namely. 


a 

A# 


2b 


a 

~s 


(28a) 


a 

~s 


(2.8b) 


where  this  expansion  Is  strictly  valid  only  for  large  k's  .  These 
plane  wave  modes  are  referred  to  as  spin  wave  modes — hence  the  subscript 
of  0fs  in  (2.8b).  In  the  above,  ^  is  a  function  of  time  only  and 
has  the  form 


(2.9) 


which  in  general  may  be  complex,  where  the  quantity  u>  is  the  frequency 
of  the  spin  waves  and  the  vectors  are  independent  of  time. 

Since  the  spin  wave  component  of  the  magnetization  a  must  be  real. 

It  follows  that 


(2.10a) 


and  therefore 


(2.10b) 


In  general,  0 t  will  contain  pairs  of  oppositely  directed  traveling  waves, 

with  propagation  vectors  k  and  -k  .  If  nonlinear  terms  are  retained 
in  (2.7)  it  will  be  found  that  terms  proportional  to  the  uniform  preces¬ 
sion  0^  will  couple  pairs  of  spin  waves  a^  and  a*k  ,  and  that  this 
coupling  will  lead  to  exponential  growth  of  pairs  of  spin  waves.  For 
each  pair  it  will  be  found  that 


(2.11) 


so  that  the  unstal  le  modes  are  standing  waves. 


B.  MAGNETIC  FIELD  COMPONENTS 

The  total  magnetic  field  g  comprises  the  following  components: 
the  applied  dc  magnetic  field  g^  ,  the  pump  rf  magnetic  field  , 

the  anisotropy  field  ,  the  dipolar  field  g^^  >  "the  demagnetizing 

field  g^  ,  and  the  exchange  field  .  Figure  2.1  shows  the  coordinate 

system  employed  and  gives  the  direction of  the  applied  magnetic  fields 
and  the  propagation  vector  k  .  The  applied  dc  magnetic  field  is  given  by 


(2.12) 
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The  pump  rf  magnetic  field  which  is  normal  to  the  applied  dc  magnetic 
field  and  is  linearly  polarized  is  expressed  as 


(2.13) 


where  h^.  *  b  cos  uipt  ,  h  is  the  maximum  amplitude,  and  aip  is  the 
frequency  of  the  pump. 

The  anisotropy  magnetic  field  arises  from  the  crystalline  anisotropy 

energy,  which,  for  a  hexagonal  qrystal,  has  its  origin  in  magnetic  dipole - 

dipole  interaction.  From  crystal  symnetry  considerations,  the  anisotropy 

o 

energy  of  a  hexagonal  crystal  is  found  to  be  of  the  form 


Ea  -  sin2  8'  +  Kg  sin1*  8'  +  Kj  sin6  8'  +  sin6  8'  cos  60'  , 

(2.14) 

where  0'  is  the  angle  the  magnetization  makes  with  the  "hard"  z-axis, 
and  <h'  is  the  angle  in  the  basal  (x,y)  plane.  For  Zn^Y  the  negative 
constant  in  (2.l4)  is  predominant.  The  energy  then  assumes  the 

uniaxial  form,  and  when  written  in  terms  of  the  direction  cosines  of  a 
becomes 

Sa  -  Kj_  (1  -  C?)  .  (2.15) 


Utilizing  (2.7),  the  anisotropy  "magnetic"  field  is  defined  by 


(2.16) 
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or 


with  *  -2K^/Mg  defined  as  the  an(Lsotropy  field  along  the  "hard" 
c-axls. 

The  demagnetizing  and  dipolar  fields  are  obtained  from  the  quasi - 
static  approximation  to  Maxwell’s  equation  written  as 

V  •  (g  +  4*Me  a)  -  0  (2.17a) 

V  X  g  *  0  .  (2.17b) 

The  quasi-static  approximation  is  fulfilled  to  a  good  approximation  when 
the  dimensions  of  the  ferrite  sample  are  much  less  than  the  free  space 
wavelength.  The  dipolar  magnetic  field,  which  is  the  demagnetizing  field 
for  spin  waves,  is  obtained  by  finding  the  solution  to  (2.17)  for  the 
spin  wave  modes  (2.8b)  which  satisfy  the  periodic  boundaxy  conditions. 

The  result  obtained  is 


V1  &  *  SL 

ttdip  *  -  1  s  —*r  •  •  <2l8> 

k  * 

The  demagnetizing  field  for  the  uniform  (k  *  0)  mode  is  obtained  from 
(2.17),  together  with  the  exact  boundary  conditions.  For  an  ellipsoid 
in  a  uniform  magnetic  field  the  demagnetizing  field  may  be  represented 
by  the  general  expression 


H  OL 
•'O 


(2.19) 


-  n  - 


where  N  is  the  demagnetizing  tensor.  For  the  special  case  when  the 

applied  magnetic  fields  are  along  the  major  ellipsoidal  axes  of  the  sample, 

** 

N  is  diagonal  with  trace  equal  to  unity.  The  field  may  then  be  expressed 
as 


a. 


Ox 


4xM 

s 


a. 


Oz 


(2.20) 


The  remaining  term  of  g  is  the  exchange  field  ,  which  arises 

from  the  quantum  mechanical  exchange  interaction  between  neighboring 
atomic  spins.  The  result  for  a  crystal  with  uniaxial  symmetry  is 


with 


D1  -  Hex||  a2 

D2  *  Vi/  ' 
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where  Hex||  are  scalar  quantities  of  the  order  of  the  Veiss  field 

and  a,c  are  the  lattice  constants  of  the  unit  cell.  It  is  apparent 
from  (2.21)  that  the  exchange  field  becomes  especially  significant  as 
the  wavelength  of  the  propagating  wave  approaches  lattice  dimensions. 

C.  UNIFORM  PRECESSION  AND  SPIN  WAVE  MODES 

The  equations  of  motion  are  obtained  by  inserting  the  total  magnetic 
field  H  ,  whose  components  were  developed  in  the  last  section,  into 
(2.6).  The  resulting  equations  of  motion  may  then  be  expanded  into 
uniform  precession  and  spin  wave  mode  components  by  employing  (2.8). 

The  uniform  mode  is  separable  from  the  spin  wave  modes  since  only  the 
latter  contain  the  spatial  dependence.  The  uniform  mode  equations  become 


“<*  -  *  l<\  +  >A  “ox  +  <»z  -  *z^x 

“oz  *  H0  +  V<“y  *  Hx>  “Ox1  %  *  %  >  <2-22b> 


where 


Also,  the  equations  for  the  spin  vave  mode  are  found  to  be 


'  -  “H0  “kz  +  “kx  *  \  <“ox  “xz  +  °bz  “kx5 

+  <Hx  “Ox  “kz  *  Nz  “oz  “kx> 

(2.23a) 

+  5  (Voz  *  kz  V  <kx  “kx  +  ky  “ky  +  kz  %z> 

*  k2  “Oz  “kx  ->D2  ^  “ox  “kz 

\z  *  “H0  V  7hy  °W  +  “m  (Ny  %  “kx  -  Hx  “ox  “ky> 

+  5  <ky  “Ox  -  kx  V  (k*  “kx  +  ky  %  +  **  <2'23b) 

+  J*i  k2  (“ox  “ky  -  “Qy  V 


Here  nonlinear  terms  containing  the  product  of  the  uniform  precession  and 
the  spin  vave  components  of  Of  are  retained,  but  terms  higher  than  first 
order  in  0^  are  neglected.  This  is  a  good  approximation  since  the 
unstable  spin  waves  must  build  up  from  thermal  amplitudes.  It  is  through 
the  terms  in  (2.23)  containing  the  uniform  precession  that  energy  is 
coupled  from  the  uniform  mode  to  the  spin  vave  modes.  The  latter  eventu¬ 
ally  will  be  found  to  go  unstable  when  the  uniform  precession  angle 
reaches  a  critical  amplitude. 
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Before  the  uniform  mode  (2.22)  and  spin  wave  mode  (2.23)  equations 
can  be  solved.  It  is  necessary  to  utilize  (2.6a)  to  eliminate  the  . 
x-component  dependence  of  the  magnetization.  By  inserting  relation  (2.6a) 
Into  (2.6a),  and  retaining  only  first  degree  terms  in  the  spin  vave 
amplitudes  as  above,  the  following  relations  are  obtained  for  the  uniform 
and  spin  vave  mode  components  of  the  magnetization  along  the  x -direction: 


(2.24a) 


°kx  *  -  <%  “ky  +  °02  •  <2-2),b> 


Since  +  0^  is  less  then  unity,  then  may  be  expanded  in  a 

binomial  series  in  powers  of  0^  and  .  For  purposes  of  solving 

the  spin  vave  equations,  the  expanded  form  of  (2.24)  is  substituted  into 

(2.23 ) ,  thus  obtaining  the  spin  vave  mode  equations  of  motion  for 

and  o^z  .  The  resulting  equations  are  linear  in  the  spin  vave  direc¬ 
tion  cosines,  but  contain  terms  which  are  modulated  by  the  uniform 
precession  direction  cosines  -.  Such  equations  are  of  the  Natthieu-Hill 
type,  and  instability  thresholds  will  be  found.  Since  substitution  of 

(2.24)  into  (2.23)  leads  to  the  appearance  of  higher  powers  of  0^  and 
aQz  ,  instabilities  of  higher  than  first  order  are  to  be  expected. 


D.  SOLUTION  OF  UNIFORM  PRECESSION  EQUATION 

The  solution  of  the  uniform  precession  mode  is  necessary  to  determine 
the  relation  between  the  driving  magnetic  field  and  the  precession  angle 
of  the  magnetization.  The  simple  solution  is  found  by  letting  «  1 
in  (2.22),  giving 

60y  -  *  -i  a0*  +  ’"**  (2.25*) 

*  °»2  %  *  ' 

.  15  . 


(2.25b) 


where 


“1  *  v  -  x +  <NZ  -  V 

A  0 

“2  “  %  +  “fc  (Hy  ‘  V 

and  the  ferromagnetic  resonance  frequency  is  given  by 


01, 


(o^  o>2) 


1/2 


For  the  special  case  of  a  disk  lying  in  the  easy  plane  of  the  crystal, 
the  dmugnetizing  factor  in  the  plane  of  the  disk  is  given  by 


* 


and  the  solution  to  (2.25)  is  found  to  be 


*“i  *i 


h  cos 


Olpt 


(2.26a) 


a0z  *  h  sin  aipt  ,  (2.26b) 

*0 
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where 


For  operation  on  ferromagnetic  resonance  it  would  be  necessary  to 
include  phenomenological  damping  in  these  equations  by  letting  o>q  ♦  Ho  • 
But  with  the  pump  frequency  significantly  below  ferromagnetic  resonance, 
damping  has  negligible  effect  on  either  the  magnitude  or  the  phase  of  Of  . 

In  fact,  it  can  be  shown  that  the  magnitude  of  the  precession  angle  off 
resonance  (o^  «  o»q)  is  1/Qq  times  that  on  resonance  for  the  same  rf 
field  strength.  This  factor  is  related  to  the  phenomenological 

damping  parameter  by  •V2^  • 

The  ratio  of  the  magnitudes  of  %  to  V  gives  a  measure  of  the 
confinement  of  the;  motion  of  the  magnetization  to  the  easy  plane.  This 
ratio  is  called  the  inverse  elllptlcity  of  the  uniform  mode  and  becomes 


-i  ,  _  > 

'V  ’ 


(2.27) 


Using  typical  experimental  values  of  the  parameters  for  ZOgY  in  (2.27)* 
a  value  for  the  inverse  ellipticlty  of  0.044  is  obtained.  This  indicates 
that  the}  motion  of  the  magnetization  is  confined  almost  exclusively  to 
the  easy  plane.  This  is  not  surprising,  considering  the  large  anisotropy 
field  of  Zn^Y  .  To  the  same  approximation  that  the  motion  of  the 
magnetization  is  in  the  easy  plane,  the  value  of  k^  may  be  taken  to 
be  unity.  Thus  from  (2,26a)  it  is  seen  that,  to  first  order,  Q^, 
becomes  Just  the  ratio  of  the  applied  rf  and  dc  magnetic  fields. 
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Since  Otpy  is  by  far  the  dominant  component  of  the  magnetization, 
an  approximate  higher  order  solution  for  the  uniform  precession  is 
obtained  by  letting  cxQz  «  0  in  (2.22).  Thus  for  the  special  case  of 
a  disk  in  the  easy  plane,  the  relation 


V“ox  •  VH0  -  (h/H0)c°s  V 


(2.28a) 


is  obtained  as  a  solution  for  the  uniform  mode.  The  ratio  defined  by 


Ap  *  h/H0  (2.28b) 

will  hereinafter  be  referred  to  as  the  pumping  angle. 

In  order  to  solve  (2.28a)  for  0^  ,  the  ms  of  (2.28a)  may  first 

be  expended  in  a  power  series  of  0^  ,  since  <  l  .  The  method  of 

successive  approximation  demonstrates  that  may  be  expressed  as  a 

Fourier  series  of  odd  harmonic  terms  in’  the  pump  frequency  oip  .  The 
result  of  the  second  iteration  is 


cos  ojpt  +  cos  3ujpt  +  A^  cos  5u>pt  +  . . . , 


(2.29a) 


where 


A?  + 


8  p  128  P 


(2.29b) 
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Since  the  pumping  angle  is  usually  very  small  (Ap  «  l)  ,  then 

A  and  A  are  given  approximately  by  the  first  term  in  the  expansion 

“ 0  ,  “l 
for  each  case. 

E.  DERIVATION  OF  EXPRESSIONS  FOR  THE  INSTABILITY  THffBRHOrjnfi 

The  equations  from  vbich  the  thresholds  of  spin  vave  instabilities 
are  to  be  calculated  are  obtained  from  (2.23)  and  (2.24).  At  this  point 
the  assumption  is  made  that  the  motion  of  the  uniform  precession  of  the 
magnetization  is  confined  to  the  x-y  or  easy  plane  (see  Fig.  2.1).  This 
amounts  to  neglecting  in  comparison  vlth  .  Assuming  also 

that  a  «  1  ,  it  is  then  possible  to  express  h  by  a  power  series 

My  y 

in  [see  (2.28)].  The  equations  of  motion,  greatly  simplified  by 

these  assumptions,  may  be  written  as 


with  coefficients  Sq 
precession  amplitude 


■  a0  ^  +  “l  (2.30«) 

^  -  bO<\z  +  bl%r  '  (2.30b) 

,  a^  ,  , .  b^  expanded  in  powers  of  the  uniform 

OL  as  follows : 

Oy 


*00  +  aQ2  °Qy  +  a04  aQy  +  *  *  * 


~  “  *  a10  +  *11  a0y  *  ®12  °^y  +  ®l4  CQy  *  *  *  *  &-30c) 


^1 


b10  +  bll  aQy  +  b12  °^y  +  bl4  °0y  + 
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The  coefficients  in  the  expansions  (2.30c)  are 

a00 

2  2 

•  -  (x  +  d  +  p  sin  0  sin  jfl) 

aG2 

*  ^04  "  |(-Hjj  +  d+  p+  sin2  ©  sin2  f) 

**10 

1  2 

*  -  2a^2  *  -  80^  *  "  2  sin  ®  sin 

an 

«  ^  sin  2©  sin  $ 

b10 

2  2 

»  x  +  d  +  sin  ©  cos  $ 

bn 

*  -  sin  2©  cos  f> 

b12 

-i  "J  P  1  P  P 

«x  +  ^Njj  +  ^-d  +  cos  ®  “  2  sin  ®  cos  $ 

bl4 

or  12  12  2 

*  x  +  gd  +  gN||  +  ^  cos  ©  -  g  sin  ©  cos  $ 

vhere 

H0  -  Hll  kMs  H 

4xM  4xM 

s  s 

ha 

P  *  - 

4nM 

s 

Dk2 

d  *  ■  -  ■  — ■ 

4xM 

(2.3oa) 


(2.30e) 
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Equations  (2.30)  show  that  the  coefficients  in  the  spin  vave  equation 
of  motion  are  functions  of  time.  Eliminating  <2^  from  (2.30a)  and  (2.30b) 


and  utilizing  the  relation  b. 


®1 


results  in 


(2.31) 


If  vere  eliminated  instead  of  0^  ,  then  the  differentiated  (dtotted) 

coefficients  would  be  different.  However ,  it  will  presently  be  shown  that 
these  dotted  coefficients  may  be  neglected.  In  this  case  the  differential 
equation  has  the  same  form  for  either 
and  (2.10)  yield* 


Equations  (2.9; 


lot  .  *  -lot 

Sk  e  +  a-k  e 


(2.32) 


The  uniform  mode  amplitude  appearing  in  the  coefficients  of  the 

differential  equation  (2.31)  is  given  by  (2.29).  If  only  the  dominant 
first  harmonic  pump  frequency  component  of  0^.  , 


(2.33) 


is  used,  (2.31)  is  a  Matthleu  equation. 

From  (2.30),  (2.31),  and  (2.33)  it  is  seen  that  terms  like  a^  Q^. 
or  b^  aQ  0^  are  proportional  to  ,  while  terms  like 

or  are  proportional  to  inyi^  a^  .  The  latter  terms  may  be 

neglected  since  they  are  smaller  by  the  factor  oup/o^,  «  1  as  well  as 
being  in  quadrature  with  the  former  terms.  The  term  containing  is 

proportional  to  o>  o>p  and  may  also  be  neglected  since  o>  oip  «  . 

This  is  not  a  loss  term,  since  two  time  derivatives  are  involved.  With 
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the  above  negligible  terms  onjitted,  and  with  the  inclusion  of  phenomeno¬ 
logical  dancing,  (2.31)  may  now  be  written  as 


°fcy  +  5  +  CCV 


(2.3*0 


where 


C 


*“l  -  bl  a0 


From  (2.30c)  it  is  possible  to  express  C  as  a  series  in  powers  of 
as  follows: 


C0  +  C1  %  +  °2  +  c3  °0y  +  ••• 


(2-35) 


whc  re 
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fl 
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aio  “  bio  aoo 


2  aio  an  ~  bn  aoo 


‘  2al0  a12  ‘  all  b12  S00  "  b10  a02 
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"to 
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Spin  wave  thresholds  are  now  determined  by  substituting  (2.32)  and  (2.33) 
into  (2.3*0.  Taking  spin  waves  at  the  frequency 


CD 


(2.36) 


and  neglecting  off -resonant  frequency  terms ,  a  secular  relation  for 
and  <2*^  is  obtained,  and  from  this  the  instability  threshold  is  found. 

This  procedure  may  be  followed  for  every  order  n  .  In  doing  this,  only 
that  term  in  (2.35)  which  has  the  sane  degree  as  the  order  of  the  instability 
is  found  to  be  relevant. 

Consider  as  an  example  the  calculation  of  the  first  order  threshold. 

By  ranking  the  above  substitutions  and  keeping  only  constant  and  funda¬ 
mental  .  pump  frequency  terms  in  C  ,  then  (2. 3*0  becomes 


( 


+ 


(2-37) 


and  its  adjoint  (complex  conjugate  with  k  -*•  -k),  subject  to  condition 
(2.36)  for  n  *  1  .  Terms  from  the  uniform  precession  amplitude  up  to 

fourth  degree  are  shown  in  order  to  indicate  the  approximate  nature  of 
the  threshold  calculations.  The  threshold  is  determined  by  setting  the 
determinant  of  the  coefficients  of  the  spin  wave  amplitudes  (®jy  »  a*^) 
in  (2.3(7)  and  its  adjoint  equal  to  zero,  giving 


A 

„  “0  ' 
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(2.38) 


The  resonance  condition  in  (2.38), 


or 


°0  + 


cAj2  + 


(2.39) 


If  the  uniform 


is  seen  to  depend  upon  the  uniform  precession  amplitude. 

precession  amplitude  is  small  (the  usual  practical  situation),  the  KHS 

2 

'0  ”  '' 

quency  from  linear  theory.  Then  the  first  order  threshold  is  found  from 
(2.38)  to  be  approximately 


of  (2.39)  “ay  be  reduced  to  CA  «  ,  vhere  is  the  spin  wave  fre- 


(A  )W 

v  p'cr 


min 


icj 


(2.4o) 


subject  to  the  condition 


0) 


here  and  hereafter  the  order  of  the  threshold  will  be  indicated  by  the 
superscript  on  A 

P 

The  higher  order  thresholds  are  determined  in  a  similar  manner.  The 
thresholds  of  even  order  n  «  2, 4  are  given  by 
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To  find  the  spin  wave  with  the  lcafeAt  threa^o^d^  (2 . 4l )  must  be  minimized  with 
respect  to  k,©,$  subject  to  the  condition  (2.36).  Each  order  threshold 
must  be  considered  separately  in  the  minimisation  process,  since  Cq  as 
well  as  ^  is  different  for  each  case.  In  determining  thresholds  of 
odd  order  higher  than  the  first,  it  is  necessary  to  take  into  account 
the  effect  of  higher  odd  harmonic  terms  in  the  uniform  precession  amplitude 
upon  the  modulation  of  the  spin  wave  frequency  <n^  .  The  even  order 

thresholds  are  unaffected  by  these  odd  harmonic  terms,  however,  because 
they  contribute  only  off -resonant  frequency  terms,  which  are  negligible. 

The  threshold  expression  for  the  third  order  then  becomes 

•(#)*] / 

lc3  -  Cj/21 


subject  to 


o> 


The  higher  odd  order  thresholds  are  similarly  affected  by  the  correspond 
ing  odd  harmonic  terms  in  the  uniform  precession  amplitude  (2.29)  •  How¬ 
ever,  only  thresholds  through  fourth  order  are  being  considered  here. 


F.  OPTIMIZATION  AND  CALCULATION  OF  THE 


I  OLDS 


In  the  optimization  of  spin  wave  thresholds  it  will  be  assumed  as  an 
approximation  that  the  phenomenological  damping  parameter  cd/q  is  inde¬ 
pendent  of  k,Q  and  0  .  Hence,  for  each  threshold  calculation  it  will 

be  assumed,  that  Q  is  a  constant. 


1.  Minimum  Thresholds  at  Spin  Vave  Resonance 


It  is  clear  from  the  threshold  expressions  obtained  in  the  previous 
section  that  the  lowest  thresholds  occur  when  the  condition 


«  ai  (2 .  h3) 

is  satisfied.  This  means  that  the  pumping  of  the  instability  at  o>  occurs 
at  the  spin  wave  resonant  frequency  ,  i.e.,  within  the  spin  wave 

manifold . 

l/2 

The  resonant  frequency  of  the  spin  wave  modes  given  by  ui^  *  C^'  may 
be  expanded  with  the  aid  of  (2.30)  and  (2.35)  to  yield 


[x'(x'  +  p)  +  sin2  Q(x'  +  p  cos2  0)]1/2 


.Recall  that 


(2.44) 


x'  =  x  +  d  , 


and 

x  «  H/4«Ms 
d  .  D  k2/4*Ms 
p  .  haA*m6 

from  (2.31).  When  plotted  as  a  function  of  k  for  a  given  value  of  x 
as  shown  in  Fig.  2.2a,  (2.44)  describes  a  continuum  of  resonant  spin  wave 
frequencies,  the  spin  wave  manifold.  It  may  be  seen  from  (2.44)  that  the 
dependence  of  oj^  upon  ft  is  due  to  the  presence  of  anisotropy.  The 
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*1  *2  X3  X4  x  -  H/Urti. 


FIG.  2.2 — a.  Sketch  of  the  spin  vave  manifold,  shoving  dispersion 

curves  in  terms  of  the  propagation  angles  0  and  . 

b.  Sketch  of  the  k  *  0  intercepts  of  the  spin  vave 
manifold  as  a  function  of  normalized  dc  field. 

The  dashed  lines  in  both  figures  represent  the 
frequency  conditions  (2.3 6), 


2' 


dashed  lines  Indicate  the  required  frequencies  for  instabilities  of  various 
orders.  For  the  case  shown,  a  resonant  instability  of  order  one  is  forbidden 
and  the  k  values  for  higher  order  instabilities  are  found  from  the  inter¬ 
sections  of  the  dashed  lines  with  the  dispersion  curves.  Equation  (2.44) 
shows  that  the  position  of  the  a^anifold  shifts  with  the  dc  field.  This 
dependence  is  shown  in,  Fig.  2.2b,  as  a  function  of  the  normalized  dc  field 
x  ,  for  the  k  «  0  intercepts  of  the  manifold.  In  this  figure  the 
abscisjaa  x  *  x^  represents  the  normalized  field  in  Fig.  2.2a,  where  the 
second  order  resonant  instability  is  Just  on  the  edge  of  the  forbidden 
zone.  As  Fig.  2.2b  shows,  the  number  of  resonant  instabilities  varies 
with  the  dc  field.  The  critical  fields  at  which  the  various  instabilities, 
"come  in"  at  k  *  0  are  found  from  the  intersections  of  the  dashed  and 
solid  curves.  For  the  general  case  of  spin  waves  with  k  *  0  ,  these 

critical  points  are  found,  with  the  aid  of  (2.36)  and  (2.44)  to  be 

n  q  «  x'(x'  +  p)  +  (x'  +  p  cos  0)  sin  0  ,  (2.45a) 


where  the  spin  wave  resonance  condition  (2.43)  has  been  invoked,  and  the 
symbol 


has  been  defined  for  convenience.  Since  the  pump  frequency  is  far  below 
ferromagnetic  resonance,  then  q  «  1 

Since  x'  >  0  and  p  >  0  ,  (2.45)  requires  that 


sin  0  cos  <  nq  p"^^ 


(2.45b) 
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which  places  some  restrictions  on  the  permissible  values  of  Q  and  $  . 

The  RHS  of  (2.45)  is  much  less  than  unity  for  the  parameters  to  be  used. 
Hence,  either  0  is  near  zero  or  else  0  is  near  x/2  .  The  reason 
for  this  is  apparent  from  Fig.  2.2a,  which  applies  to  the  case  q  «  1  , 


p  »  1  . 

Assuming  a  pump  frequency  of  (Op/2*  *1.32  kHc  and  the  parameter 
values  for  Zn^Y  at  29°C  ,  given  in  Appendix  A,  the  maximum  allowable 

values  of  x'  may  be  determined  from  (2.45)  by  setting  0*0  .  Ihe 


Recall  that 


these  maximum  values 

of  x '  are : 

for  n  =  1 

x^  *  0.00250 

n  «  2 

*2  -  0.00998 

n  *  3 

x'  «  0.022k 

n  *  4 

-  0.0396 

✓ 

X 

*  x  +  d  , 

(2.46a) 


where  x  =  H/4itM  is  the  normalized  internal  magnetic  field  and 
d  «  Dk  /4kM  is  the  normalized  exchange  field.  It  is  apparent  from 
Figs.  2.2a  and  2.2b  that,  within  the  spin  wave  manifold,  x'  (and  hence 
the  internal  magnetic  field)  is  small  because  the  pump  frequency  is  near 
the  bottom  of  the  spin  wave  manifold.  For  the  k  *  0  intercepts  of  the 
manifold  (shown  in  Fig.  2.2b)  the  maximum  values  of  x'  *  x^  given  above 
become  maximum  values  of  x  «  x^  .  The  corresponding  values  of  the 
applied  (external)  magnetic  field  Hq  then  become 


H0n  '  <xn  +  NJ|>  4**s 


(2.46b) 
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These  fields  will  be  "passed  through"  during  the  experimental  measurements 

of  instability  thresholds  in  order  to  see  if  any  spin  wave  instability 

effects  are  present  (see  Chapter  Vi).  Thus,  as  x  is  increased  (by 

increasing  the  dc  magnetic  field  HQ)  above  a  certain  xn  ,  the  n 

order  spin  wave  instability  is  forced  off  resonance.  The  x  then,  cor- 

th  n 

responds  to  the  dc  field  at  which  the  n  order  spin  wave  instability 
first  becomes  allowed  on  resonance.  Since  these  maximum  values  of  x' 
are  so  snail,  it  will  be  quite  permissible  later  to  neglect  x  and  d 
in  any  term  containing  p  ,  since  p  >  1 

Each  threshold  will  now  be  considered  separately  and  optimized  with 
respect  to  the  four  parameters  oi,  9  >  f  ,  k  .  The  first  order  threshold 
at  spin  wave  resonance  (u>  =  cn^)  in  terms  of  these  four  parameters  is 
found  from  (2.40)  to  be 

2 

(Vcr  '  "lnSQ-  (7Vp)  cos  f  sin  26  (2'U7) 

subject  to  (2.45)  for  n  *  1  .  To  minimize  (2.47)  as  it  stands  would 
require  $  *  0  and  0  *  jt/4  .  However,  the  frequency  condition  (2.45) 
precludes  this  combination  of  $,©  by  specifying  that  either  ©  or 
*/2  -  j  t  is  small .  Since  the  denominator  of  (2.47)  is  doubly  periodic 
in  ©  ,  only  values  0  <  ©  <  x/4  ought  to  be  considered.  Moreover, 

from  (2.45)  it  is  clear  that  a  smaller  ©  permits  a  large  value  of 
cos  $  ;  or,  conversely,  a  small  value  of  x/2  -  permits  a  large  value 

of  sin  2©  .  Minimization  then  reduces  to  consideration  of  the  following 

two  cases:  (l)  ©  small,  «  0  and  (2)  ©  *  n/4  ,  $  near  jr/2 

The  former  yields  the  lowest  threshold  and  is  given  approximately  by 


(A 

^p'cr 


[q2  -  x'(x'  +  p)]1^2 


(2.48) 
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It  is  evident  from  (2.1*8)  that  the  actual  minimum  threshold  occurs  when 
iL  *  x  +  d  is  smallest,  that  is, when  d  0  .  Therefore,  those  spin 

waves  which  go  unstable  must  be  those  of  lowest  k  m  0  .  This  tends  to 
violate  the  original  plane  wave  assumption  for  1  surge  k  .  The  method 
is  justified  in  practice,  though,  since  measured  thresholds  have  been 

3 

found  to  agree  with  the  theory. 

The  second  order  threshold  at  spin  wave  resonance  is  obtained  from 
(2.4l)  by  proceeding  in  a  similar  manner.  The  result  is 


(A  }(*) 
v  p'cr 


min 


16  cl 


pfd  4-  x-1 


-  sin2©  COS20(p*d+l)  4  (Hjj4x)(m4d4p/2)4l  4  — ^ 


(2.49) 


Since,  by  (2.45b),  sin  ©  cos  is  a  very  small  quantity,  then  clearly  the 
minimum  threshold  occurs  for  ©  «  0  .  Inserting  ©  »  0  into  (2.49)  gives 


(A  )^ 
'Ap'cr 


(2.50) 


where  x  ,  d  terms  have  been  neglected  since  x  4  d  <  0.01  .  For  ©  »  0  , 
it  is  seen  from  (2.45a)  that  d  may  be  large  when  x  «  x^  .  This  means 
that  the  unstable  spin  waves  may  have  large  k  ,  by  contrast  with  the 
first  order  case.  From  (2.50)  the  threshold  is  independent  of  k  .  This 
result  for  the  on-resonance  second  order  threshold  for  spin  waves  is 
compatible  with  the  on-resonance  threshold  result  obtained  by  Schlomann, 
et  al.,1^  for  the  case  of  an  infinitely  thin  disk  (Hj|  »  0) 
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The  third  order  threshold  at  spin  wave  resonance  is  found  from  (2.42) 

to  be 


(V«  ■  Bln 


72  q 


Q  sin  2©  cos  0  ( p  +  d  + 


(2.51) 


The  denominator  has  the  same  form  as  the  first  order  threshold.  Hence, 
the  minimum  occurs  at  0  =  0  ,  ©  small,  and  d  0  --  giving,  with  the 

usual  assumptions. 


(3) 

cr 


36  q4 


Qp1/2  [9<r  -  X 


172 


o 


(x  +  p)] 


(2.52) 


The  fourth  order  threshold  at  spin  wave  resonance  is  determine^  from 
(2.41)  to  be 


,i/ii 


(A  )(lt) 
p  cr 


512  q 


min 


Q(Nii  +  x) 


l  +  ^(p+d)  -  sin2©^"^  +  cos20^ 


(2.53) 


where  second  order  product  terms  in  d,  x,  Njj  have  been  neglected. 
Employing  the  same  type  of  reasoning  as  used  for  the  second  order 
threshold  above,  the  minimum  of  (2.53)  is  found  to  occur  at  0=0  , 

©  »;  0  ,  and  d  large.  The  resulting  fourth  order  threshold  becomes 

approximately 
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(2.5*<) 


2.  Minimum  Thresholds  For  Off -Resonance  Spin  Waves 

Recall  from  the  last  section  that  when  x'  is  increased  above  a 
certain  x^  (by  increasing  the  dc  magnetic  field),  the  n^  order  spin 
wave  instability  is  forced  off  resonance.  In  this  case,  for  a  given  n  , 
the  relationship  between  the  frequency  to  at  which  the  instability  is 
being  pumped  and  the  resonant  spin  wave  frequency  as^  it 


®k 


>  as 


(2.55) 


The  first  order  threshold  for  off -resonance  spin  waves  may  be  found 
from  (2.40)  by  substituting  the  value  for  as  was  done  for  the  on- 

resonance  thresholds.  The  result  is 


(1) 

cr 


min 


r  /  £ 

A2?/2 

2 

.2  2\2  ( as 

(id  +  0^)  +  ^ 

2  ^ 

<^,(x'  +  p)  cos  0  sin  20 


(a.  36) 


subject  to  the  parametric  frequency  condition 


<n 


/ 


where  the  off -resonant  spin  wave  condition  (2.55)  must  be  satisfied.  T3ie 
frequency  condition  (2.55)  places  only  a  lower  bound  on  ;  and  (2.5 6) 

is  to  be  minimized  with  respect  to  k,©,$  ,  subject  to  this  side  condi¬ 

tion.  From  (2.44),  condition  (2.55)  will  be  satisfied  for  all  spin  waves 
if  it  is  satisfied  for  the  0*0  ,  k  *  0  spin  waves.  Therefore  (2.56) 
may  be  minimized  independently  with  respect  to  k,0  ,  and  ^  .  From 

(2.46)  the  condition  on  x  corresponding  to  (2.55)  is 

x  >  0.0025 

for  the  first  order  process  and  the  parameters,  used  before. 
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(2.5 (T) 


Given  that  p  »  x'  ,  the  only  substantial  dependence  of  (2 .56) 
upon  k  occurs  in  the  numerator  through  0^  in  the  resonance  term. 
From  Eq.  (2.44)  or  Fig.  2.2a  it  is  apparent  that  (2.56)  is  a  minimum 
for  k  •*  0  .  To  simplify  the  calculation,  the  assumption  will  now  be 

made  that  oi  is  sufficiently  far  off  resonance  (below  01^)  that 


(2.58) 


Thus  (2.56)  may  be  written  approximately  as 


(A 

v  p'cr 
off  res 


*  ml 


2  +  sin^ 0  ^  ~  +  cos  $  jjj 


cos  <f)  sin  20 


(2.59) 


The  mLnimum  of  (2.59)  occurs  for  $  »  0  and  for  0  small,  with 


0 


opt 


(X  -  q  /p) 


l/2 


(2.60) 


Inserting  (2.60)  into  (2. 59)  yields  the  follow'  ig  minimum  threshold  ex¬ 
pression  for  frequencies  far  off  spin  wave  resonance: 


*  29opt  -  2<x  -  92/p)1/2  •  (2.61) 

off  res 


The  second  order  threshold  off  spin  wave  resonance  may  be  found  from  1 
(2.4l)  by  substituting  the  value  of  ,  as  was  done  in  the  case  of  [ 


the  on -resonance  thresholds  of  second  order.  The  result  obtained  is 


where  the  offuresonance  spin  nave  condition  (2.55 )  waatJb*  satisfied.  Tram 

(2.46)  the  condition  on  x  is  now 

x  >  0.00998  .  (2.63) 


As  in  the  case  of  the  first  order  process,  (2.62)  may  be  Minimized  inde¬ 
pendently  with  respect  to  the  parameters  k,  0,  0  . 

The  numerator  of  (2.62)  is  a  minimum  when  0*0  and  k  *  0  ,  since 

at  these  parameter  values  cjc^  is  a  minimum.  The  denominator  (discussed  in 
the  preceeding  section)  is  a  maximum  at  0*0  and  is  insensitive  to  the 
value  of  k  .  Therefore,  the  minimum  occurs  for  0*0  and  k  **  0  . 

Neglecting  terms  in  the  denominator  which  do  not  contain  p  ,  (2.62) 
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where 


2  2 

(\  «  uijj  (p  +  x)x 


If  the  unstable  spin  wave  is  far  off  resonance,  Eq.  (2.58)  applies,  and 
(2.64)  becomes 


(A  )(S) 

^  p'cr 


x  -  4q 2/p 
1  +  N|,/2 


(2.65)  • 


which  is  independent  of  the  spin  wave  Q  .  It  is  interesting  that  this 
result  has  the  same  form  as  (2.6l).  It  dif feres  in  the  allowed  values  of 
x  ,  which  are  given  by  (2.63),  rather  than  by  (2.58).  Consideration  of 
off -resonance  thresholds  for  processes  of  higher  order  will  be  shown  to 
be  unnecessary. 


It  has  been  assumed  as  an  approximation  in  each  of  the  above  mini¬ 
mizations  that  the  phenomenological  damping  parameter  o/Q  is  the  same 
for  all  spin  waves — i.e.,  the  same  for  all  k,O,0  .  The  relaxation 

frequency  is  related  to  the  phenomenological  damping  parameter  by 


o>/2Q 


(2.6©  ? 


This  parameter  is  usually  defined  in  terms  of  the  spin  wave  linevidth 
as  follows: 


r^Hj/2 


'2.67)1 
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Thus  from  (2 .6|)  and  (2.  «f>  an  expression  for  the  spin  wave  Q  becomes 


Q  *  ai/bWj  •  (2.68) 

For  polycrystalline  yttrium  iron  garnet  with  aluminum  substitutions 
(YA1IG)  ,  the  spin  wave  linewidth  is  found  to  be  proportional  to  fre¬ 
quency  for  ©  =  0  spin  waves. ^  The  spin  wave  linewidth  in  single  crystal 
YIG  at  the  higher  frequencies  is  also  proportional  to  frequency;  but  at 

the  lower  frequencies,  due  to  splitting  processes,  it  is  no  longer  , 

12  13  14  15 

so.  *  *  ’  Each  of  the  thresholds  derived  in  the  previous  sections 

for  Zn^Y  have  minima  for  small  k  ,  and  small  0  or  0  **  0--  since 
instabilities  occurred  at  frequencies  near  the  bottom  of  the  spin  wave 
manifold  (see  Fig.  2.2a).  For  such  spin  waves,  splitting  processes  of 
relaxation  are  not  possible.  Confluence  processes,  however,  are  possible 
but  not  likely  because  k  «  0  .  It  is  therefore  to  be  expected  that  the 

linewidth  of  ©  **  0  spin  waves  in  single  crystal  YIG  will  be  proportional 
to  frequency,  even  at  low  frequencies.  The  spin  wave  linewidth  of  ZngY 
at  1.32  kMc  is  not  known,  but  on  the  basis  of  a 'superficial  comparison 
with  YIG  ,  it  might  be  expected  that  the  spin  wave  linewidth  0?  Zi^Y  - 
for  ©  0  is  also  proportional  to  frequency  at  low  frequencies .  If  this 

is  true,  then  the  spin  wave  Q  of  (2.6$)  will  be  a  constant  independent 
of  frequency. 

Before  the  Q  can  be  determined  it  is  necessary  to  obtain  some 
estimate  for  the  spin  wave  linewidth  of  ZngY  .  Dixon,  et  al.^" 
measured  a  spin  wave  linewidth  of  ^1^  *  1.24  Oe  (for  ©  *  x/2  ,  f  ■  0 
spin  waves)  in  a  parallel  pump  experiment  performed  at  17.2  kMc  for  single 
crystal  samples  of  Zn^Y  with  substituted  manganese.  They  obtained  a 
ferrimagnetic  resonance  linewidth  of  3*8  Oe  at  9  kMc  for  the. .minimum  line- 
width  case.  The  ferrimagnetic  resonance  linewidth  of  the  single  crystal 
ZnpY  sample  used  in  the  experiments  of  Chapter  VI  was  about  35  Oe  at 
6.5  kMc.  Schlomann,  et  al.,  found  for  YIG  that  the  ratio  of  the 
resonance  to  the  spin  wave  linewidths  remains  constant*  Assuming  this 
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relationship  holds  for  Zn^Y  and  that  is  proportional  to  frequency, 

the  spin  wave  linevidth  at  I.32  kMc  is  estimated  to  be 


AHk  »  2.44  Oe 


(2.69) 


Assuming  as  an  approximation  that  the  spin  wave  linevidth  is  independent 
of  ©  ,  then  (2.69)  may  be  assumed  for  all  spin  waves.  The  spin  wave 
Q  may  now  be  calculated  from  (2.68),  with  the  following  result: 

Q  *  193  ,  (2.79) 

independent  of  frequency,  and  propagation  direction.  It  is  also  assumed 
that  Q  is  independent  of  k  . 

4.  Calculation  of  Thresholds 

Prom  the  ferrite  properties  given  by  Appendix  A  and  the  spin  wave  Q 
given  by  (2.7Q),  the  spin  wave  thresholds  uay  now  be  calculated.  The 
actual  parameter  values  employed  are 

P  -  3.22 

q  =  O.O898  , 


and 


N|j  -  0.099 

is  the  demagnetizing  factor  in  the  plane  of  the  ferrite  disk  used  in 
the  experiments  reported  in  Chapter  VI. 

Figure  2.3  shows  curves  of  the  spin  wave  thresholds  through  fourth 
rder  derived  in  this  chapter.  Thresholds  for  both  the  on -Resonance 
(ux.  »  to)  and  off -resonance  (co^  >  to)  cases  are  shown.  The  even  order 
resonant  thresholds  are  nearly  constant  over  the  permissible  range  of  x  , 
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■resonance 


—  off-resonance 
•  -  Off -resonance 
Parameter  values: 


V2*1  *  1-32  Wfc 

Q  *  193 

See  Appendix  A  for. 

°f »  —  , .f*^r 


0.01 


0.02 


ferro^^eJ?feSholds  t}>r°ugh  fourth  ord«  *  *  m  H/U 

3ff-rSa»ce  ”i°n“ce'  Thresholds  ^Iwhfr*llu“ciei'  f«-  ' 
spin  «aves  are  shoanT  b°th  0“-™»onance  a* 


and  it  has  been  shown  that  the  unstable  spin  waves  propagate  along  the 

direction  of  the  dc  magnetic  field  (0i  =  0)  over  this  range  of  x 

The  odd  order  resonant  thresholds  vary  significantly  with  x  due  to  the 

dependence  of  each  threshold  upon  0  .  At  the  bottom  of  the  spin  wave 

manifold,  the  value  of  ©  is  zero  and  the  value  of  x  is  a  maximum  [see  (2.46)]. 

At  x  =  0  ,  ©  is  a  maximum — with  0  *  0.05  for  n  *  1  and 

*  max 

©  =0.15  for  n  =  3  .  All  of  the  resonant  thresholds  show  discon- 

tinuities  at  the  maximum  values  of  x  corresponding  to  unstable  waves 
at  the  edge  of  the  spin  wave  manifold.  Hence,  for  each  order  n  ,  the 
spin  wave  instability  of  order  n  is  no  longer  allowed  on  resonance 
above  a  maximum  value  of  x  [determined  by  (2.46)]. 

For  both  the  even  and  odd  order  instabilities,  the  off -resonance 
thresholds  exist  of  x  greater  than  these  maximumr  values .  However,  for 
values  for  x  slightly  less  than  the  critical  value  xn  ,  the  off -resonance 
threshold  dominates  the  on-resonance  threshold.  The  reason  for  this  is 
that  the  odd  order  on-resonanoe  thresholds  become  infinite  at  the  value  of 
x  corresponding  to  the  edge  of  the  spin  wave  rrpnifold,  and  the  finite 
off -resonance  threshold  takes  over  before  this  point  is  reached  in  each 
case.  This  behavior  is  shown  in  Fig.  2.3  for  the  first  order  threshold. 

The  second  order  off -resonance  threshold  is  also  shown  in  the  figure.  The 
third  and  fourth  order  off-resonance  thresholds  are  similar  in  shape  to  the  , 
plots  of  the  first  and  second  order'  thresholds,  but  have  not  been  shown 
because  they  are  considerably  higher. 

The  spin  wave  thresholds  calculated  for  each  order  n  above,  were 
determined  on  the  assumption  that  the  instability  for  a  given  order  n 
was  the  dominant  one.  However,  it  is  clear  from  Fig.  2.3  that  the 
instability  with  the  lowest  threshold,  regardless  of  its  order,  will 
dominate  the  rest.  Thus  for  the  spin  wave  Q  *  193  *  the  first  order 
threshold  is  dominant  at  the  lowest  values  of  x  .  Where  the  off- 
resonance  first  order  threshold  intersects  the  second  order  threshold, 
the  latter  "takes  over."  The  third  and  fourth  order  thresholds  are  never 
reached,  since  the  second  order  off -resonance  threshold  dominates  them  1 
even  at  the  higher  values  of  x  .  If  the  spin  wave  Q  is  actually 
larger  than  the  estimated  value  of  193#  then  the  thresholds  will  have 
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a  behavior  similar  to  those  shown  in  Fig.  2.3.  This  is  due  to  the  inverse 
dependence  of  the  thresholds  upon  the  n  root  of  Q  >  where  n  is  the 
order  of  the  threshold.  Hence  for  Q’s  larger  than  193>  the  thresholds 
of  lowest  order  will  be  reduced  the  most.  If  the  spin  wave  Q  is 
actually  smaller  than  193>  then  the  lowest  order  thresholds  will  be 
increased  the,  most  due  to  the  change  in  Q  .  In  such  a  situation  it  is 
possible  that  the  second  order  off -resonance  threshold  will  intersect  the 
third  order  resonance  threshold  and  the  latter  will  dominate  for  a  small 
range  of  x  .  It  is  seen  in  Fig.  2.3  that  this  intersection  is  nearly 
possible  with  a  Q  of  193* 

It  is  significant  that  with  such  a  small  range  of  x  it  is  possible 
to  pass  through  so  many  orders  of  spin  wave  thresholds.  This  is  possible 
only  because  the  pump  frequency  is  so  low  (1.32  kMc).  Another  advantage 
of  having  a  low  pump  frequency  in  an  investigation  of  phonon  instabilities 
is  that  spin  wave  thresholds  are  high;  and  it  is  likely  that  phonon 
thresholds  will  be  observed.  An  investigation  of  the  level  of  phonon 
thresholds  is  the  subject  of  the  next  chapter. 

G.  SUMMARY 

In  this  chapter  the  uniform  precession  and  spin  wave  mode  solutions 
to  the  equation  of  motion  of  the  magnetization  were  obtained  by  assuming 
transverse  pumping  of  a  Zn^Y  disk  with  the  eapy  plane  lying  in  the  plane 
of  the  disk.  The  pump  frequency  of  the  rf  magnetic  field  is  assumed  to 
be  far  below  ferromagnetic  resonance.  The  motion  of  the  magnetization 
in  the  uniform  precession  mode  is  shown  to  be  confined  almost  exclusively 
to  the  easy  plane.  Owing  to  the  nonlinearity  of  the  equation  of  motion 
of  the  magnetization,  spin  wave  modes  are  coupledjto  the  uniform  mode 
and  thresholds  of  spin  wave  instabilities  are  determined  up  to  fourth 
order.  These  spin  wave  instabilities  occur  at  multiples  of  one-half 
the  pump  frequency.  The  threshold  expression  for  each  order  instability 
is  minimized  separately  for  on-resonance  and  off -resonance  spin  waves. 
Assuming  the  relaxation  frequency  is  the  same  for  all  spin  waves,  an 
estimate  was  then  obtained  for  the  spin  wave  Q  .  Having  the  spin  wave 
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Q  and  the  parameters  of  the  ferrite  disk,  the  optimum  thresholds  were 
then  calculated  as  a  function  of  the  internal  dc  magnetic  field  of  the 
ferrite.  The  results  show  that  the  dominant  spin  wave  thresholds  are 
the  first  and  second  order  on-  and  off-resonance  spin  wave  thresholds. 

The  first  order  threshold  is  dominant  only  for  very  small  internal  dc 
fields.  The  third  and  fourth  order  thresholds  are  never  reached  for  the 
value  of  Q  assumed. 


CHAPTER  III 


FIRST  ORDER  PHONON  INSTABILITY  THRESHOLDS 

A.  COUPLING  MECHANISMS  FOR  PHONON  INSTABILITIES 

Chapter  II  showed  that  spin  wave  instabilities  result  from  nonlinear 
coupling  of  the  uniform  mode  of  the  magnetization  to  spin  wave  modes 
through  the  free  magnetic  energy  of  the  ferrite.  Similarly,  nonlinear 
coupling  of  the  uniform  mode  to  elastic  wave  modes  takes  place  through 
the  magnetocrystalline  energy.  These  elastic  wave  modes  will  then  like¬ 
wise  experience  instabilities  when  the  uniform  mode  reaches  a  certain 
threshold.  These  elastic  wave  or  phonon  thresholds  will  be  observed 
if  they  are  lower  than  the  spin  wave  thresholds.  Some  of  the  reasons 
why  the  phonon  thresholds  are  expected  to  be  lower  than  spin  wive  thresh¬ 
olds  were  given  in  Chapter  I. 

The  magnetocrystalline  energy  referred  to  above  is  the  free  energy 
of  the  crystal  which  depends  upon  the  magnetization  and  crystal  strain. 
This  energy  is  due  to  the  dipolar  and  pseudo-dipolar  interaction  and  the 
dependence  of  such  interaction  upon  lattice  distortion.  Since  the  strain 
is  always  many  orders  of  magnitude  less  than  unity,  the  magnetocrystalline 
energy  E  is  given  accurately  by  a  Taylor  expansion  with  respect  to 
strain,  as  follows: 
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'let* 


(3-1) 


where  the  strains  e 


ij 


represent 


to  the  unmagnetized  state  and  a 


Each  of  the  subscripts  i,  i 
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distortions  of  the  crystal  with  respect 
is  the  normalized  magnetization  vector, 
-t  ,  etc.,  takes  on  the  coordinate 
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values  x  ,  y  ,  z  .  Since  the  magnetocrystalline  energy  must  be  invariant 
to  spatial  inversion  of  the  magnetization,  an  expansion  of  the  coefficients 
of  (3*1)  must  contain  only  even  powers  of  the  magnetization  direction 
cosines.  The  development  of  the  anisotropy  energy  term  E  in  this  way 
leads  to  the  anisotropy  constants --independent  of  the  crystal  strains. 

The  anisotropy  energy  for  a  hexagonal  crystal  was  given  by  (2.l4).  The 
second  term  in  (3*1)  is  called  the  magnetostrictlve  energy  term.  The 
coefficients  of  the  lowest  order  terms  in  <2  for  this  term  are  called 
the  magnetoelastic  coupling  constants.  The  third  term  in  (3*1)  is  some¬ 
times  referred  to  as  the  intrinsic  energy  term,  and  the  lowest  order 

coefficients  in  the  expansion  of  this  term  are  called  the  magnetoelastic 

7 

stiffness  constants. 

Regardless  of  the  magnetic  state  of  the  crystal,  the  elastic  energy 
Ee  is  given  by  a  Taylor  expansion  with  respect  to  strain,  as  follows: 


E 

e 


^Cijk-L  €ij  €ki  +  cijk£ran  €t$  €$£  €mn  +  *''  ^  * 

(3.2) 


where  c^^  8X6  elastic  stiffness  constants  and  c 
the  coefficients  of  the  lowest  order  enharmonic  elastic  energy  term. 

Consider  now  that  the  crystal  is  magnetized  due  to  the  application 
of  a  sat  mating  dc  magnetic  field.  This  causes  linear  saturation 
magnetostriction,  and  the  resulting  equilibrium  strains  which  will  be 
present  are  determined  by  the  relations 


V. d-  (E  +  E  )  «  0  ,  (3«3) 

d€ij  k  e 

thereby  minimizing  the  total  strain -dependent  part  of  the  crystal  energy. 
These  equilibrium  (static)  strains  establish  an  "operating  point"  about 
which  dynamic  strain  deviations  can  take  place  due  to  applied  dynamic 
stresses  and  elastic  wave  oscillations.  For  the  following,  dynamic 
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strain!  and  strain  waves  are  of  particular  interest,  since  they  are 
capable  of  giving  rise  to  instabilities  under  certain  conditions.  It 
will  be  useful  then  to  resolve  each  strain  component  into  two  parts: 
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(3.4) 


where  represent  the  static  or  equilibrium  strains  determined 

by  (3.3 )>  and  the  new  e  .  represent  the  dynamic  strains  of  the  magnetized 
crystal. 

The  strain  given  by  (3*4)  will  now  be  inserted  into  the  expression 
(3.1)  for  the  magnetocrystalline  energy.  The  anisotropy  term  is  unaf¬ 
fected  by  the  strain.  A  typical  example  of  a  smgnetostrictive  energy 
term  would  be 


b  (Of)  +  Cjj  )  ,  (3.5) 

where  b  is  a  magnetoelastic  constant.  Since  the  energy  must  still 
contain  only  even  powers  of  a  ,  then  (a)  must  also  contain 

only  even  powers  of  a  ,  besides  being  very  small.  Thus  the  first  term  of 
(3.3)  contributes  only  a  small  higher  order  correction  to  the  anisotropy 
term.  Ihe  second  term  of  (3-5 )  has  the  same  form  as  the  original  msgneto- 
strictive  energy  term.  It  is  this  term  which  is  responsible  for  what 
is  called  the  magnetoelaatic  effect.  Similarly,  a  typical  example  of 
an  intrinsic  energy  term  would  be 


a  “m“n  <e«  +  £y '(2}  +  <£f)  £«  +  ‘kt5  ' 

(5-6) 

for  the  reason 


.(o). 


(9) 


given  above,  the  must  contain  only  even  powers  of  a  ,  besides 


where  d  is  the  magnetoelastic  stiffness  constant 

(0) 

'ij 

being  very  small.  Thus  the  first  term  in  (3*6 )  contributes  a 
higher  order  correction  term  to  the  anisotropy  energy*  The  second  mad 
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third  terms  contribute  a  higher  order  correction  to  the  magnetos trictive 
energy.  The  last  term  in  (3*6)  is  the  only  one  which  has  the  form  of 
the  intrinsic  energy — and  this  fora  is  the  same  as  appeared  in  the  original 
energy  expression  (3-1).  It  is  this  last  term  which  is  responsible  for 
what  is  called  the  intrinsic  effect.  It  is  apparent  from  (3*6)  that  this 
effect  does  not  depend  upon  strain  produced  due  to  the  orientation  of 
the  magnetization.  From  the  above  discussion  it  is  therefore  concluded 
that  the  magnetocrystalline  energy  given  by  (3*1)  may  be  accurately  used 
in  its  present  form  to  express  the  energy  of  the  dynamic  strains. 

The  strain  (3-M  will  next  be  inserted  into  the  expression  (3*2) 
for  the  elastic  energy.  A  typical  example  of  the  first  term  in  (3-2) 
would  be 
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Using  the  previous  argument  regarding  e 
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(5-7) 


...  (a),  it  is  seen  that  the 

ij  ~ 

first  term  in  (3-7)  is  a  higher  order  correction  to  the  anisotropy 
energy.  The  second  and  third  terms  are  higher  order  contributions  to 
the  magnetostrictive  energy.  Hie  last  term  in  (3.7)  represents  the 
elastic  energy  as  given  in  the  original  expression  (3*2).  A  typical 
example  of  the  anharmonic  energy  term  in  (3*2)  would  be 
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The  last  term  in  (3*8),  which  is  an  anharmonic  term  in  the  dynamic 
strains,  is  neglected.  The  first  and  second  terms  are  higher  order 
corrections  to  the  anisotropy  and  magnetostrictive  energies,  respectively. 
The  third  term  has  the  same  form  as  the  dynamic  intrinsic  energy  term 
in  (3«6)  and  is  responsible  for  what  is  called  the  morphic  effect. 


Even  though  the  energy  of  the  intrinsic  and  norphic  effects  have  the 
same  form,  there  is  a  possibility  of  distinguishing  between  these  two 
effects.  The  energy  of  the  morphic  effect  [third  tern  in  (5*8)],  unlike 
that  of  the  intrinsic  effec  c,  depends  upon  a  distortion  of  the  lattice 
in  response  to  the  direction  of  the  magnetization. 

Therefore,  if  the  magnetization  is  changing  direction  with  time, 
it  is  likely  that  the  inertia:  of  the  crystal  will  be  a  factor  in  de¬ 
termining  the  magnitude  of  the  morphic  effect.  It  has  been  suggested 
that  the  morphic  effect  should  probably  be  the  most  effective  at  tens 

of  megacycles  where  crystals  of  a  reasonable  size  exhibit  well-defined 
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elastic  resonances.  Hence,  at  microwave  frequencies,  it  is  likely  that 
the  morphic  effect  will  be  negligible  in  comparison  with  the  intrinsic 
effect. 

From  the  above  discussion  it  has  been  shown  that  the  magnetocrystalline 
energy  expression  (3.1)  and  the  elastic  energy  expression  (3*2)  may  be 
used  in  their  original  form  to  express  the  energy  contained  in  the  dynamic 
strains  when  a  saturating  magnetic  field  is  present. 

From  the  energy  expressions  above  it  is  apparent  that  the  energy 
associated  with  the  morphic  and  intrinsic  effects  is  a  factor  of  the 
order  of  the  elastic  strain  smaller  than  the  energy  associated  with  the 
magnetoelastic  effect.  It  is  not  surprising,  then,  that  experiment*  in 
single  crystal  nickel  show  the  dependence  of  the  propagation  velocity 
of  acoustic  waves  to  be  well  explained  on  the  basis  of  the  magnetoelastic 
effect  alone.  Even  for  large  bias  fields,  the  intrinsic  effect  was 
negligible  in  comparison  with  the  magnetoelastic  effect.  Therefore, 
the  magnetoelastic  effect  may  reasonably  be  assumed  to  be  the  dominant 
mechanism  in  the  investigation  of  phonon  instabilities  in  7* y,Y  .  It 
will  be  found  later  from  the  experiments  described  in  Chapter  VI  that 
this  assumption  appears  to  be  substantiated. 

Physically,  the  magneto elastic  effect  described  above  causes  a 

reactive  loading  and  damping  of  the  elastic  waves  by  motion  of  the  mag- 
19 

netization.  It  is  further  to  be  shown  in  chapter  that,  when  the 
magnetization  is  pumped  by  the  application  of  a  sufficiently  large  rf 
magnetic  field  (in  the  presence  of  a  dc  magnetic  field),  these  elastic 
waves  exhibit  instabilities  when  the  uniform  precession  amplitude 
exceeds  a  threshold  value. 
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B.  EIASTIC  AMD  MAGNETOELASTIC  ENERGY  CONTRIBUTIONS  TO  EQUATIONS  OF  MOTION 


In  Chapter  II  it  was  indicated  why  Zn^Y  may  be  be  regarded  as 
a  uniaxial  crystal  to  a  first  approximation.  It  has  been  shown  from 

symmetry  arguments  that  the  elastic  energy  for  a  uniaxial  crystal  is 
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where  the  subscripts  of  the  elastic  constants  refer  to  the  transformed 

variables  J  *  x  +  iy  ,  t)  =  x  -  iy  ,  z  . 

Similarly,  it  has  been  shown  that  the  magnetoelastic  energy  for  a 
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uniaxial  crystal  is  given  by 
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where  the  transformed  variable  subscripts  have  been  suppressed  to  simplify 
the  expression  and  to  enable  it  to  be  written  in  an  analogous  form  to  (3*9)* 
Since  the  magnetoelastic  energy  contains  both  strain  and  direction 
cosines  of  the  magnetization,  it  will  contribute  both  to  the  equation 
of  motion  of  the  magnetization  (2.6)  as  well  as  to  the  elastic  equation 
of  motion.  It  is  through  the  magnetoelastic  term  that  energy  is 
transferred  from  spin  waves  to  elastic  waves. 


1.  Spin  Wave  Equations 

The  contribution  of  the  magnetoelastic  energy  to  the  equations  of 
motion  (2.6),  and  ultimately  to  (2.23),  is  determined  by  inserting  (3.10) 
into  the  relation  (2.7)  to  give  the  effective  magnetic  field,  and  then 
inserting  the  latter  into  (2.6).  The  following  results  are  obtained: 
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M  ;  x  y'  yy  xx7  ;5V  x  y  xy 
s 

+  Bi.ql (at  €  -  a  e  )] 

4z'xyz  y  xz 


(3-llc) 


2.  Elastic  Wave  Equations 

The  fundamental  elastic  wave  equation  of  motion  is  given  by 


po*i  *  I  ?<1+ V 


d2S 


(3-12) 
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where  i  ,  j  *  1  ,  2  ,  5  ;  is  the  elastic  displacement  in  the  di¬ 

rection  of  the  Cartesian  coordinate  x.  ;  p  is  the  density;  and  E 

1  c  ^ 

is  the  free  energy  of  the  crystal.  The  factor  ^-(l  +  5^  )  is  needed 
because  symmetrical  physical  strain  terms  in  the  free  energy  have  been 
®mbined.  The  strains  combined  in  this  way  are  referred  to  here  as 
pseudo-physical  strains  since  the  trans formation  properties  of  the  physical 
strain  tensor  are  still  retained;  no  distinction  is  made  between  equivalent 
strain  components.  Each  component  of  the  strain  written  in  tensor  form 
is  given  by 


2  V&Xj  dxt J 


(5.13) 


The  only  effective  contributions  to  the  free  energy  E  in  (3 . 12)  will 
be  due  to  elastic  and  magnetoelastic  energies,  since  other  terms  of  the 
free  energy  are  independent  of  the  strain.  The  expanded  form  of  the 
elastic  equations  of  motion  are  obtained  by  the  substitution  of  (5-9) 
and  (3-10)  into  (3*12)  with  the  subsequent  utilization  of  (3.13)*  The 
result  is  written  symbolically  as  the  sum  of  the  elastic  and  magneto- 
elastic  contributions 


P0Ri  "  (p0Ri^e+  (pORi\n  ' 


(3.14a) 


with  the  elastic  terms  given  by 


We  *  2(2C‘ 


a2R 

+  x 


d2R 


a2R 


+  2C 


+  2C 


*n*n’  ''ssnn'  "  einn  ^y2  £znz 


dzc 


+  4c, 


a2R 


32r 


..  .  _ -  +  2(C.  +  Ct  )  - - 

dxdy  dxdz 


( 3 • l4b ) 
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and 


2 

o  R 

(p nR  )  =  2(2 C  „  +  C  )  - J+2C 

oy'e  m  en  linn  ^y2  « 


32r 


l  inn 


-r^+ sci 


a2R 


+  4c 


in  in 


a2R 


3yd  x 


+  2(c  +  c  ) 

£  r\zz  |  zrjz 


3x 
32r 

z 

3yd  z 


32R  32R 

(p  ri R  )  =  2  ( C  +  C  )  (  -  -~X"  +  — — i 

0  z  e  ft ZTJZ  ^zz  V  3x3z  3y35 


32R  32r 


+  2C. 


z  +  - -  I  +  c 


ft  ZT>2  V  3x2  ’  3y2 


►) 


a2B 


|ZT)2  dz2 


zzzz  ^2  * 


(3.14c) 


(J.UkI) 


and  the  magnetoelastic  terms  by 


(p_R  )  =  r-  [B0(oP  +  C?)  +  Bx(a?  -  0?)] 

v  0  x7m  3x  2fc  x  y  3  x  y7 


+  r^(tV»V+^(Wz> 


(p.R  )  =  [B_ (oP  +  a?)  +  B,  (C?  -  Q?)] 

0  y  m  3y  2V  x  y  3  y  x 


(3.l4e) 


+  I  If  (^V  +  I  57  (Vy“z> 


(3-l4f) 


(p  r  )  =  i  I-  (b.  a  a  )  +  \  (B.  a  a  ) 
v  O  z'm  2  5x  v  4  x  z'  2  d y  4  y  z7 


+  li  (BiQ^) 


(3.l4g) 


Since  (3*14)  is  invariant  under  a  transformation  to  any  rotated  axes  x',y'  , 
the  magnetoelastic  behavior  is  isotropic  in  the  x-y  plane.  Outside 
the  easy  plane  the  behavior  is  anisotropic. 


4 
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As  in  the  case  of  the  spin  wave  modes,  the  assumption  of  periodic 
■boundary  conditions  is  valid  for  acoustic  wave  propagation  because 
2n/k  «  sample  dimensions.  Hence  a  plane  wave  solution  for  the  dis¬ 
placement  R  of  (3.8)  is  in  order  and  is  given  by  the  Fourier  series 
expansion 

Zik  •  r 

Rke“~  .  (3.15) 

k 


The  displacement  R  contains  the  time  dependence  of  the  wave,  and 
has  the  form 


lost 

e 


+  r' 
~k 


-ioit 

e 


(->.16) 


'which  in  general  may  be  complex,  where  the  quantity  o>  is  the  frequency 
of  the  elastic  waves  and  the  vectors  r^  ,  r^  are  independent  of  time. 
Since  the  elastic  displacement  R  must  be  real,  then  the  elastic  wave 
amplitudes  are  related  by 

S-k  -  £k  ’  (5.17a) 


and  therefore 


£k  -  £*-k  •  (5.17b) 

In  general,  R  will  contain  pairs  of  oppositely  directed  traveling  waves, 
with  propagation  vectors  k  and  -  k  .  As  in  the  spin  wave  case,  the 
uniform  precession  will  be  found  to  couple  pairs  of  elastic  waves  r^ 
and  r*  ,  and  this  coupling  will  lead  to  exponential  growth  of  pairs 
of  elastic  waves.  For  each  pair  it  will  be  found  that 

I rk I  *  lr_J  '  (3-18) 

so  that  the  unstable  modes  are  standing  waves. 
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3*  Solutions  to  the  Elastic  Wave  Equation 


An  examination  of  (3*1*0  reveals  that  the  simplest  solutions  of 
the  form  of  a  single  term  in  (3*15)  occur  when  the  propagation  vector 
k  is  along  one  of  the  crystal  axes  x  ,  y  or  z  .  Since  the  magnetoelastic 
equations  are  isotropic  in  the  x-y  or  easy  plane,  it  may  he  arbitrarily 
assumed  that  the  propagation  vector  k  is  along  the  x-axis.  The  obvious 
independent  solutions  to  (3* 1*0  for  elastic  waves  propagating  in  the 
x-y  plane  are  therefore  of  the  form 


Ry  (t,  X) 

Rx  (t,  x),  (3.19a) 

and 

H  (t,  x) 

2  • 


For  elastic  waves  propagating  along  the  z-axis  the  obvious  independent 
solutions  to  (3* 1*0  have  the  form 


Rx  (t,  z) 

By  (t,  z)  (3-19b) 

and 

R  (t,  z)  . 
z  7 


Because  of  the  isotropy  of  the  x,y  plane  the  first  two  are  equivalent, 

and  the  second  solution  will  not  be  considered  further.  The  solution 

to  (3*1*0  with  wave  propagation  outside  of  the  easy  plane  and  not  along 

the  z-axis  is  much  more  complicated  than  are  the  above  solutions.  Such 

21 

waves  are  called  quasi-longitudinal  and  quasi -transverse.  The  name 

quasi  is  applied  to  these  waves  because,  in  general,  the  displacement 
vector  R^  is  neither  parallel  nor  normal  to  the  propagation  vector  k  . 
To  avoid  unnecessary  complications,  the  simple  solutions  to  (3*1*0  will 
first  he  considered  and  phonon  instability  thresholds  will  be  calculated 
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from  these.  It  might  be  expected,  due  to  the  large  anisotropy  field, 
that  parametric  excitation  of  elastic  waves  propagating  or  polarized 
outside  of  the  easy  plane  will  je  greatly  suppressed.  If  the  thresholds 
of  such  waves  having  the  simple  forms  given  above  are,  indeed,  substantially 
higher  than  the  thresholds  of  waves  propagating  and  polarized  in  the 
easy  plane,  then  it  will  be  assumed  that  the  thresholds  of  the  quasi¬ 
longitudinal  and  quasi -transverse  waves  mentioned  above  are  also  sub¬ 
stantially  higher. 

C.  EQUATIONS  OF  MOTION  IN  A  ROTATED  COORDINATE  SYSTEM 

The  direction  of  propagation  of  acoustic  waves  propagating  in  the 
easy  plane  will  not  necessarily  be  in  the  same  direction  as  the  applied 
dc  magnetic  field.  Therefore  it  becomes  necessary  to  consider  acoustic 
wave  propagation  at  an  arbitrary  angle  with  respect  to  the  applied  dc 
field.  Hence,  rotation  of  the  x  ,  y  axes  through  an  arbitrary  constant -angle 
relative  to  the  field  is  required  (see  Fig.  3*1)*  Any  vector 
component  Ai  in  the  original  system  is  transformed  to  the  rotates 
system  A^ ,  .  Thus,  the  transformation  become  ' 

A  ,  =  A  cos  +  A,  sin  r 
x  x  T0  y  T0 

V“  -  Ax  sin  *0  +  Ay  c°s  *0 

A  .  -  A  ,  (3.20) 

Zi  Lt 

where  A^  can  stand  for  0^  ,  k^  or,  in  particular,  any  component  of 
the  right  hand  side  of  the  equation  of  motion  for  the  magnetization  (2.6) 
(see  Fig.  3.1).  In  making  the  transformation  it  is  desirable  to  retain 
the  identity,  so  to  speak,  of  the  displacement  components  R^  in  the 
fixed  system,  so  it  will  be  possible  to  identify  the  simple  modes.  Hiis 
will  be  done  explicitly  in  the  individual  cases.  Where  there  are  only 
terms  in  a  ,  as  in  the  spin  wave  equations  of  motion  before  the 


-  - 


X 


FIG.  3.1 — Orientation  of  rotated  coordinate  system  with  respect  to 

fixed  system,  showing  applied  magnetic  fields  in  the  rotated 
system. 
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magnetoelastic  effect  was  considered,  every  term  in  the  new  system  has 
the  sane  form  as  in  the  old  and  may  simply  he  primed.  The  magnetoelastic 
terns,  containing  components  of  R  as  well  as  a  ,  are  transformed  to 
the  rotated  system  by  (3.20).  Then  a  is  referred  to  the  rotated 
system  by  using  the  inverse  of  (3.20 ).  The  displacement  JR  ,  on  the 
other  hand,  is  left  in  the  fixed  system  in  order  to  retain  the  identity 
of  the  elastic  wave.  Adding  the  magnetoelastic  terms  then  yields 

\r'  *  aoakz'  +  +  (W  (5-2U) 

'  b0\z'  +  +  >  (3*21b ) 


where  the  primes  on  the  coefficients  mean  that  uniform  precession  am¬ 
plitudes  contained  therein  are  also  primed. 


D.  DERIVATION  OF  FIRST  ORDER  THRESHOLDS 


Only  the  first  order  pumping  of  phonon  instabilities  will  be  of 

interest,  since  the  first  order  thresholds  are  lowest.  Taking  only 

first  order  terms  in  QL  /  and  neglecting  Ot  ,  «  a  ,  as  before, 

cy  z  y 

the  coefficients  given  in  (3*21)  become 


✓ 


j 

“k 


* 


(3.22) 
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where 


Each  type  of  elastic  wave  propagation  will  now  be  considered  in  the 
order  shown  in  (3  *19).  For  waves  propagating  in  the  x-y  plane,  it  is 
assumed  that  k*r  «  kx  in  (3.15)*  The  first  such  wave  R  (t,x)  is  for 
a  wave  polarized,  or  having  particle  displacement,  in  the  x-y  or  easy 
plane.  This  case  will  be  treated  in  detail  and  will  serve  as  a  basis 
for  the  other  cases.  The  relevant  components  of  (3*ll)»  which  contain 
only  terns  in  ,  become 


<y« 


a  a 
x  z 


a 

z 


x  y 


(5.23) 


Here  the  strain  is  related  to  the  elastic  displacement  by 


xy 


Ik  p  itat 

r  V* 


(5.2*) 
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which  follows  from  (3.13)  and  (3*15).  Each  component  (£*  )  of  (3.23) 

1  m 

must  be  transformed  into  the  rotated  coordinate  system  in  accordance 
with  (3.20),  and  each  direction  cosine  must  be  transformed  by  the 
inverse  of  (3*20).  Since  it  is  desired  to  retain  the  identity  of  the 
elastic  mode,  the  components  of  the  displacement  £  are  referred  to 
the  fixed  coordinate  system.  Recalling  from  (2.8)  that  each  direction 
cosine  (and  hence  has  uniform  precession  and  spin  wave  com¬ 
ponents,  and  retaining  the  spatially -dependent  terms,  ,  only  to 

first  order  gives 


*  <V4  -  0 

2iriL 

*  \  _  ■  .  J 

«V»  M  <eo*  2*0  -  2<V  sin  s*0)  k  V  e  •  (5.25) 


Substitution  of  (3.25)  into  (3»2l)  and  elimination  of  Qt^  /  yields  the 
following  equation: 


(a,  ,  +  a>  a  ,  *  K  R. 

v  ky  m  xy  mp  Tsy 


(3.26a) 


where 


4  *  4  +  P)  (x'  +  sin2  *0  *  sin  2  *0  aQy'^  (5.26b) 


K 


mp 


2iykB,ui 

- - -Z2.J.  (x'  +  p)  (cos  2 if  - 


M 


2  sin  2*0O^,) 


.  (3.26c) 


The  frequency  cjd^  is  simply  the  resonant  spin  wave  frequency  [see  also 
(2.44)],  and  K  is  the  coupling  coefficient  between  the  spin  waves 
and  the  elastic  waves  arising  from  the  magnetoelastic  coupling.  It  has 
been  assumed  in  deriving  (3.26)  that  the  magnitude  and  direction  of  the 
propagation  vector  k  is  the  same  for  elastic  waves  as  for  spin  waves. 
This  is  reasonable  since  maximum  mutual  coupling  between  the  spin 
and  elastic  waves  requires  the  magnitude  and  direction  of  k  to  be  the 
same. 
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The  elastic  wave  equation  for  the  case  under  consideration  Is  the 
one  for  which  1  »  y  in  (3.12).  Utilizing  the  transformation  (3*13) 
for  the  direction  cosines  and  retaining  spatially  dependent  terms  only 
to  first  order ,  the  aagn etoe las t 1 c  contribution  (3.l4f)  becomes 

(PoVm  *  2iB5  (cos  ^0  ’  ,ln  ** 0 )Lk  *V  eifaC  '  (5*27) 

k 

Adding  (3*27)  to  the  strictly  elastic  tern  (3-l^c )  gives  the  equation 


+  a> 


(3.28a) 


where 


2ik  B, 

. 3  (.cob  2(r 

P0 


20^,  tin  2*0) 


(3.28b) 


(3-28c) 


where  again  the  spin  wave  and  phonon 
The  equations  of  notion  for  the 
are  (3*26)  and  (3*28).  Losses  are 
these  equations,  with  the  result 


k's  must  be  equal. 

combined  magnetic -elastic  system 

introduced  phenomenologically  into 


v +  2\  V +  4  V  ■  ■wv 

(3.29») 

V  +  2"pV  +  <0pV  *  s(- V  » 

(3-29b) 

where  the  relaxation  frequencies  tj^  -  a^/2Q 
defined  in  terms  of  the  respective  resonance 


>  1,  *  v*S>  *” 

frequencies  and  Q  factors. 


-  59  - 


A  physical  interpretation  of  (3*29)  is  found  by  assuming  zero 
pumping  field  (i.e.,  0^*0).  First  insert  (2.9)  primed  and  (3.16)  into 
(3.29).  A  nontrivial  solution  exists  only  vhen  the  secular  determinant 
is  sex  equal  to  zero.  If  the  frequency  is  assumed  to  be  complex. 


-»  u)  +  iq 


(3.30) 


the  solution  may  be  expressed  as 


(oP  -  aP  )(aP  -  aP )  =  K  K 

'  p  m  pm  mp 


(3.31a) 


,(<4  -  > +  y  4 

222 

u)  +  a)  -  2a> 
m  p 


(3.31b) 


where  it  has  been  assumed  that  q  «  o>  .  Note  the  symmetry  of  the 

solutions  with  respect  to  a>  ,  o>  .If  the  roots  of  (3 -31a)  are 

®  P 

plotted  as  a  function  of  k  ,  as  in  Pig.  3*2,  the  usual  dispersion  curve 
is  obtained.  The  interaction  is  greatest  at  the  cross-over  frequency, 
where  o>  *  o>  .In  the  system  considered  here,  however,  it  will  be 


required  that  <b  «  o> 

JL 

proxiaated  as 


Under  these  conditions,  (3*31)  can  be  ap- 


2 

<0  - 
=  P 


(5-J2a) 


V  \  J!*riSL 

*  /i\ 


(3.32b) 


This  applies  to  the  lower  part  of  the  lower  branch  in  Fig.  3*2,  where 
the  mode  is  substantially  an  elastic  wave.  The  phonon-like  nature  of 


the  resonance  (3* 32a)  is  denoted  by  the  subscript  p  on  00 


(3.32a)  indicates  the  dependence  of  the  phonon  frequency 


PO  * 


Equation 
upon  the 
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spin  wave  frequency  uim  ,  as  well  as  upon  the  spin  wave -phonon  coupling 


parameters  K _  ,  K 

BJP  pPm 

that  K  K _ «  of  . 

pm  lip  p 


For  typical  values  of  the  parameters  it  is  found 

is 


2  2 

Then,  since  o>m  »  oj^  ,  the  frequency  oi  ^ 


nearly  equal  to  the  unperturbed  phonon  frequency  a> 


Also,  the  re¬ 


laxation  frequency  tj  is  nearly  equal  to  the  phonon  value  r4  ,  and 
is  relatively  insensitive  to  the  value  of  . 

Equation  (3- 31)*  and  hence  (3*32),  were  derived  under  the  assumption 
of  zero  pumping  field  (i.e.,  aQy'  *  0)  •  In  determining  the  phonon  in¬ 
stability  threshold,  the  modulation  of  cu  and  q  through  the  time- 
varying  quantity  ®Qy'  be  considered.  The  threshold  is  determined 

most  simply  by  obtaining  the  differential  equation  for  alone.  First, 

writing  (3*29a)  in  operator  form  results  in 


(D2  +  2V>  +  4)  V 


(3-33) 


where  D  symbolizes  the  partial  time  derivative  operator.  Assuming 
operation  fair  below  the  spin  wave  manifold,  the  first  term  in  (3. 33) 


may  be  neglected  by  comparison  with  (x> 


m 


and 


°ky 


✓  may  be  expressed  as 


,-l 


mi 


2q 


1  + 


m 


m 


V 


(3-3*0 


m 


Expanding  the  operator  expression  and  taking  only  the  first  two  terms 
(assuming  2\/<um  *  <<:  l)  *  and  substituting  (3.3*0  into  (3*29) 

gives  the  result 

V  +  &1V  +  “p2V  *  0  ’  <3-35) 

where  R.  has  the  form  (3*16)  and  the  values  of  q  and  a/  are  of 
x  p 

the  same  form  as  (3*32)  for  the  case  0^,  »  0  .  However,  in  the  derivation 

of  (3*35)  <i>m  ,  K  and  are  functions  of  the  pump  or  uniform  precession 

amplitude  * 
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It  follows  that  the  phonon  frequency  will  be  nodulated  at  the  pump 
frequency,  and  to  first  order  in  aQy'  v*11  give 


o> 


*2 


%0  <1  +  G 


(5.36) 


where  from  (2.29),  a  ,  =  A  cos  out  to  the  first  approximation.  Hie 

Qy  P  r 


Oy  p 

modulation  parameter  Q  is  given  by 
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1 

a  4?  . 
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K  K 
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ba  ; 
oy 

“  ~2 
0) 
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d<V 

J 

.  (3-37) 


The  damping  parameter  q  will  also  be  modulated  by  the  pump  amplitude 

a_  ,  through  the  dependence  of  q  upon  a)  and  the  coupling 
ty  m 

coefficients.  However,  q  is  nearly  equal  to  the  phonon  damping  parameter 
,  as  was  previously  discussed,  and  if  the  phonon  is  fairly 

large,  then  q^  «  .  Therefore,  the  effect  of  modulation  of  the 

damping  parameter  q  upon  the  instability  threshold  will  be  negligible, 
as  can  be  seen  by  comparing  the  contributions  of  the  second  and  third 
terms  in  (3*35)  in  the  following  procedure. 

Substituting  (3*16)  into  (3*35  )>  and  assuming  the  frequency  condition 


u> 


ay's 


(3.38) 


for  a  first  order  instability  yields  the  equation 


p  + 


io>  0) 


S 


2 

“to 


■)r»* 


GA 

2 


2  # 
*pO  r-ky 


(3.39) 


and  its  adjoint  (cpmplex  conjugate  with  k  -»  -k)  ,  where  nonresonant 
frequency  terms  have  been  neglected.  ^  ^  Since  a  nontrivial  solution  to 


^  ^Chapter  IV  gives  a  more  complete  treatment  of  the  solution  to 
(3*35 ),  including  operation  above  threshold  where  the  solution  is  unstable. 
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(3*38)  exists  only  when  the  secular  determinant  is  equated  to  zero, 
then  the  minimum  pumping  threshold  becomes 


cr 


min 


^  2  ,2 
+  V 


+ 


1/2 

.  (3.*0) 


It  is  appa'  it  from  (3*40)  that  the  minimum  threshold  will  obtain  when 
the  resonance  condition 


oj  .  =  01 

pO 


(3-41) 


is  satisfied.  Since  the  second  tern  in  (3* 52a)  is  small. 


05 


PO 


is  given 


approximately  by  oi^  of  (3.28).  The  minimum  phonon  threshold  on 
resonance  then  becomes 


(Vcr  _  ’  (3,42) 

subject,  of  course,  to  (3.38)  and  (3*4l).  Conditions  (3*38)  and  (3.4l) 
determine  the  value  of  the  optimum  propagation  constant  k  at  which 
the  minimum  threshold  occurs.  The  threshold  (3.42  )  must  then  be  further 
minimized  relative  to  the  angle  ,  which  appears  as  a  variable  in  the 
modulation  parameter  G  .  For  the  spin  wave  threshold  calculations,  a 
similar  condition  to  (3.4l)  included  a  dependence  upon  the  direction  of 
propagation  of  the  waves,  and  the  value  of  k  was  not  specified  by 
the  frequency  condition  alone.  However,  in  the  phonon  case  k  is  de¬ 
termined  by  (3*41 ).  As  a  result,  the  normalized  exchange  parameter  d 
is  a  constant. 

Returning  to  the  special  case  R^  (t,x)  by  inserting  parameters 
given  in  (3.2 6)  and  (3.28)  into  (3*37)  and  finally  into  (3*42),  the 
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threshold  of  instability  is  found  to  be 


(A  ) 
pyx 


(lwMs)2  “Wx'  +  sl“2  *o> 

"  /  cos  • 

“  S  *5  “  % 


(3-43) 


The  coordinate  subscripts  of  here,  as  well  as  hereafter,  refer  to 
the  direction  of  polarization  and.  propagation,  respectively,  of  the  un¬ 
stable  elastic  waves.  It  is  noteworthy  that  the  second  term  in  the  de¬ 
nominator  of  (3.43)  is  the  most  significant  in  lowering  the  threshold, 
because  x'  is  much  less  than  unity  and  the  angle  is  of  the  order 

of  several  degrees  (to  be  shown  in  the  next  section) — also  very  small. 
This  term  comes  from  the  modulation  of  the  spin  wave  mode  frequency 
o>m  of  (3.26),  as  can  be  shown  from  (3*37).  The  first  term  in  the  de¬ 


nominator  of  (3.43 )  is  not  completely  ignorable,  however,  and  comes 
from  modulation  of  the  coupling  between  the  spin  wave  and  acoustic 
wave  modes. 

The  procedure  for  determining  the  pumping  threshold  for  the 
remaining  types  of  elastic  waves  in  (3*19)  is  basically  the  same  as 
described  above  for  the  R  (t,x)  case.  The  frequency  of  the  spin  wave 

y 

mode  0^  is  the  same  for  all  cases  and  is  given  by  (3.26b).  The 

velocity  of  propagation  (o^/k)  of  the  elastic  wave,  however,  depends 

upon  the  mode  and  direction  of  propagation.  Hie  magnetoelastic  coupling 

parameters  K  ,  K  also  depend  upon  the  mode  and  direction  of  propa- 
mp  pm 

gation  of  the  elastic  wave,  and  hence  are  different  for  each  case 
considered. 


The  pump  threshold  angle  for  longitudinal  waves  R^(t,x)  propa¬ 
gating  in  the  x-y  plane  If 


(““V2  <2Ctn  sn  +  W(x'  +  **  io> 


(A  ) 

V  p  XX 


2x  ^  B”  sin  4  f  0  (  4  • 


(x'  +  sin  co*  3*, 


(3.44) 
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In  this  case,  modulation  of  the  spin  wave  mode  frequency  has  less  effect 
than  does  the  coupling  parameter  modulation.  The  reason  is  that  the 
magnitude  of  the  magnetoelastic  coupling  is  weaker  (for  small  angles) 
for  longitudinal  elastic  waves  than  it  is  for  transverse  elastic  waves. 

It  will  be  found  that  for  the  first  two  cases  the  lowest  threshold  occurs 
for  elastic  waves  propagating  nearly  along  the  applied  dc  magnetic  field 
(i.e.,  will  be  small). 

The  remaining  case  for  waves  propagating  in  the  easy  plane  is 
Rz(t,x)  ,  which  represents  a  transverse  wave  polarized  (i.e.,  with 
particle  displacement)  in  the  z-direction.  The  result  for  the  pump 
threshold  angle  is 


zx 


4Ct  ZT1Z  <•*»,?  (*'  *  P)  _ 

*  <ip  sin  a»c 


(J.^5) 


Here,  the  weak  magnetoelastic  coupling  as  well  as  the  influence  of  the 
anisotropy  "magnetic "  field  greatly  increases  the  threshold.  There  is 
only  one  term  in  the  denominator  of  (3»^3 )>  because  the  modulation  ef¬ 


fects  of  the  coupling  parameters  K 


and  the  spin  wave  frequency 


to  upon 
m 


mp  pm 

oo  ~  of  (3.32a)  combine  algebraically  into  a  single  term. 


PO 


Now  consider  the  two  cases  of  elastic  waves  propagating  along  the 
c-  or  z-axis.  First  take  Rx(t,z)  ,  which  is  a  transverse  wave  polarized 
in  the  easy  plane.  The  threshold  pumping  angle  is  given  by 


,  %  N,)  d+x'  +  p) 

(A  )  =  iZtZ - 2 - 


p  xz 


*  B4  sin  2 if 


0 


(3-U6) 


This  result  is  slightly  larger  than  (3*^3)  due  to  the  dependence  of  the 
spin  wave  frequency  upon  direction  of  propagation.  For  propagation  of 
spin  waves  normal  to  the  applied  dc  magnetic  field,  the  induced  dipolar 
demagnetizing  field  is  maxi mum --and  hence  the  spin  wave  frequency  is 
maximum  and  constant.  Therefore,  there  is  only  one  term  in  the  denominator 
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of  (5*^0  because  only  the  coupling  parameters  contribute  to  the  modu¬ 
lation  of  0)  -  .  For  the  thresholds  (A  )  and  (A  )  ,  however, 

pu  P  yx  p  ** 

the  modulation  of  w>  „  was  due  to  variations  in  oj  as  well  as  in  the 

pO  m 

coupling  parameters. 

The  last  case  is  for  R  (t,z)  ,  a  longitudinal  elastic  wave  propa- 

z 

gating  along  the  c-axis.  The  threshold  angle  for  this  case  is 


(3. *7) 


Actually  this  result  obtains  only  Insofar  as  is  neglected  in  com¬ 

parison  with  unity — a  condition  well  satisfied,  since  Ot ^  «  Q^y  due 
to  the  inherent  ellipticity  of  the  uniform  precession  mode.  Elastic  waves 
of  this  type  cannot  be  pumped  unless  there  is  some  deviation  of  the  uni¬ 
form  motion  of  the  magnetization  permitted  in  the  direction  of  the  par- 
ticle  displacement.  And  of  course  if  such  deviation  is  neglected,  then 
the  threshold  becomes  infinite. 


E. 


OPTIMIZATION  AND  CALCULATION  OF  PHONON  INSTABILITY"  1 


IQLDS 


1.  Minimization  of  Thresholds 


Each  pump  threshold  angle  given  by  relations  (3*^1)  through  (3*^) 

for  the  various  cases  is  a  function  of  ,  the  angle  between  the 

x-axis  and  the  applied  dc  magnetic  field.  To  detexmine  the  ml  nimiM 

threshold  angle  it  is  necessary  to  optimize  the  pumping  threshold  A^ 

with  respect  to  ,  subject  to  the  condition  (3*^1)  which  detexmlnes 

the  value  of  k  (and  hence  d).  The  optimum  propagation  angle  (t0^ppt 

is  found  for  each  kind  of  instability,  and  the  corresponding  optimum 

threshold  angles  (A  )  are  evaluated.  By  direct  co^>aiison,  the  mode 

p  cr 

having  the  lowest  threshold  is  then  found. 

Starting  with  the  simplest  case,  R  (t,x),  it  is  apparent  from  in- 

spection  that  (f  )  ,  =  n/U  ,  and  that  the  minimum  pumping  threshold 

u  op  v 

angle  is  therefore 


<Vz* 


cr 


V* 
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The  minimum  threshold  for  the  case  H  (t,z)  has  a  similar  form  and  is 
given  by 


(A  ) 
p  xz 


S«.«  (I“V 


cr 


SBi 


(3-^9) 


Optimization  for  the  R  (t,x)  and  R  (t,x)  cases  is  not  quite  so  obvious, 

X  y 

but  the  optimum  angle  nay  be  obtained  by  following  the  procedure  given 
in  the  first  paragraph  of  this  section.  A  rough  plot  of  (3.4l)  and  (3.^2) 
as  a  function  of  reve&ls  that  this  optimum  occurs  for  small  angles. 

Using  the  1 1  angle  approximations  for  the  trigonometric  functions 

in  (3*^3)  *nd  (3.44),  differentiation  shows  the  optimum  angle  to  be  given 
to  first  order  by 


)  as 

^O'opt 


(3.50) 


where  it  has  been  assumed  that  x'  «  1  .  Inserting  (3*50)  into  the 

small  angle  «  l)  approximations  to  (3*^3)  and  (3*44)  yields  the 
following  minimum  pimp  thresholds: 
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v  p  xx 
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(5-52) 


Note  that  these  two  minimum  thresholds  are  the  only  ones  which  depend 
significantly  upon  the  parameter  %'  . 
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2.  Comparison  of  Different  Fhonon  Threshold* 


A  comparison  of  the  above  phonon  thresholds  will  now  be  marie. 
Knowing  p  =  3*24  for  ZngY  ,  and  x'  «  1  ,  it  is  apparent  front  (3*1*8) 
and  (3*49)  that  (Ap)EX  <  (Ap)Xg  •  Assuming  the  factor  C/Q^B2  is 
approximately  the  same  for  each  threshold,  and  that  x'  «  1  ,  the 

following  approximate  ratios  between  the  other  thresholds  are  obtained: 


(A  ) 
v  p  zx 


(A  ) 

x  p  XX 


(AJvx  =  :  l  :  4x'. 
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pyx 
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(5-53) 


From  (5*53)  It  is  seen  that  the  highest  threshold  occurs  for  R  (t.jx) 
elastic  waves.  The  lowest  threshold  occurs  for  B  (t,x),  which  represents 

y 

transverse  waves  propagating  and  polarized  in  the  easy  plane.  As  was 
suspected  initially,  thresholds  are  substantially  higher  for  elastic 
waves  propagating  or  polarized  outside  of  the  easy  plane.  Bence  the 
thresholds  for  the  more  complex  quasi -transverse  and  quasi -longitudinal 
elastic  waves  will  not  be  considered. 


3.  Calculations  for  the  Lowest  Pumping  Thresholds 

The  magnetoelastic  constant,  the  elastic  constant,  and  the  phonon 
Q  of  Zn^Y  are  not  known.  However,  reasonable  estimates  can  be  made. 

Perhaps  the  best  estimate  can  be  made  for  the  elastic  constants, 
since  they  do  not  change  so  much  from  one  material  to  another.  The 
elastic  energy  expressed  in  terms  of  engineering  strains,  rather  than 
in  terms  of  the  pseudo-physical  strains  employed  above,  is  given  by 


6 

V?  I  'y'i'j  ■  '  <5-54.) 

M=i 


where  the 


are  the  elastic  constants,  and  the  strains 
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ere  related  to  the  tensor  strains  by 


22 


■  C11 

for 
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Expanding  (3.54a),  substituting  the  tensor  strains  given  by  (3»54b), 
and  caapar lng  the  resulting  elastic  energy  expression  in  terms  of  the 
c^s  with  (3»9)  gives  the  following  relations: 


C  =  cl..  C  *  id.. 

zzzz  53  |  tjzz  2  13 

ff  (S11+  ^  W  -  Hi*  (3.55) 

c«nn  =  £  ^cu  ' 

The  elastic  constant  C  will  be  estimated  by  taking  the  value  of 

1 p**  nn  /  5 

C. .  -  C-3  *  1.37  X  10  erg/ cnr  for  Zn  ,  which  hats  a  slightly  higher 
-U-  22 

density  than  Zn^Y  and  is  in  the  same  crystal  class.  From  (3.  35) 

the  elastic  constant  CL.  may  then  be  found  to  be 


C,{  -  3*3  x  1011  erg/cm3  .  (3.56) 


The  velocity  of  transverse  phonons  propagating  in  the  easy  plane  may 
now  be  found  from  (3*56)  and  (3*28b).  The  result  is  given  by 

1/2 

,  (3-57) 
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where  the  x-ray  density  given  in  Appendix  A  has  been  used  for  p 


-  70  - 


The  magneto elastic  constant  of  the  diagonal  strain  components  for 

*  2* 

cubic  crystal  YIG  is  ^  *  5  x  10  erg/car  .  The  msgnetoelaatic 
constants  for  hexagonal  Zn^Y  will  be  assumed  to  be  the  same  order  of 
magnitude  as  this,  giving  the  result 

Bj  **  5  x  10^  erg/ a?  .  (3. 58) 

6  2k 

The  acoustic  Q  of  Zn^Y  is  of  the  order  of  10  at  1  He. 

Assuming  the  acoustic  Q  varies  inversely  as  the  frequency,  an  acoustic 

Q  at  0)  /2x  *  0.66  kttc  is  estimated  to  be 
P 

*  1500  .  (3-59) 


Inserting  (3*59 )i  (3.58)  and  (3*56)  into  (3.52 ),  the  following 
numerical  expression  for  the  lowest  phonon  threshold  is  obtained: 


(apV* 


cr 
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(3.60) 


This  result  is  plotted  in  Fig.  3*3*  Figure  3*3  also  shorn  the  relationship 

of  h  to  (A  )  ,  which  is  the  same  here  as  in  Chapter  II. 

p  cr 

Using  the  estimated  constants  above,  the  assumption  made  that 
a^0  =  u>p  in  Equation  (3*4l)  is  found  to  be  valid  for  x'  >  0.001. 

even  before  this 


However,  the  theory  breaks  down  for  another 

is  reached.  In  deriving  the  above  phonon  thresholds, 

a)  approaches 


value  of  x 
2  2 

a>  «  o>  was  assumed.  At  the  lover  values  of  x 
p  m 

the  phonon  frequency  ,  and  this  assumption  no  longer  holds.  In 
particular,  when  cd^  =  3  the  actual  phonon  threshold  is  about 

2 Ofo  lower  than  the  value  given  by  (3»6o).  At  »  2 the  actual 
threshold  is  about  4c$  lower.  Belov  the  dc  magnetic  field  corresponding 
to  this  point  it  is  suggested  (belov)  that  spin  vave  thresholds  dominate 
the  phonon  threshold. 
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F.  COMPARISON  OF  THEORETICAL  THRESHOLDS  FOR  PHONON  AND  SPIN  WAVE 
INSTABILITIES 

The  minimum  phonon  threshold  shown  In  Fig.  3»3  is  a  function  of 
x'  [*(H  +  Dl?y4*M  ]  .  The  spin  wave  thresholds,  however,  shown  in 

Fig.  2.3  are  functions  of  x(»  H/WM  ),  rather  than  x'.  Consequently, 
before  a  comparison  of  the  phonon  and  spin  wave  thresholds  can  be  side 
it  is  necessary  to  know  the  value  of  the  normalized  exchange  parameter 
d  =  Dl^/WM  relating  x  and  x'(by  x'  *  x  +  d).  The  value  of  d  will 
be  found  experimentally  as  described  in  Chapter  VI.  Assuming  for  the 
useful  range  of  x  that  d  is  negligible  (i.e. ,  x  x'),  the  phonon 
threshold  of  Fig.  3*3  may  then  be  compared  directly  with  the  spin  wave 
thresholds  of  Fig.  2.3.  It  is  seen  from  this  comparison  that  the  mag¬ 
nitude  of  the  second  order  spin  wave  threshold  is  comparable  with  the 
magnitude  of  the  phonon  threshold  at  values  of  x  below  about  Xg  , 
where  o>m  =  2o^  *  o^,  .  On  the  other  hand,  if  d  is  not  negligible 

it  is  possible  for  a  field  region  to  exist  where  the  second  order  spin 
wave  threshold  will  dominate  the  phonon  threshold.  The  reason  for  this 
is  that  the  effect  of  a  finite  (positive)  value  of  d  is  to  increase 
the  phonon  threshold  for  a  given  field  value  x  .  This  effect  appears 
to  have  been  observed,  and  it  will  be  discussed  in  Chapter  VI. 

G.  SUMMARY 

Various  mechanisms  which  can  couple  the  uniform  precession  mode 
of  the  magnetization  to  elastic  waves  are  discussed;  namely,  the  magneto 
elastic  effect,  the  intrinsic  effect,  and  the  morphic  effect.  In  the 
calculations  of  this  chapter  it  is  ass  vised  that  the  magnetoelastic 
effect  is  the  dominant  mechanism  in  the  pumping  of  elastic  waves  by 
the  uniform  precession.  The  spin  wave  equations  of  motion  of  Chapter  II 
were  modified  to  include  contributions  from  the  magnetoelastic  energy 
which  couples  energy  from  the  uniform  precession  to  elastic  waves. 
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The  elastic  vave  equations  follow  from  the  elastic  energy,  together  with 
the  magnetoelastic  energy  expression.  Hie  simple  solutions  to  the  elastic 
wave  equations  are  the  ones  considered  in  the  elastic  wave  instability 
calculations. 

In  order  to  take  into  account  propagation  of  elastic  waves  in  the 
easy  plane  at  an  arbitrary  angle  with  respect  to  the  direction  of  the 
dc  magnetic  field,  the  equations  of  motion  are  expressed  in  terms  of 
a  rotated  coordinate  system.  The  elastic  displacement  components  were 
left  in  the  fixed  coordinate  system  in  order  to  maintain  the  identity 
of  the  simple  elastic  modes.  For  each  simple  mode  case  the  elastic 
and  spin  wave  equations  were  combined,  and  a  Matthieu  type  of  equation 
follows.  For  each  case  the  first  order  elastic  wave  (or  phonon) 
threshold  was  determined. 

The  lowest  threshold  was  found  to  exist  for  transverse  waves  pro¬ 
pagating  and  polarized  in  the  easy  plane  of  the  ferrite.  After  making 
estimations  of  the  elastic  constant,  the  magnetoelastic  constant,  and 
the  phonon  Q,  ,  and  knowing  values  for  the  ferrite  parameters,  threshold 
values  were  calculated  as  a  function  of  the  effective  internal  dc 
magnetic  field  (i.e.,  internal  plus  exchange  fields).  The  threshold 
was  found  to  increase  more  rapidly  than  linearly  with  increasing  internal 
field.  Assuming  the  exchange  contribution  to  the  effective  internal 
field  is  small,  a  comparison  of  the  lowest  phonon  threshold  was  made 
with  the  spin  wave  thresholds  of  Chapter  II.  From  this  comparison  it 
was  expected  that  the  phonon  threshold  would  dominate  the  spin  wave 
thresholds  at  internal  fields  greater  than  that  corresponding  to  the  point 
where  the  second  order  off -resonance  spin  wave  threshold  begins.  Below 
this  point  the  phonon  and  spin  wave  thresholds  were  expected  to  be 
comparable. 
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CHAPTER  IV 


THEORETICAL  POKER  ABSORPTION  BY  PHONON  INSTABILITIES 

A.  DIVISION  OF  POKER 

Suppose  microwave  (rf)  power  P^  is  Incident  upon  an  electromagnetic 

resonant  system  (dielectric  resonator)  which  contains  a  ferrite.  Of  this 

incident  pump  power  P.  a  certain  portion  P  is  absorbed  by  the  reso- 

nant  system.  This  absorbed  power  sustains  the  rf  magnetic  field  in  the 

dielectric  resonator  by  overcoming  losses  in  the  dielectric,  the  ferrite, 

and  in  the  waveguide  walls.  For  a  transverse  pumping  experiment,  the 

ferrite  loss  comprises  mainly  the  loss  associated  with  damping  of  the 

magnetic  precession.  This  loss  may  be  due  to  tvo-mognon  pit  scattering, 

25 

rapidly-relaxing  inpurity  or  Kasuya-LeCrav  confluence  processes. 

As  the  rf  power  level  is  increased,  additional  ferrite  loss  will 
occur  due  to  spin  wave  or  elastic  wave  instabilities.  In  such  a  case 
a  transient  or  irregularity  can  be  observed  on  the  reflected  pump  pulse 
when  the  absorbed  power  (or  rf  magnetic  field  amplitude)  has  reached  a 
critical  level.  In  Chapter  II  it  was  shown  that  spin  wave  instabilities 
might  occur  when  the  pump  frequency  is  within  the  spin  wave  manifold, 
and  phonon  instabilities  should  occur  when  the  puap  frequency  is  below 
the  spin  wave  manifold.  This  chapter  is  concerned  with  the  prediction 
of  the  transient  loss  due  to  the  growth  of  phonon  instabilities. 

B.  GROWTH  OF  PHONON  INSTABILITIES  FROM  THERMAL  LEVEL 

In  Chapter  III  it  was  shown  that  elastic  waves  are  parametrically 
pumped  by  the  uniform  precession  mode.  Nhen  the  portion  of  the  pump 
power  going  into  the  growing  elastic  waves  becomes  equal  to  the  power 
being  lost  through  the  elastic  wave  damping,  then  the  elastic  wave  or 
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phtonon  threshold  is  reached.  This  threshold  was  calculated  in  Chapter  III. 
It  is  clear  that  below  the  phonon  threshold  elastic  wave  growth  will  be 
damped  out.  Above  the  phonon  threshold,  however,  there  is  a  net  power 
going  into  the  elastic  waves  and  the  latter  become  unstable.  Finally, 
at  a  certain  level  above  threshold,  the  power  going  into  the  growing 
elastic  waves  will  reach  a  sufficiently  high  level  that  it  can  be  observed 
as  a  transient  or  irregularity  in  the  reflected  pump  pulse  waveform. 

.  This  occurs  at  the  critical  power  level  referred  to  in  Section  A  above. 

1.  Solution  to  the  Elastic  Wave  Equation 

It  was  found  in  Chapter  III  that  elastic  waves  having  the  lowest 
threshold  are  transverse  waves  propagating  and  polarized  in  the  easy 
plane  and  at  a  frequency  of  one  half  the  pump  frequency.  The  basic 
differential  equation  governing  the  motion  of  these  elastic  waves  was 
given  by  (3.35)*  For  pumping  angles  at  or  below  the  phonon  threshold, 
the  solution  had  the  form  given  by  (3*16)  and  (3.1Tb).  In  order  to  take 
into  account  the  growth  of  unstable  elastic  waves  above  threshold,  (3.16) 
is  modified  to  include  the  exponential  growth  parameters  as  follows: 


(iaj+s)t 

ky  S 


+ 


*  (-iu*s)t 
-ky  e 


(4.1) 


This  expression  for  the  elastic  displacement  will  be  inserted  into  the 
elastic  wave  equation  (3*35)  to  determine  the  growth  parameters  for 
operation  above  threshold;  but  first  (3*35)  will  be  modified  to  include 
the  time  dependence  of  the  rf  magnetic  field  amplitude. 

In  deriving  the  phonon  thresholds  of  Chapter  III  it  was  assumed  that 
the  pump  angle  A^  was  constant.  This  implies  the  power  pulse  incident 
upon  the  ferrite  system  had  the  ideal  pulse  characteristics  of  perfect 
flatness,  zero  rise  apd  fall  times,  and  no  frequency  modulation.  Such 
an  ideal  pulse,  sketched  in  Fig.  4.1a,  could  not  be  obtained  in  practice 
chiefly  because  the  required  pulse  length  was  of  the  order  of  100  psec. 
In  general,  it  is  difficult  to  obtain  a  flat  pulse  of  this  length  with 
existing  techniques,  particularly  at  the  frequepcy  (1.32  kMc)  and  power 
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a  b 

FIG.  4.1 — Sketch  of  tine  dependence  of  the  incident  puqp  power  for  the 
ideal  (a)  and  actual  (b)  pulses . 
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FIG.  4. 2- -Sketch  of  tine  dependence  of  the  rf  negnetic  field  anplltude 
for  the  ideal  (a)  and  actual  (b)  pulse  cases. 


FIG.  4.3 — Sketch  of  tine  dependence  of  the  rf  nagnetic  field  amplitude 
when  the  apparent  threshold  is  reached  for  the  Ideal  (a)  and 
actual  (b)  pulse  cases.  The  apparent  critical  field  h  end 
build-up  tine  t  indicate  the  initial  appearance  of  ffie  power 
absorption  due  to^grovth  of  phonon  instabilities. 
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level  (~  lOW  peak)  desired.  (Refer  to  Chapter  Y  for  details  on  the 
technique  employed  to  obtain  the  100  jisec  pulse;  refer  to  Chapter  VI 
for  photographs  of  the  actual  power  pulse  waveforms  used  in  the  threshold 
experiments.)  13ie  incident  power  for  the  actual  pulse  is  sketched  in 
Fig.  4.1b  in  terms  of  the  approximate  representation 


(4.2) 


where  t  *  0  specifies  the  time  when  the  pulse  has  "settled  down"  to  a 
frequency  within  the  bandwidth  of  the  resonant  system;  Piffl  is  the 
incident  power  at  t  *  0  ;  and  t  is  the  time  constant  during  the 

useful  length  of  the  pulse,  which  exbludes  the  frequency  modulation 
at  the  end  of  the  pulse.  Since  the  reflected  power  level  is 

related  to  the  incident  power  P.^  by 


where  T  is  the  reflection  coefficient,  then  the  absorbed  power  Pa  is 
given  by 


Pa  -  (i  -  |r|2)  ^  . 

Chapter  V  shows  that  the  resonator  magnetic  field  amplitude  h  is  related 

to  the  absorbed  power  P  as  follows: 

gl 


h 
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Using  the  above  relations  between  h  ,  P  ,  r  and  P.  ,  then,  the  tiae 

ft  1 

dependence  of  the  field  amplitude  is  found  with  the  aid  of  (4.2)  to  be 

h  «  h^  e-t/2T  ,  (4.3) 

where  h^  is  the  initial  Mariana  Magnetic  field  amplitude  corrt  aponding 
to  Pia  ,  and  r  has  been  assuaed  to  be  constant  during  the  pulse. 

The  ideal  and  actual  field  amplitude  pulse  shapes  reaching  the  ferrite 
sample  are  sketched  in  Pigs.  4.2a,b.  Bie  actual  pulse  approaches  the 
ideal  pulse  in  flatness  as  the  tiae  constant  2t  becoaes  infinite.  The 
actual  pulse  shape  is  considered  here  rather  than  an  agproxiaation  _ 
thereto  because  it  will  be  found  that  aore  inf oraation  regarding  the 
nature  of  the  phonon  absorption  can  be  obtained  theoretically  fron  the 
actual  pulse  shape. 

Rie  puaplng  angle  was  defined  by  (2.28b)  as 

Ajp  -  h/H0  -  (*.<0 


In  Chapters  II  and  III,  where  h  was  assuaed  constant,  the  threshold 
pumping  angle  or  critical  angle  was  given  by 


<Vcr  *  hcA>  • 

where  hcj>  is  the  critical  or  threshold  rf  field  aaplitude  and  HQ  is 
the  dc  field  at  the  ferrite  saaple.  In  this,  chapter  where  rf  field 

V-  •-  • 

amplitudes  greater  than  h^  are  considered,  the  threshold  angle  is 
given  by 
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9 


(*.6) 


'*««  <VaP  “d  h  refer  to  the  apparent  threshold  angle  and  ap¬ 
parent  critical  field,  respectively,  at  which  the  power  absorption  due 
to  the  growth  of  instabilities  is  sufficiently  high  to  be  seen  as  an 
irregularity  in  the  reflected  pimp  pulse  (Fig.  4.3). 

For  the  actual  pulse  the  apparent  critical  field  was  determined  by 
Measuring  the  power  level  of  the  incident  pulse  at  the  time  t  where 
the  instability  first  appeared  on  the  reflected  pulse.  'Rius  the  apparent 
critical  field  for  the  actual  pulse  is  given  by 


-t  _/2t 

(actual  pulse)  h  *heF  ,  (4.7) 

ap  op 


where  h^  is  the  corresponding  value  of  hffl  .  For  the  ideal  pulse, 
where  t  ,  the  apparent  critical  field  becomes 


(ideal  pulse)  h  *  h  ,  (4.8) 

ap  mp 

and  the  instability  would  first  be  observed  at  the  end  of  the  pulse, 
t  *  T  .  The  relationships  between  these  quantities  expressed  in 
(4.8)  and  (4.7)  for  the  ideal  and  actual  pulse  shapes  are  shown  in 
Figs.  4.3a  and  b.  The  dashed  lines  indicate  the  irregularities  due  to 
the  instability. 

Inserting  the  time  decaying  pumping  angle  (4.4),  with  h  given 
by  (4.3),  into  the  elastic  wave  equation  of  motion  ( 3 -  35 )  gives 


PO 


COS 


vK"° 


(4.9) 
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A  trial  solution  (4.1)  with  s  «  s(t)  Is  now  inserted  into  (4.9)*  giving 


2s  ♦  ts  +  (s  +  ts)2  +  2tj(b  +  ts)  +  od^q  -  o>2 

(4.10) 

+  21»(s  ♦  ts  +  n)J  +  r*^  -  0  , 

and  its  adjoint  (complex  conjugate,  with  k  -*■  -k)  ,  where  the  frequency 

condition 

oi  *  oip/2  (4.11) 

has  been  Invoked.  The  secular  equation  for  this  will  require  that 
|r^(  *  jr  ^1  ,  which  is  the  requirement  for  a  standing  wave  solution 

as  mentioned  in  Chapter  III  /especially.  (J.lB). 

For  reasons  to  be  discussed  later,  only  resonant  instabilities  will 
be  considered  in  calculating  the  power  absorbed  by  the  pimp.  The  real 
part  of  the  coefficient  for  r^  in  (4.10)  is  therefore  assumed  to  be 
very  nearly  zero  at  all  times,  giving 


2hu(s  +  ts  +  ♦  <*3^0  r!*y  *  0  *  (4-^) 


The  secular  equation  of  (4.12)  and  its  adjoint  in  this  case  yields 


ts  +  8  +  ij 


c.(+  \ 


(4.13) 
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when  written  in  terms  of  the  critical  field  h  ,  specified  in  general 
by  (4.5)  and  (3.42)  and  written  as 


(k.lk) 


where  [(f  is  the  optimum  value  of  the  modulation  parameter 

G(t0)  ,  and  the  phonon  is  related  to  the  relaxation  fi*e(iuency 

n  by  tj  *  ^po/2^  *  If  is  a  slowly  varying  function  of  time, 

the  MS  of  (4.13)  W  be  expanded  in  a  Taylor  series.  A  series 
solution  for  the  exponential  growth  parameter  s  is  thereby  found  to  be 


a(*0  .  t)  -  -  n  +  n 


G(*J  h'(t) 


fit) 


cr 


(4.15a) 


where 


h'<t)  .  h  £  — - —  (g\* 

“„to  <d+1>d!  V  / 


m 


‘•«T* 


h(i)' 


(4.15b) 


is  an  effective  rf  field  amplitude.  The  series  converges  rapidly  since 
t  «  2t  ,  and  only  two  or  three  terms  are  needed.  For  small  t/T  , 
h'(t)  decreases  with  a  time  constant  approximately  one  half  that  of 
h(t).  . 

2.  Power  Absorbed  by  the  Unstable  Modes  (Rionon  Absorption) 

Assuming  that  all  of  the  pump  energy  which  reaches  the  elastic  wave 
system  goes  into  pumping  the  unstable  phonon  modes,  the  total  elastic 
energy  may  then  be  represented  (within  a  constant)  by  the  energy  contained 
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only  In  the  unstable  modes .  Thus  the  elastic  energy  density  for  trans¬ 
verse  acoustic  vares  propagating  and  polarized  in  the  easy  plane  is  found 
from  (3-9)  to  be 


X  -  *Ctv  t  , 

p  A  inn  *y  ' 


.  .  6) 


vhere 


1  dR 

V  ;  if 


Utilizing  the  expansion 


■v-X 


(4.17) 


vhere 


R-ky  “  \f 


the  time  average  stored  energy  in  the  sample  is  found  for  standing  wave 


modes  to  be 


C  -  /  S  dV  - 
P  J  P 


cunn  TX  k2||r*/  +  lr-^.|2j  -  <4'l8> 


vhere  a  different  grovth  factor  has  been  allowed  for  each  mode,  and  V 
is  the  volume  of  the  sample. 
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Since  the  unstable  nodes  will  begin  growing  exponentially  from  thermal 
lerel  it  t>0  ,  the  are  rage  thermal  energy  in  the  modes  must  be  found  in 

order  to  evaluate  the  energy  at  later  times.  For  microwave  frequencies 
to  Is  much  less  than  k^T  (where 

classical  statistics  may  be  used.  Thus  the  average  thermal  energy  con 
taineA  in  each  mode  is 


kg  is  Boltzmann's  constant ), 


The  total  average  energy  (4.l8)  may  be  written 


(4.19) 


(4.20) 


since  there  are  two  standing  wave  modes  (sin  kx  apd  cos  kx)  for  each 

value  of  k  .  The  energy  will  be  used  to  evaluate  the  power 

absorbed  from  the  pmg>.  Equation  (4.20)  should  therefore  be  summed  only 

over  the  mode#  which  are  growing  in  time. 

In  determining  the  values  of  k  which  will  contribute  significantly 

to  the  suMsatioo  in  (4.20)  it  is  helpful  to  utilize  the  k-space  concept. 

Figure  4.4  shows  a  representation  of  a  portion  of  k- space  for  periodic 

28 

boundary  conditions.  Eae  distance  separating  modes  along  a  given  co¬ 
ordinate  direction  1  in  k-space  is 


dL  -  2*/*!  ,  (4.21) 


where  a^  is  the  dimension  of  the  sample  in  the  direction  1  For  a 
ferrite  disk,  the  mode  spacing  may  be  approximated  by 


(4.22) 
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in  the  plane  of  the  disk,  vhere  d  is  the  disk  diameter s  and  by 


2*/\ 


(4.23) 


normal  to  the  disk,  vhere  v  is  the  disk  thickness. 

In  Chapter  III  it  was  assumed  that  there  vas  only  one  unstable  mode, 

resonant  at  oip/2  .  This  mode  is  represented  by  a  single  point  in 

k-space,  the  magnitude  of  the  wave  vector  k^  being  determined  by  the 

frequency  condition  (3*41)  and  the  orientation  of  k^  by  the  direction 

of  propagation  relative  to  the  dc  magnetic  field.  In  Chapter  III  an} 

optimum  value  vas  found  for  the  angle  betveen  k  and  the  dc  field. 

For  convenience,  the  axes  in  k-space  vill  be  chosen  so  that  k  cor- 

responds  to  this  optimum  direction,  and  the  k  ,k  plane  corresponds 

2c  y 

to  the  plane  of  the  disk  sample  (see  Fig.  4.4).  Although  the  assump¬ 
tion  of  a  single  mode  is  a  valid  one  in  the  determination  of  the  phonon 
threshold  angle  or  critical  rf  magnetic  field,  it  is  not  valid  for 
operation  above  threshold.  The  reason  is  that  other  modes  vill  be  pumped. 
This  case  vill  nov  be  considered. 

The  first  restriction  on  the  modes  vhich  vill  contribute  to  (4.20) 
is  that  <bp/2  should  be  vithin  the  mode  bandvidth.  Otherwise  they  vill 
grow  only  very  slowly,  if  at  all.  This  bandvidth  is 


Vs 


(4.24) 


and  from  (3-28b)  the  corresponding  “wavenumber  bandwidth"  is 


*4  '  k<A>«k0 


(4.25) 
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This  shows  that  the  relevant  modes  must  lie  within  a  thin  spherical  shell 
of  radius  kQ  and  thickness  Ak  in  k-space  (see  Fig.  4.4).  In  order  to 
compare  the  “wavenumber  bandwidth”  with  the  mode  spacing,  typical  para¬ 
meters 


cttj/4x  »  0.66  kNc 

d  -  0.054  in. 

3  4 

Qp  *  10-5  -  10 

are  chosen.  It  can  then  be  verified  that 


Ak  « 


> 


as  shown  in  Fig.  4.4.  The  radius  of  the  shell,  k^  ,  is  not  shown  to 
scale.  Using  the  parameters  assumed  above  results  in 


1^  -  150  |i  -  (103  -  104)  dk 


(*•26) 


and  the  thickness  of  the  shell  is  very  small  compared  with  its  radius. 

A  more  significant  restriction  on  the  modes  relates  to  the  dependence 
of  the  modulation;  parameter  G  on  the  angle  .  This  limits  the 

contributing  modes  to  a  portion  of  the  above  mentioned  spherical  surface 
in  the  neighborhood  of  the  axes  k^  .  Figure  4.5  shews  a  plot  of  the 
threshold  angle  (A  )  vs  with  a  typical  experimental  value  of 

*  w  ^  v 

x  *  0.05  and  parameter  values  as  specified  by  Chapter  III.  The  angle 

(Vopt  indlcated  calculated  fro*  the  s»ll  angle  approximation 

(3*50)*  It  is  seen  to  agree  closely  with  the  actual  minimum  of  (Ap)^  • 

The  threshold  becomes  infinite  at  *  s/4  .  The  magnitude  of  the 

threshold  (a  i  as  a  function  of  is  shown  in  Fig.  4.6a.  Due  to 
p  u 
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FIG.  4.5 — Threshold  (A  )  vs.  rotation  angle  for  resonant  modes 

p  ^ 

(cn  *  CDpQ  -  0ip/2)  .  -  86  - 


FIG.  4.6a — Sketch  shovinr  one  cycle  variation  In  the  Magnitude  of  the 
resonant  threshold  angle  |  (A  )  |  vs  rotation  angle  ♦ 
The  function  Is  repetitive  vi$hr$eriod  ac  .  J 


FIG.  4. 6h- -Sketch  shoving  one  cycle  variation  of  the  aodulation 
parameter  Magnitude  |g(#q) |  corresponding  to  the 
threshold  variation  of  Fig.  4.6a. 


the  inverse  dependence  of  A  upon  sin  ,  A  has  singularities 

nx  *  ^  * 

at  ♦q  *  $  vhere  n  is  any  integer.  Figure  4.6b  shows  the  modula¬ 

tion  parameter  G(<Tq)  »  which  varies  inversely  with  A^  ,  in  accordance 
with  (3*42).  Note  the  narrow  range  of  over  which  G  is  large. 

This  means  that,  when  the  rf  pump  field  is  only  slightly  greater  than 

h  ,  the  growing  modes  in  the  k  ,k  plane  of  Fig.  4.4  lie  within  a 
cr  x  y 

range  of  angles  about  %>oPt  • 

The  threshold  angles  shown  in  Figs.  4.5  and  4.6  apply  only  to  trans¬ 
verse  inodes  with  the  wave  vector  and  polarization  in  the  plane  of  the 

disk  (k  ,k  in  Fig.  4.4).  Modes  lying  outside  the  k  ,k  plane  in 
x  y  x  y 

Fig.  4.4  should  also  be  considered  in  evaluating  (4.20).  Furthermore, 
there  will  be  two  transverse  and  one  longitudinal  mode  associated  with 
each  point  in  Fig.  4.4,  and  these  should  all  be  included  in  the  summa¬ 
tion.  It  was  shown  in  Chapter  III,  however,  that  the  threshold  angle 
is  very  much  larger  when  the  wave  vector  is  normal  rather  than  parallel 
to  the  plane  of  the  disk.  On  the  basis  of  this  result,  and  the  wider 

spacing  of  the  modes  along  the  k  axis  (Fig.  4-4),  Eq.  (4.20)  will 

z 

only  be  summed  over  the  k  ,k  plane.  Since  it  was  also  shown  in 

x  y 

Chapter  III  that  the  transverse  wave  polarized  in  the  plane  has  the  . 
lowest  threshold,  the  restriction  of  the  summation  in  (4.20)  to  modes  of 
this  kind  will  be  retained. 

From  the  above  considerations  it  is  seen  that  the  summation  in  (4.20,) 

will  be  taken  only  over  modes  lying  within  the  intercept  of  the  spherical 

shell  upon  the  k  ,k  plane  in  Fig.  4.4.  Within  this  region  the  modes 

x  y 

are.  distributed  with  an  approximately  uniform  angular  spacing, 


2„ 

—  <  — 4  «  l 


radian 


kod  k0 


(from  4.26);  and  there  will  be  only  one  mode  at  each  angular  position. 
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since  Ak  is  such  less  than  2x/d  .  The  summation  in  (4.20)  aay  therefore 
be  approximated  by  an  integral 


♦a^Vopt 


'*BTk0 


26lr(Hf0,k,t)t 

e 


■(+0>opt 


(lt.27) 


It  is  not  necessary  to  consider  negative  values  of  k^  ,  since  tvo 

standing  wave  modes  for  each  value  of  k  have  been  allowed  for  in  (4.20). 

In  obtaining  (4.27)  from  (4.20),  a  factor  of  2  has  been  taken  out  of  the 

summation  to  account  for  phonon  modes  which  occur  for  k  near  «  -(♦_)  , 

u  o  opt 

as  well  as  near  t  •  Figure  4.6  shows  that  these  values  of 

are  the  ones  near  which  G(tg)  #  and  hence  s(tg)  ((4.15a)],  has 
its  maximum  values. 

Since  all  modes  considered  in  (4.27)  are  close  to  resonance,  the 
growth  parameter  for  a  resonant  mode  s(t0,kQ)  will  be  used  in  evaluat¬ 
ing  the  integral.  Also,  the  first  term  of  a  Taylor  series  expansion 
about  (^o^opt  wi^-  ^  1186(1  for  the  modulation  parameter 
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and 


8k0(VLo't) 


n  -  i  + 


—  U-±< 

\r  L  \ 


♦o  '  <*o>oPt 


r 


•  (^-30) 


Using  (4.27)  and  (4.30)  gives 


♦2  =sm h'  ft 

%T  kQ  2tjt(-l+»7hcr)  ^  2hcr  x'  (f0 


^O^opt 


r 

) 


a*. 


:il 


(4.31) 


For  t  greater  than  l/tj  and  ♦0“^O^opt  of  the  order  fy^opt  1 

the  integrand  of  (4.31)  becomes  vanishingly  small,  and  the  integration 

limits  may  be  extended  to  infinity  without  significantly  affecting  the 

value  of  the  integral — in  accordance  with  the  usual  procedure  for  the 

27 

method  of  steeppst  descent.  The  integral  reduces  to 


2«x'  h 


M2 


cr 


9nt  fa' 


(^•32) 


From  (4.20),  (4.31),  and  (4.32)  the  total  stored  energy  in  the 
growing  modes  is 


BkgT  e 


2nt(-i  +  h7hcr) 


(**•33) 
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where 


,  (t  >  l/n)  ,  (it.3^) 


is  the  effective  number  of  growing  modes.  Equation  (4.34)  shows  that 
the  effective  number  of  growing  modes  decreases  with  time,  even  when  the 
pump  field  amplitude  h'  is  constant.  Prom  (4.33)  the  effective  nuriber 
of  modes  is  defined  relative  to  the  marl—  m  growth  rate.  As  time  increases, 
the  modes  with  small  growth  rates  become  less  and  less  significant  by 
comparison  with  the  rapidly  growing  modes,  and  the  number  of  effective 
modes  decreases.  When  the  condition  t  >  l/tj  is  not  satisfied,  the 
method  of  steepest  descent  is  not  valid.  Fortunately,  the  experimental 
effects  are  observed  at  times  satisfying  this  condition. 

Knowing  the  phonon  energy  as  a  function  of  time,  the  total  power 
Pp  going  into  the  phonon  system  may  be  determined  from 

ae 

p_  -  — £  ♦  2iie  .  (**-35) 

p  dt  p 


The  first  term  in  (4.35)  represents  the  increase  in  stored  energy.  If 
the  stored  energy  were  not  increasing  there  would  still  he  a  power  loss 
in  the  phonon  system.  This  additional  power  loss,  represented  by  the 
second  term  of  (4.35),  is  associated  with  the  phonon  ftp  .  Inserting 
(4.33)  into  (4.35)  yields 
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Since  h'  is  a  slowly  varying  function  of  time,  it  may  be  taken  as  a 
constant  in  (4. 36)  except  in  the  exponential  factor  of  6^  .  Assuming 
that  t  is  greater  than  a  few  jisec  [(l/4ijt)  «  1]  ,  the  phonon  power 
absorption  is  given  approximately  by 


where  N(t)  is  approximately 


H(t.) 


2x'  h 


cr 


JCTjt  h. 


m 


(4.37a) 


(4.37b) 


The  damping  parameter  n 


is  related  to  the  phonon 


%  by 


11  *  *Pc/2S> 


(4.37c) 


C.  PHONON  ABSORPTION  IN  TERMS  OF  MEASURABLE  PARAMETERS 


It  is  now  necessary  to  relate  the  phonon  absorption  P^  to  the 
transient  observed  in  the  reflected  pump  pulse.  From  simple  conservation 
of  power,  the  reflected  power  Pp  is  given  by 


P 

r 


P 

a 


(4.38) 
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where  P^  is  the  incident  power  and  P&  is  the  absorbed  power.  Given 
that  P.^  is  constant,  the  differential  relation  which  follows  from  (4.38), 
namely 


(*•39) 


shows  that  any  increase  which  occurs  in  P  (due  to  phonon  absorption) 
will  be  accompanied  by  a  corresponding  decrease  in  P^  . 

The  power  P  absorbed  by  the  dielectric  resonator  system  is  given 
in  terms  of  the  stored  energy  U  by  the  familiar  expression 


P 

a 


(4.4o) 


where  oip  is  the  pump  frequency  and  is  the  effective  Q  of  the 

resonator  system.  Since  the  absorbed  power  is  a  function  of  the  stared 
energy  as  well  as  of  the  system  Q  ,  the  following  differential  relation¬ 
ship  is  obtained  from  (4.40): 


AP 


fid 


+  AP. 


(4.41a) 


where  represents  the  additional  power  loss  required  to  sustain 

the  Increased  energy  of  the  electromagnetic  field,  and  is  given  by 


Oip 

APfl.  -  ~  au  .  (4.4lb) 

<*0 
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The  increase  in  ferrite  power  absorption  AP^,  is  given  by 


(4.4lc) 


This  power  absorption  takes  place  only  when  the  of  the  system  is 

changed  by  introducing  additional  losses  into  the  system — for  instance, 
through  growth  of  spin  wave  or  elastic  wave  instabilities  in  the  ferrite. 
For  an  overcoupled  resonator  system,  the  external  Q  is  given  by 

«ext  *  Vr0  ' 

where  rQ  is  the  voltage  standing  wave  ratio  (VSWR)  at  resonance. 

For  a  given  position  of  the  movable  short,  ^e3C^  is  constant  (see 
Chapter  V).  Therefore  (4.4lc)  becomes 

AP  .  -Pa—  .  (4.1*3) 

ro 

It  is  desired  to  obtain  the  ferrite  absorption  in  terms  of  the  reflected 

power  deviation  AP  .  The  reflected  power  P  is  related  to  the 

r  r 

incident  pcwer  by  the  familiar  expression 


P 

r 


(4.44) 


« 


4 


I 


Where  the  magnitude  of  the  reflection  coefficient 


|r|  is  defined  by 


ro  +  1 


(4. 1*5) 


f 
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By  taking  the  differential  of  (4.4k)  with  respect  to  ,  it  can  be 

shown  that 


(4.46) 


Combining  this  result  with  (4.43)  yields  the  following  for  the 
change  in  the  power  absorbed  by  the  ferrite: 


-AP 

r 

Irl 


Consider  now  that  the  change  in  ferrite  absorption  AP^  represents 
the  increasing  transient  power  loss  due  to  the  growth  of  instabilities 
within  the  ferrite.  Since  the  growth  of  phonon  instabilities  is  assumed 
to  be  the  dominant  transient  loss  mechanism,  it  may  be  concluded  that 


APF  -  Pp  .  (4.48) 

It  is  seen  from  (4.47)  that  AP„( *P  )  Is  proportional  to  ,  the 

r  p  r 

observable  change  in  reflected  power,  provided  the  reflection  coefficient 
T  does  not  change  appreciably.  Since  the  resonant  system  is  overcoupled, 
-APr  will  be  positive  and  increasing  as  the  ferrite  losses  increase. 

Hence  the  reflected  power  waveform  will  show  a  decrease  with  increasing 
time  due  to  the  phonon  absorption. 

The  ideal  situation  from  the  standpoint  of  observation  of  instability 
thresholds  is  to  have  a  marl  mum  relative  change  in  reflected  power.  Thus 
combining  (4.44)  with  (4.47)  and  (4.1j8)  yields 


<*.J»9) 
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The  increase  in  ferrite  power  absorption  AP^,  is  given  by 


?F  *  ^ 


■CO 


(4.4lc) 


This  power  absorption  takes  place  only  when  the  Qq  of  the  system  is 
changed  by  introducing  additional  losses  into  the  system — for  instance, 
through  growth  of  spin  wave  or  elastic  wave  instabilities  in  the  ferrite. 
For  an  overcoupled  resonator  system,  the  external  Q  is  given  by 


(4.42) 


where  rQ  is  the  voltage  standing  wave  ratio  (VSWR)  at  resonance. 
For  a  given  position  of  the  movable  short,  Qext  is  constant  (see 
Chapter  V).  Therefore  (4.4lc)  becomes 


(4.43) 


It  is  desired  to  obtain  the  ferrite  absorption  in  terms  of  the  reflected 


power  deviation  AP 


The  reflected  power  P  is  related  to  the 

r 


incident  power  by  the  familiar  expression 


r|2  p. 


(4.44) 


Where  the  magnitude  of  the  reflection  coefficient  |r|  is  defined  by 


r0  ~  1 
r0+  1 


(4.45) 
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it  can  be 


By  taking  the  differential  of  (4.44)  with  respect  to  r q  , 
shown  that 


(4.46) 


Combining  this  result  with  (4.43)  yields  the  following  for  the 
change  in  the  power  absorbed  by  the  ferrite: 


-AP 

ap_  -  — p  .  (u.vr) 

F  lr| 

Consider  now  that  the  change  in  ferrite  absorption  APp,  represents 
the  increasing  transient  power  loss  due  to  the  growth  of  instabilities 
within  the  ferrite.  Since  the  growth  of  phonon  instabilities  is  assumed 
to  be  the  dominant  transient  loss  mechanism,  it  may  be  concluded  that 


4Pf  -  pp  •  (1».1t8) 

It  is  seen  from  (4.1*7)  that  <2iP_(«P  )  is  proportional  to  AP  ,  the 

r  P  r 

observable  change  in  reflected  power,  provided  the  reflection  coefficient 
T  does  not  change  appreciably.  Since  the  resonant  system  is  overcoupled, 
-APr  will  be  positive  and  increasing  as  the  ferrite  losses  increase. 

Hence  the  reflected  power  waveform  will  show  a  decrease  with  Increasing 
time  due  to  the  phonon  absorption. 

The  ideal  situation  from  the  standpoint  of  observation  of  instability 
thresholds  is  to  have  a  maximum  relative  change  in  reflected  power.  Thus 
combining  (4.44)  with  (4.47)  and  (4.46)  yields 


P 

r 


|r| 


(4.49) 
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for  the  relative  change  in  reflected  power.  The  incident  power  is 

determined  by  the  threshold  magnetic  field  amplitude  at  which  the  phonon 
absorption  reaches  an  observable  level.  Thus  only  |r|  remains 

to  be  chosen.  It  is  clear,  therefore,  that  the  maximum  relative  change 
in  reflected  power  will  occur  for  the  smallest  reflection  coefficient 
compatible  with  the  minimum  detectable  reflected  power. 


D.  APPARENT  CRITICAL  FIELD 


The  relationship  between  the  phonon  power  absorption  P^  and  the 

apparent  critical  magnetic  field  amplitude  h  will  now  be  determined. 

ap 

1.  Idc^l  Pulse 


’  00 


since  the  ideal  pulse  is  assumed  to  be 


The  phonon  absorption  P^  for  the  ideal  pulse  may  be  found  from 
(4.37a)  by  lettinp  t 
perfectly  flat. 

tion  when  saturation  of  the  reflected  pump  pulse  first  appears,  and 
recalling  that  h  *  h  (and 


On  defining  P  *  (P  )  as  the  observed  ohonon  absorp- 
p  p  ap  ■* 


h  _  =  h _ )  for  the  ideal  pulse,  an  ex- 

m  ap  mp  7 

pression  may  be  obtained  from  (4.37a)  for  the  apparent  critical  field 


ap 


as  follows : 


4*  (Vap 

haP  2’>N(tap)  V 


(*>■50) 


The  h  which  appears  in  the  log  argument  has  only  a  second  order 

e/iect  on  the  actual  value  of  h  .  From  (4.50)  it  is  evident  that 

ap 

the  longer  the  build-up  time  t  of  the  phonon  instabilities,  the 

ap 

closer  h  approaches  h  .  Thus,  ideally,  the  critical  threshold 
ap  cr 

would  be  observed  with  an  incident  pump  pulse  of  infinite  length.  Since 

the  pulse  must,  practically,  have  a  finite  length,  then  it  follows  that 

the  lowest  apparent  critical  field  h  will  be  observed  at  the  end  of 

ap 

the  pulse--at  which  point  the  instabilities  have  had  the  longest  time  to 

grow.  Hence  the  build-up  time  t  is  equal  to  the  pulse  length  in  the 

ap 
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ideal  pulse  case.  It  is  also  noteworthy  that  the  smaller  the  observable 


change  in  the  phonon  power  absorption  P 


the  closer  h 


approaches 


p  '  ap 

the  critical  value.  In  practice,  the  critical  field  h  is  never  reached, 

vi 

and  thus  the  factor  h  is  always  greater  than  unity  for  observable 

phonon  instability  thresholds.  Nevertheless,  it  is  possible  to  extrapolate 
back  in  order  to  find 

cr 

field  amplitude  h  above  the  critical,  value. 


hcr  by  measuring  growth  rates  as  a  function  of 


2.  Actual  Pulse 

The  actual  pulse,  with  its  finite  decay  constant  t  ,  is  considered 
next.  The  phonon  absorption  in  this  case  is  obtained  from  (4.37®) 

with  t  finite.  Solving  (4.37a)  for  hffl  in  the  same  way  as  for  the 
ideal  pulse  case  yields  the  following  result: 


<P  ) 

c?.,  .ElSE..  , 

*>■<  V  V 


where  (P  )  is  again  the  observed  phonon  absorption  level  when  satur 
P 

ration  of  the  reflected  puap  pulse  firirt  appears,  and  i a  the  value 

of  hm  at  this  threshold  absorption  level  (see  Fig.  4.3). 

It  would  now  be  possible  to  write  the  apparent  critical  field,  h^ 
(Fig.  4.3)  in  terms  of  (4.51).  At  this  poipt,  however,  it  is  easiest 
to  work  with  the  field  amplitude  ii  ,  and  the  results  below  will  be 
given  in  terms  of  this  quantity. 

The  value  of  t  is  the  chief  difference  between  the  actual  end 
ap 

the  ideal  pulse  cases.  It  can  be  seen  by  an  examination  of  (4. 51)  and 

(4.50)  that  if  t  is  smaller  for  the  actual  pulse  case,  then  the 

ap 

'threshold  will  be  higher.  If  the  theory  for  the  actual  pulse  is  to  be 
useful  in  practice,  some  criterion  must  be  established  for  determining 
t  .  From  the  expression  (4.37a)  for  the  phonon  power  it  is  noted 
that  the  dominant  time -dependent  factor  is  the  exponential.  Initially, 


as  the  time  increases  the  exponential  grows  in  magnitude.  However,  at 
sufficiently  large  times  the  effect  of  the  pump  field  amplitude 


h'  *  bj1  “  VT) 


becomes  significant  and  the  absorbed  power  goes  through  a  maximum. 

The  effect  of  the  less  significant  time -dependent  factor  N  is  to  cause 
this  maximum  to  appear  sooner  than  it  otherwise  would. 

The  assumption  is  now  made  that  this  maximum  value  of  P^  is  equal 
to  (P  )  occurring  at  time  t  (Fig.  4.3).  The  time -dependent  portion 
of  (4.37a)  is 


2^t  |~-1  +  E2-  (l  - 

L  cr  '  'J 


(4.52) 


Differentiating  (4.52)  to  obtain  the  condition  for  a  maximum  value  of 

P  yields 
p 


(*.53) 


as  the  apparent  build-up  time  of  the  instability.  This  result  shows  that 
the  longer  the  pulse  time  constant  t  f  the  longer  is  the  time  t 

ap 

before  the  phonon  absorption  reaches  a  maximum — provided  the  pulse  length 

is  sufficiently  long.  If  t  approaches  infinity,  then  t  must  approach 

ap 

the  pulse  length,  as  in  the  ideal  pulse  case. 

An  assumption  is  now  made,  which  is  justified  in  practice  for  the 
range  of  the  parameters  employed,  that 


(4.54) 
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Using  this  approximation,  (4.53)  becomes 


(4.55) 


Substituting  (4.55)  into  (4.51)  and  solving  for  h ^/h^  yields  the 
following  result: 


(4.56a) 


where  the  definition  has  been  made  that 


l 


1 

- 4n 

2»F 


(P  )  h 

.  P.  »g  ■;» 


2n»(t#p)  y  h 


mp 


(l».56b) 


In  determining  h1J^/hcr  *  the  higher  order  effect  of  h >5/bcr  in  the 
argument  of  the  logarithm  may  be  neglected  to  first  order,  and,  if  neces¬ 
sary,  may  be  accounted  for  later  by  means  of  successive  approximation. 

From  the  above  analysis  for  the  actual  pulse  shape,  estimates  will 
now  be  made  for  the  ratios  h and  bap/bcr  as  vail  as  for  the 
build-up  time  tftp  .  The  estimated  value  for  the  phonon  (■*  1500) 


from  Chapter  III  will  be  used  here. 

t  «  309  psec  and  (Pj_  «  0.74  If 

P  *P 

(Chapter  Vi),  the  following  results 
(4.55): 

(experiment  i) 


The  pulse  length  was  120  psec. 


Taking  the  experimental  parameters, 

,  associated  with  experiment  I 
are  obtained  from  (4.56),  (4.77)  end 

1.202 

1.014  (4.57) 

104  psec  . 
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In  a  similar  way,  on  taking  the  experimental  parameters  t  *  124  jisec 

and  (P  )  «  0.0 6V  associated  with  experiment  II  (Chapter  VI )  we 

P  fiP 

obtain  the  following  results: 

(experiment  II ) 


h 

mp 


1.319 


=  1.036 

h 

cr 

t  =60  usee 
ap  H 


(^.58) 


In  this  case  the  pulse  length  is  90  nsec.  The  build-up  time,  t  , 

ap 

is  shorter  in  experiment  II  chiefly  because  the  pulse  decay  time  constant 

t  is  shorter.  Thei  more  sensitive  detection  used  in  experiment  II  also 

aided  in  reducing  t  .  The  amplitude  ratios  hn?/hcr  and  kap/hcr 

are  higher  for  experiment  II  because  the  shorter  the  pulse  decay  time 

constant  the  higher  must  be  the  power  level  (or  field  level)  to  produce 

a  given  time  average  power  level.  A  comparison  of  the  calculated  and 

observed  values  of  t  will  be  made  in  the  following  section  and  in 

ap 

Chapter  VI. 


E.  PHONON  ABSORPTION  EFFECTS 

The  expression  for  the  theoretical  phonon  absorption  as  a 

function  of  time  is  given  by  (4.37a).  Since  the  effects  of  phonon 
absorption  were  observed  in  greater  detail  in  experiment  II,  the 
theoretical  nature  of  the  phonon  absorption  for  the  parameters  of 
experiment  II  will  be  investigated  further.  Taking  the  parameter 
values  for  experiment  II  and  the  ratio  ti/h^  [given  above  by  (4.58)] 
for  which  the  minimum  detectable  phonon  power  absorption  takes  place, 
(4.37a)  may  be  plotted  as  a  function  of  time.  This  plot  is  shown  in 


-  ice 


Fig.  4.7.  The  peak  absorption  occurs  at  t  *  60  ^sec  ,  as  predicted 

ap 

by  (4.55)*  The  curve  is  not  quite  symmetrical  about  t  ■  t  ,  but  is, 
in  fact,  slightly  steeper  for  t  less  than  t  . 

Next,  the  effect  of  this  absorption  upon  the  reflected  pump  pulse 
will  be  predicted.  Combining  (4.2)  with  (4.44),  the  reflected  pump  wave¬ 
form  in  the  absence  of  Instability  is  given  by 


where 


* 


P 

rm 


('*•59) 


The  reflected  pump  waveform  (4-59}  Is  plotted  In  Fig.  4.8a  for  the  expert 
ment  II  value  of  T  and  for  an  arbitrary  power  level  below  threshold. 

At  threshold  power  levels,  (4.59)  la  incomplete  because  the  effect  of 
the  phonon  absorption  upon  the  total  reflected  power  must  be  considered. 
Then  the  total  reflected  power  (pr)^ot  becomes 


(P  ),  .  »  P  +  AP  ,  (4.60) 

'  r'tot  r  r  9 

where  is  given  by  (4.59),  and  the  change  in  reflected  power  AP^ 
is  found  from  (4.47)  and  (4.48)  to  be 


(4.61) 


for  the  overcoupled  case.  Since  AP  is  small,  it  will  be  assumed  that 

r 

jrj  does  not  change  significantly  during  the  growth  of  P^(t) 
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Lo'o 


FIG.  4.7 


Normalized  phonon  absorption  curves  Tor: 

(a)  <  1.319  ,  -S*  <  1.036 

h  h 

cr  cr 


M 

-2®  .  1.319 

,  -SB  .  1.036 

h 

h 

cr 

cr 

t  *  60  (isec 
ap 

h 

h  ^ 

(c) 

.  1.327 

,  -22  «  1.036 

h 

cr 

h 

cr 

P*P 


0.06  K 


p'ap 


t  *  6l.2  (INC 

ap 


FIG.  4. 8- -Normalized  reflected  power  vs  time  shoving  tbs  predicted  phonon 
absorption  for  the  parameter  of  experiment  II  and  *  1300 
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Assume  now  that  the  ratio  of  the  minimum  observed  change  in  reflected 
power  to  the  reflected  power  at  the  apparent  threshold  is  given  by 


|  *  0.1  .  (4.62) 

ap 

Recall  that  the  phonon  absorption  build-up  to  and  decay  from  the  apparent 
threshold  peak  was  plotted  in  Fig.  4.7  for  this  case.  Through  (4.62), 
(4.60),  and  Fig.  4.7/  the  effect  of  the  phonon  absorption  upon  the  re¬ 
flected  pump  pulse  is  determined,  and  the  result  is  plotted  in  Fig.  4.8b. 
The  reflected  pulses  for  cases  (a)  and  (b)  have  been  normalized  in  order 
to  show  best  the  deviation  in  the  pulse  waveform  due  to  the  phonon 
absorption . 

At  a  still  higher  ratio  of  AP  to  P  of 

r  r 

|  -  0.5  ,  (4.63) 

ap 


a  similar  result  obtains  and  is  shown  in  Fig.  4.8c.  The  curve  is  again 

normalized  as  above.  Note  the  build-up  time  t  is  slightly  longer 

ap 

in  this  case--in  accordance  with  (4.55)-  Values  of  the  parameters  for 
each  of  the  above  three  cases  are  shown  in  Fig.  4.8. 

The  above  theory  is  based  upon  the  assumption  that  |r|  *  constant. 
However,  the  fact  that  |r|  varies  by  40$  during  the  pulse  [as  in  case  (c)] 
is  not  enough  to  alter  significantly  the  nature  of  the  absorption  during 
the  pulse.  This  is  due  to  the  sharpness  of  the  absorption  peak.  If  the 
variation  in  |r|  becomes  much  larger,  however,  the  phonon  absorption 
peak  will  appear  broadened  due  to  two  effects.  The  first  follows  from 
(4.47),  where  the  effect  of  a  smaller  |r|  is  to  make  the  change  in 
reflected  power  smaller,  i.e.,  a  given  phonon  absorption  appears  smaller. 
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The  second  is  due  to  the  effect  of  a  smaller  |r|  in  decreasing 
in  (4.4la)  [see  also  (4.47)  and  (4.39) ]>  and  hepce  the  rf  magnetic  field 
amplitude.  The  latter  has  the  effect  of  reducing  the  growth  rate  of  the 
absorption  [through  reduction  of  h .^h  In  (4. 37a)  ] . 

The  phonon  absorption  shown  in  Fig.  4.8  has  been  observed  experi¬ 
mentally^  and  the  actual  waveforms  of  the  pump  pulse  will  be  discussed 
in  Chapter  VI.  The  build-up  times  observed  experimentally  were  con¬ 
siderably  shorter  than  those  calculated  [(4.57)  and  (4.58)].  This  may 
imply  that  the  phonon  is  actually  lower  than  originally  estimated. 

Following  a  process  almost  the  reverse  of  that  used  above,  a  value  of 
Qp  can  be  determined  which  will  correspond  to  the  observed  build-up 

time  (t  )  .  For  instance,  knowing  (t  )  _  ,  then  h  /h  is 

found  from  (4.55)  to  be 


Having  a  result  for 


the  result 


l  - 


(4.64) 


(4.65) 


follows  from  (4.56a).  The  value  of  may  them  be  determined  from 

the  expression  (4.56b)  for  -t  .  Taking  the  observed  build-up  time  of 
(t  )  «  32  usee  for  experiment  I  with  other  parameter  values  given 

previously,  the  value  of  Q  obtained  is 


(experiment  I) 


(4.66) 
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with 


Vh=r  *  1'055 

haj/her  ’  ' 

Taking  the  observed  build-up  time  of  (^ap^xp  *  ^8  jjisec  for  experiment 
II  with  previously  given  parameter  values,  the  Q  in  this  case  becomes 


(experiment  II ) 


(4.67) 


with 

*  1.127 

*  O.96I 


^mp^cr 

h  Jh 
ap  cr 


It  is  expected  on  the  basis  of  the  linear  theory  that  the  two  above 
values  of  should  be  equal.  However,  if  nonlinear  damping  is  present 

this  need  not  be  required.  In  such  a  case  the  above  Q^’s  nay  be  spoken 
of  as  "effective”  Q^s  .  The  nature  of  nonlinear  damping  is  not  known, 
but  the  effect  of  such  damping  would  most  likely  be  to  cause  the  in 

(4.37a)  to  decrease  with  increasing  phonon  absorption.  This  is  precisely 
the  relationship  which  exists  for  the  above  results  (4.66)  and  (4.67). 

A  further  substantiation  of  the  presence  of  nonlinear  daoping  Is  seen 
in  the  reflected  pump  waveforms  of  experiments  I  and  II  of  Chapter  VI, 
where  the  build-up  time  was  found  to  decrease  with  increasing  absorption 
levels.  Such  an  observation  is  contrary  to  the  result  (4.55)  fron  the 
linear  theory  with  Q.  constant,  where  the  reverse  is  predicted.  However, 

j? 

if  the  "effective"  is  allowed  to  decrease  with  increasing  phonon 
absorption,  then  the  linear  theory  [see  (4.37)1  will  also  predict  the 
observed  behavior. 
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F.  SUMMARY 


The  various  loss  mechanisms  of  a  ferrite  are  discussed  briefly. 
Interest  here  vas  in  operation  at  a  pump  frequency  and  dc  magnetic  field 
where  the  phonon  threshold  will  dominate  other  instability  thresholds. 

The  elastic  wave  equation  obtained  from  Chapter  III  for  vaVes  having 
the  lowest  threshold  was  modified  to  include  the  time  variation  of  the 
pumping  angle  arising  from  experimental  limitations.  '  An  approximate 
„  solution  to  this  equation  was  found  for  values  of  the  pumping  angle 
above  the  critical  threshold  angle.  A  calculation  of  the  power  being 
absorbed  by  the  unstable  elastic  modes  was  then  made  by  assuming  growth 
of  these  modes  from  thermal  energy  levels.  The  effective  number  of 
growing  modes  was  found  to  be  a  function  of  time  and  depends  upon  the 
value  of  the  pumping  angle  above  the  critical  angle,  as  well  as  upon  the 
ferrite  sample  dimensions,  phonon  Q  ,  and  effective  internal  dc  field. 
The  absorption  observable  experimentally  on  the  reflected  pump  pulse  was 
found  to  be  proportional  to  the  phonon  absorption  by  the  ferrite  and 
inversely  proportional  to  the  reflection  coefficient. 

The  optimum  conditions  required  for  observation  of  the  relative^ 
change  in  reflected  power  were  also  discussed.  The  apparent  critical 
magnetic  field  was  derived  for  the  ideal  and  actual  pulse  cases  in  terms 
of  the  phonon  absorption  level  which  is  first  observable  experimentally. 
The  build-up  time  corresponding  to  this  absorption  level  was  also  derived. 
For  experimental  and  ferrite  parameters  the  phonon  absorption  and  its 
effect  upon  the  reflected  pump  pulse  were  calculated  for  the  actual 
pulse.  The  theoretical  build-up  time  to  the  phonon  absorption  peak  was 
found  to  be  longer  than  that  observed  experimentally.  It  is  suggested 
that  nonlinear  damping  is  present. 


CHAPTER  V 


EXPERIMENTAL  TECHNIQUES  AND  MICROWAVE  PARAMETER  MEASUREMENTS 

A.  MICROWAVE  CIRCUITRY  FOR  THRESHOLD  MEASUREMENTS 

The  rf  magnetic  field  required  in  order  to  observe  the  phonon  thresh¬ 
old  may  be  estimated  from  the  theoretical  results  of  Chapter  III.  From 
Fig.  3*3>  a  typical  threshold  angle  may  be  taken  to  be 


Assuming  a  dc  magnetic  field  of'  HQ  «  4C0  Oe,  the  amplitude  of  the  critical 
rf  magnetic  field  will  be 

h.  4C  Oe  . 

cr 

r 

To  achieve  such  a  high  rf  field  at  the  frequency  of  interest  (l.32  kMc), 
a  conventional  metal  ..cavity  was  first  considered.  Large  incident  power 
levels  (~  100  kW)  would  be  required  to  achieve  the  desired  rf  field 
strength.  This  is  unfortunate  because  high  voltage  breakdown  within  the 
cavity  is  likely  to  occur  at  thes  power  levels.  The  conventional  metal 
cavity  at  1.32  kMc  also  has  the  disadvantage  of  having  a  small  filling 
factor  for  a  ferrite. 

Since  the  prospect  of  breakdown  was  to  be  avoided  and  a  large  filling 

factor  was  desired,  other  methods  of  resonating  a  microwave  electromagnetic 

field  were  investigated.  Such  a  method  was  found  in  the  dielectric 
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reaonator.  '  y  A  dielectric  resonator  having  a  high  dielectric  constant 
was  considered.  In  such  a  case  the  possibility  of  breakdown  is  practi¬ 
cally  eliminated.  Also,  a  large  filling  factor  can  be  obtained  because 
the  resonator  volume  is  small.  Furthermore,  large  rf  fields  can  be 


110  - 


obtained  with  considerably  lower  incident  power  levels  than  would  be 
required  for  a  metal  cavity. 

A  dielectric  resonator  was  fabricated  from  single  crystal  strontium 

titanate.  Strontium  titanate  possesses  a  very  high  isotropic  dielectric 

constant  (s'  *  277),  a  low  loss  tangent,  and  is  commercially  available 

in  large  single  crystal  boules.  With  such  a  high  dielectric  constant, 

x 

the  resonator  dielectric  volume  at  1.32  kMc  was  only  1.28  cm  ,  and  a 
large  ferrite  filling  factor  was  obtained.  A  resonator  system  Q  of 
over  1000  was  also  obtainable.  Thus,  direct  observation  of  absorbed 
power  by  the  ferrite  was  easily  possible. 

In  order  to  achieve  the  large  rf  magnetic  fields  with  the  dielectric 
resonator,  it  was  still  necessary  to  employ  a  moderately  high  microwave 
power  source--one  having  a  maximum  power  output  of  about  50  W  peak. 
Thermal  limitations  of  the  dielectric  resonator  and  ferrite  dictated 
pulsed  rather  than  cw  operation.  The  microwave  power  source  was  pro¬ 
vided  by  a  50C  mW  signal  generator  (General  Radio  GR-1218-A)  in  con¬ 
junction  with  an  L-band  klystron  amplifier  (Raytheon  SAL-8l),  shown  in 
Fig.  5*1*  The  klystron  beam  was  pulsed  by  the  high  voltage  pulser 
modulator  (Levinthal  79-M).  The  modulator  contained  the  power  supply 
which  charged  up  a  15  kV  1  nfd  capacitor.  The  pulsed  discharge  from 
the  capacitor  furnished  the  beam  current  for  the  klystron.  The  pulse 
length,  or  length  of  each  discharge  time,  was  about  100  lisec.  This  was 
long  enough  to  enable  observation  of  the  predicted  phonon  absorption 
peak.  The  pulse  length  was  determined  by  a  multivibrator  circuit  which 
controlled  the  capacitor  discharge  through  tt  vacuum  tube  circuit.  The 
switching  characteristics  were  limited  to  a  large  extent  by  the  quality 
of  the  transformer  employed  at  the  output  of  the  multivibrator  circuit. 
For  transformers  available,  a  compromise  had  to  be  made  between  fast 
rise  and  fall  times  and  minimum  sag  for  the  100  usee  pulse,  me  effect 
of  the  transformer  sag  was,  of  course,  to  enhance  the  sag  in  the  capaci¬ 
tor  discharge  and  hence  in  the  amplified  rf  output.  The  pulae  repetition 
rate  was  controlled  by  the  pulse  trigger  external  to  the  modulator 
It  was  necessary  to  make  the  pulse  repetition  rate  as  low  as  possible 
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in  order  to  minimize  the  sag  in  the  capacitor  discharge  current.  Yet 
this  rate  had  to  be  sufficiently  high  to  give  a  continuous  display  of 
the  power  pulse  waveforms  on  the  long  persistence  screen  of  the  oscillo¬ 
scope.  This  determined  the  pulse  repetition  rate  of  10  cps. 

An  attenuator  pad  served  to  improve  the  impedance  match  of  the 
klystron  to  its  load.  After  the  attenuator  pad,  a  dual -directional 
coupler  was  inserted  for  the  purpose  of  detecting  incident  and  reflected 
power  from  the  dielectric  resonator.  The  output  of  the  dual -directional 
coupler  was  fed  to  a  precision  variable  attenuator  (Hewlett-Packard  HP-J^-A). 
The  crystal  at  the  output  of  this  attenuator  detected  the  power  level. 

All  power  measurements  were  made  at  the  convenient  minimum  signal  level 
of  2  cm  deflection  (0. IV)  on  the  oscilloscope.  Thus,  for  any  given 
power  input  to  the  load,  there  was  a  corresponding  setting  of  the  pre¬ 
cision  variable  attenuator  which  gave  a  2  cm  deflection  on  the  oscillo¬ 
scope.  The  variable  attenuator  was  originally  calibrated  with  a  thermis¬ 
tor  as  &  load.  The  thermistor  was  used  in  conjunction  with  a  power 
meter  to  determine  the  average  power  level  for  each  attenuator  setting 
over  the  desired  power  range.  Knowing  the  pulse  length,  pulse  repetition 
rate,  and  the  average  power,  the  peak  power  level  was  obtained  and  the 
attenuator  calibrated  accordingly. 

The  circuitry  connections  for  the  above  apparatus  were  made  with 
7/8  in.  and  type-N  coax*  An  adapter  to  the  standard  L-band  rectangular 
waveguide  was  employed  at  the  output  of  the  directional  coupler.  The 
waveguide  provided  a  means  of  coupling  power  to  the  dielectric  resonator. 
Sufficient  spacing  between  the  dielectric  resonator  and  the  adapter  was 
allowed  to  avoid  any  higher  order  electromagnetic  fringe  fields  arising 
from  the  transition. 

The  dc  magnetic  field  was  first  applied  across  an  air  gap  the  width 
of  the  L-band  waveguide.  The  field  was  obtained  from  an  improvised  Mag¬ 
net  fabricated  from  two  large  toroidal  field  coils  and  iron  bars.  The 
latter  provided  the  necessary  low  reluctance  flux  path.  With  the  field 
coils  located  near  the  pole  faces  (for  maxia\aa  gap  flux),  a  sufficiently 
large  and  fairly  uniform  dc  magnetic  field  was  attained.  The  Magnetic 
field  was  first  calibrated  by  a  gaussmeter  of  the  Hall  effect  type— 
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which  was  later  calibrated  by  a  nuclear  magnetic  resonance  (NMR)  probe. 

The  first  set  of  experimental  threshold  measurements  (hereinafter 
referred  to  as  experiment  i)  were  made  by  Rising  this  magnet  and  the 
setup  of  Fig.  $.1. 

Another  set  of  threshold  measurements  (hereinafter  referred  to  as 
experiment  II)  was  made  with  a  large  commercial  magnet  (Varian  12  in.). 

In  this  case  the  height  of  the  L-band  waveguide  had  to  be  reduced  to 
3  in.  to  fit  inside  the  magnet  air  gap.  For  both  magnets  a  constant 
current  regulated  power  supply  served  to  maintain  a  constant  magnetic 
field  with  a  minimum  of  drift  due  to  coil  heating.  With  the  Varian 
magnet  it  was  possible  to  obtain  a  broader  range  of  magnetic  field  and 
a  greater  magnetic  field  uniformity  at  the  sample,  and  an  NMR  probe 
could  be  used  directly  to  calibrate  the  field.  Another  improvement  in 
experiment  II  included  a  circulator,  which  increased  the  sensitivity 
of  the  critical  field  measurements. 

B.  DIELECTRIC  RESONATOR  SYSTEM 

The  dielectric  resonator  system  will  now  be  described  in  greater 
detail.  The  experimental  arrangement  is  shown  in  Fig.  5.2.  The  di¬ 
mensions  of  the  strontium  titanate  dielectric  resonator  were 
0.430  in.  X  0.431  in.  X  0.423  in.  It  has  been  mentioned  already  that 
the  properties  of  strontium  titanate  utilized  were  its  low  loss  tangent 
and  its  very  high  and  essentially  isotropic  dielectric  constant  (e'  *  277)* 
The  rod  shown  supporting  the  resonator  also  had  a  low  loss  tangent,  but 
it  had  a  low  dielectric  constant.  The  dielectric  resonator  was  located 
in  the  center  of  the  L-band  waveguide  cross-section.  A  hole  (0.06l  in.) 
was  cut  in  the  resonator  as  shown  and  a  ferrite  disk  (0.054  in.  diam 
x  O.OO76  in.)  supported  by  a  glass  rod  was  inserted.  This  resonator 
operated  in  the  dominant  H  mode  at  I.32  kMc.  In  the  field  configuration 
of  this  mode  the  maxinum  H  field  occurs  at  the  center  of  the  resonator 
and  is  polarized  in  the  direction' of  the  glass  rod  (see  Fig.  5*3)*  Since 
the  E  field  is  weakest  where  the  H  field  is  strongest,  the  hole 
distorts  the  field  pattern  very  little.  The  coupling  of  the  resonator 
H  mode  to  the  standing  waveguide  TE^q  mcx*e  was  opiiioized  by  position¬ 
ing  the  movable  short. 
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MOVABLE  SHORT- 


FERRITE 


FREQUENCY 

TUNING 

PLUNGER 


/ 

!  / 
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GLASS 

ROD 


STRONTIUM  TITANATE 
DIELECTRIC  RESONATOR 


TO  HIGH  POWER 
KLYSTRON  SOURCE 


H0»dc  MAGNETIC 
FIELD 


h9rf  MAGNETIC 
FIELD 


FIG.  5. 2- -Setup  for  ferrite  transverse  punping  experiment  using 
strontium  titanate  dielectric  resonator. 
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Electric  field  lines 
Magnetic  field  lines 


FIG.  5»3 — Dielectric  resonator  shown  with  coordinate  system  and 
a  few  of  the  electric  and  magnetic  field  lines  for 
the  dominant  H-mode . 
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The  properties  of  the  ferrite  and  of  the  dielectric  resonator  were 
found  to  be  significantly  t emperatur e -dependent .  The  ambient  temperatures 
of  the  resonator  and  ferrite  sample  were  determined  by  mounting  a  thermo* 
couple  on  the  post  supporting  the  dielectric  resonator.  The  dielectric 
constant  of  the  resonator  had  'the  greatest  sensitivity  to  temperature 
changes  as  far  as  the  actual  threshold  measurements  were  concerned, 
since  the  resonator  frequency  was  thereby  affected.  A  frequency  pertur¬ 
bation  plunger  was  employed  to  compensate  for  this  effect.  The  tempera¬ 
ture  coefficient  of  the  resonant  frequency  was  about  2.5  Mc/°C.  The 
plunger  was  capable  of  changing  the  resonant  frequency  over  a  range  of 
about  25  Me.  With  the  three  cavity  klystron  adjusted  for  1.J2  kMc,  it 
was  then  possible  to  resonate  the  dielectric  resonator  over  an  ambient 
temperature  range  of  about  10°C.  It  was  found  desirable  to  perform  ex¬ 
periments  at  ambient  temperatures  where  the  plunger  could  be  maintained 
as  far  as  possible  from  the  resonator  and  still  provide  adequate  frequency 
tuning  during  the  data  runs.  This  minimized  distortion  of  the  resonator 
rf  field,  auid  reduced  the  dependence  of  the  unloaded  upon  the 

plunger  position. 

C.  HEATING  OF  THE  DIELECTRIC  RESONATOR  AND  FERRITE 

Heating  of  the  dielectric  resonator  and  ferrite  due  to  microwave 
power  absorption  will  now  be  investigated.  First  consider  the  dielectric 
resonator  without  the  ferrite.  The  average  microwave  power,  pft  , 
absorbed  by  the  dieleitric  resohator  systeah  id  “dissipated  through  di¬ 
electric  loss  in  the  strontium  titanate  resonator,  the  dielectric  sip- 
porting  rod,  the  Ducco  cement  contact  between  the  rod  and  resonator, 
finger  contacts  of  the  movable  short,  and  waveguide  walls.  The  power, 

Pd  ,  actually  dissipated  in  the  vicinity  of  the  resonator  may  be 
given  by 


(5.1) 
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Assuming  -that  the  support  rod  and  the  cement  contact  have  a  somewhat 
lower  loss  than  the  dielectric  resonator,  it  is  estimated  that  ^  *  0.5  . 

The  energy  dissipated  within  the  resonator  volume  is  related  to 
the  temperature  increase  £1  by 

Ed  ■  Wi*  '  (5-2) 


where  c^  is  the  specific  heat,  is  the  density,  and  is  the 

volume  of  the  dielectric  resonator.  The  factor  f .  is  a  factor  greater 

d 

than  unity  and  is  used  to  account  for  other  energy  losses  in  the  vicinity 
of  the  resonator  due  to  thermal  conduction,  convection,  and  radiation. 
Equating  (5*2)  to  the  energy  input  P ^At  yields  the  following  for  R^  , 
Ihe  rate  of  change  of  resonator  temperature  with  time: 


dT 

dt 


fdVdvd 


(5-3) 


Taking  the  typical  average  power  level  p  »  10  mW  (corresponding  to 

»  a 

1CW  peak  power  and  10“J  duty  cyclte),  and  using  the  values  of  cd  and 

p  for  strontium  titanate  given  in  Appendix  A,  and  the  resonator 

volume  =  1.28  cnr  ,  R .  is  estimated  to  be 

d  d 


R,(10  mW) 
d 


JC/min 


(5-*0 


where  f .  has  yet  to  be  evaluated. 

CL 

The  temperature  dependence  of  the  strontium  titanate  dielectric 
constant  strongly  affected  the  resonant  frequency  of  the  dielectric 
resonator.  From  the  value  of  the  temperature  coefficient  of  the  di¬ 
electric  constant  given  in  Appendix  A  it  can  be  shown  that  the  tempera¬ 
ture  coefficient  of  the  resonant  frequency  is  about  2.5  Mc/°C.  This 
dependence  was  observed  experimentally  as  described  in  the  previous 
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section.  The  time  rate  of  change  of  resonator  frequency  at  the  10  aM 

average  power  level  was  also  observed  experimentally .  From  these  ob¬ 
servations  the  time  rate  of  change  of  resonator  temperature  was  found 
to  be 


R  (10  mW)  *  0.07  °C/min  .  (5*5) 

exp 


By  comparison  with  (5*4)  It  is  estimated  that  f^  *  5*5  • 

Equation  (4.53)  expresses  the  relation  between  the  peak  pulse  ab¬ 
sorbed  power  P&  and  the  peak  pulse  stored  energy  V  in  the  dielectric 

resonator  system.  For  the  average  :  absorbed  power  p  and  the  average 

£ 

stored  energy  U  ,  the  relation  is 


(5.6) 


This  expression  applies,  of  course,  to  the  resonator  system  with  or 

without  the  ferrite  present.  The  effect  of  the  ferrite  upon  U  is  very 

small.  The  ratio  of  the  ferrite  volume  to  the  yolume  of  the  dielectric 

-4 

resonator  Is  of  the  order  of  10  .  Even  if  the  ferrite  has  a  relative 

permeability  as  large  as  100,  the  effect  of  the  ferrite  qpon  the  stored 
energy  u  will  be  negligible  (~  Vft )  .  The  effect  of  the  ferrite  vpon 
the  is  about  an  order  of  magnitude  larger.  In  fact,  a  typical 

change  in  observed  experimentally  was  about  10J&. 

Considering  only  the  effect  upon  ,  the  average  power  absorbed 

by  the  ferrite  is  found  to  be 


(5-7) 


Thus  the  ratio  of  the  total  average  power  absorbed  by  the  ferrite  to 
the  average  power  absorbed  by  the  resonator  system  is  given  approximately 
by  the  relative  change  in  . 
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The  time  rate  of  change  of  temperature  of  the  ferrite  is 


Py 

Ry  *  - - -  ,  (5.8) 

fy  Cy  Py  Vj, 

where  the  parameters  are  counterparts  of  those  used  in  (5*3  )•  The  factor 

f_  which  accounts  for  conduction,  convection,  and  radiation  losses  of 
i 

the  ferrite  is  expected  to  be  larger  than  the  factor  for  the  resonator. 
This  is  because  the  broad  face  of  the  ferrite  is  mounted  on  a  glass  rod 
which  provides  a  fairly  good  heat  sink — thereby  promoting  large  con¬ 
duction  losses.  The  dielectric  support  for  the  resonator  conducts  heat 
away  from  about  one-sixth  of  the  resonator  surface  area,  whereas  the 
gLass  rod  conducts  heat  away  from  about  one-half  of  the  ferrite  area. 

Consequently,  it  might  be  estimated  that  f_  *  •  Thus,  with  a 

-4  3  r  d 

ferrite  volume  Vy  ■  2.86  x  10  cm  and  the  values  given  in  Appendix 
A  for  the  other  quantities,  (5*8)  becomes 

Ry(lO  mW)  *  0.38  °C/sec  .  (5.9) 

Therefore,  if  it  takes  10  sec  to  determine  the  data  for  one  experimental 
point,  then  the  temperature  increase  of  the  ferrite  would  be  about 
AT  *  3*8°C  .  It  is  clear  that  measurements  must  be  taken  quickly  at 
these  high  power  levels.  Ifce  threshold  level  was  usually  determined 
as  quickly  as  possible,  but  never  in  much  less  than  about  10  sec.  After 
each  threshold  measurement  the  incident  power  was  reduced  and  the  sample 
was  allowed  to  cool  before  the  next  measurement  was  taken.  It  was 
assumed  that  the  ferrite  reached  the  ambient  temperature  during  this 
cooling  period. 

At  the  lowest  power  levels  of  about  P  *  1  W  peak,  the  estimated 

ft 

ferrite  temperature  increase  for  a  10  sec  measurement  time  would  be 
AT  *  0.38°C  .  Here,  then,  the  actual  temperature  of  the  ferrite  differs 
little  from  the  equilibrium  value.  This  equilibrium  temperature  was 
determined  experimentally  by  attaching  a  thermocouple  to  the  dielectric 
rod  supporting  the  resonator. 
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With  a  pulse  repetition  rate  of  10  cps,  the  average  temperature 
Increase  of  the  ferrite  during  each  pulse  was  a  factor  of  10  smaller 
than  the  average  temperature  increase  of  the  ferrite  in  one  second* 
Thus,  heating  effects  during  the  pulse  were  negligible  at  the  average 
power  levels  considered. 


D.  MAGNETIC  FIELD  DISTRIBUTION  IN  TEN  DIELECTRIC  RESONATOR 


The  basic  relation  between  the  stored  energy  U  In  the  rf  magnetic 
field  of  the  resonator  and  the  peak  pulse  power  absorbed  by  the 
system  is  given  by  Eq.  (4.4o).  In  this  section,  U  will  be  obtained 
in  terms  of  the  maximum  rf  magnetic  field  amplitude  hQ  * 

The  rectangular  parallelepiped'  dielectric  reaonatpr,  together  .with 
the  coordinate  system  used  to  describe  the  electric  and  magnetic  fields 
for  the  dominant  H  mode,  is  shown  in  Fig.  5*3*  The  electric  field  Is 
strongest  near  the  resonator  faces  parallel  to  the  z-axis,  while  the 
magnetic  field  is  strongest  in  the  center  of  the  resonator  and  along 
the  z  direction.  A  few  of  the  strongest  field  amplitude  lines  are 
sketched  in  Fig.  5»3* 

Open  circuit  boundary  conditions  (OCB)  require  that  the  electric 
field  be  parallel  to  and  the  magnetic  field  normal  to  the  boundary. 

For  the  dominant  H  mode  described  above  it  is  found  that  at  ths  faces 
of  the  dielectric  resonator  parallel  to  the  z-axis,  the  OCB  conditions 
are  nearly  satisfied,  since  the  relative  dielectric  constant  is 
much  greater  than  unity.  The  advantage  of  recognizing  OCB  conditions 
is  that  the  theory  is  simplified  by  terminating  ths  electric  and  mag¬ 
netic  fields  at  the  air -dielectric  boundary.  For  ths  faces  of  ths 

resonator  normal  to  the  z-axis,  the  usual  bgttndary  conditions  for 

20 

electric  and  magnetic  fields  are  used.  Using  Tee's  results  for  ths 
special  case  of  a  rectangular  parallelepiped  with  a  square  cross  section 
(a»b )  ,  the  electric  fields  based  upon  the  above  assumptions  are  given 
by 


-^0^0 

2x 


xx  xy 

sin  —  cos  —  A(z) 

a  b 


(5.10a) 
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jo»infthn  xx  xy 

— . . .  cos  —  sin  —  A(z) 

2x  a  b 


(5.10b) 


where  h^  is  defined  as  the  maximum  value  of  the  magnetic  field  (at 
the  center  of  the  resonator).  Inside  the  resonator  'A(z)  is  defined  bj 


A(z)  *  cos  £  (z  -  j)  , 


(5.10c) 


and  outside  by 


A(z)  *  Aq  e 


-t0'*  - 1 


(5.104) 


where  the  quantity  A^  is  obtained  by  matching  boundary  conditions 
at  z  »  0  or  c  .  The  quantities  £  and  £Q  are  defined  approximately 


by  the  relations 


2  2 


(5.11) 


The  magnetic  field  components  may  easily  be  determined  by  substituting 
(5>10)  into  Maxwell's  equations.  Also  from  Maxwell's  equations,  the 
following  expression  for  the  resonant  frequency  is  derived: 


|V(ifl 


(5.12) 


It  was  by  inserting  experimental  values  for  a  *  b,  c  and  the  resonant 

frequency  f  into  (5*12)  that  the  experimental  value  for  strontium 

titanate  of  € '  *  277  ,  given  in  Appendix  A,  was  obtained.  This  result 

23  go 

agrees  favorably  with  values  measured  by  others.  *  y 

The  total  stored  energy  in  the  resonator  system  is  most  simply 
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given  by  the  integral  of  the  electric  energy  density  when  the  latter  la 
a  maxi mum  in  time,  for  example 

U  »  g-  J  e  |g|2  dxdydz  .  (5»13) 

apace 


It  is  now  assumed  that  the  OCB  conditions  at  the  resonator  faces  parallel 
to  the  z-axis,  which  apply  approximately  to  the  fields  within  the 
resonator,  apply  equally  well  outside  of  the  resonator.  With  this 
assianption,  the  result  of  taking  the  integral  (3*13)  over  all  space  is 


0  *  > 


(5.1^) 


where. 


for  a  « 


f  - 


b  the  factor  f  is  given  by 


(5.14b) 


By  solving  (3*11)  and  £c  and  £Qc  ,  using  the  resonator  dimensions 
given  previously,  the  value  of  f  obtained  is  f  *  0.0 6  •  It  is 
interesting  that  this  value  of  f  is  comparable  to  the  value  of  unity 
which  would  be  obtained  if  there  were  no  z -dependence  of  the  fields 
and  if  all  of  the  energy  were  contained  within  the  resonator  volume. 

In  deriving  the  shove  results,  OCB  conditions  were  assumed  at 
four  faces  of  the  rectangular  parallelepiped  dielectric  resonator  as 
well  as  at  the  extension  of  these  faces  external  to  the  resonator.  For 
strontium  titanate^  the  relative  dielectric  constant  is  high  (c '  *  277); 
and  the  OCB  conditions  are  well  satisfied  within  the  resonator.  Homing 
it  is  not  clear  how  good  the  OCB  assumption  is  outside  the  resonator. 
This  problem  could,  of  course,  be  overcome  by  solving  for  fields  satis¬ 
fying  the  exact  boundary  conditions,  but  this  would  be  Impossible  for 
the  rectangular  parallelepiped,  except  by  numerical  methods. 
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There  is,  fortunately,  another  alternative  hy  which  the  validity 

of  the  approximate  solution  for  the  rectangular  resonator  may  be  checked. 

The  spherical  dielectric  resonator  operating  in  the  dominant  H  mode  has 

electric  and  magnetic  field  configurations  quite  similar  to  those  of 

the  rectangular  parallelepiped  operating  in  the  same  mode.  What  makes 

the  spherical  dielectric  resonator  of  special  interest  is  that  the 

electric  and  magnetic  fields  can  be  determined  exactly,  and  therefore 

an  exact  expression  for  the  energy  storage  can  be  derived.  Using  the 
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energy  relations  given  by  Yee  7  for  the  dominant  mode  of  the  di¬ 
electric  sphere  (where  the  exact  resonance  condition  for  the  sphere  has 
been  employed),  the  factor  f  defined  by  (5.14a)  is  found  to  be  f  =  O.85  • 
Thus  a  dielectric  sphere  having  the  same  volume  as  the  rectangular 

parallelepiped,  but  necessarily  a  slightly  higher  frequency  (^),  has 
about  the  same  energy  storage.  This  result  makes  more  reasonable  the 
original  assumption  of  OCB  conditions  for  the  rectangular  parallelepiped 
dielectric  resonator,  and  in  particular  the  extension  of  these  conditions 
external  to  the  resonator. 

On  combining  (4.53)  and  (5«l4a),  the  expression  for  h^  becomes 


(5.15) 


By  inserting  the  aforementioned  parameter  values,  then  (5*15)  becomes 

h0  -  0.352  (PaiJ0)l/2  •  (5.16) 

The  ferrite  is  located  in  the  center  of  the  dielectric  resonator 

in  the  actual  experiment,  and  therefore  it  is  worthwhile  to  investigate 

the  variation  of  the  magnetic  field  near  this  point.  Figure  5.4  shows 

the  variation  of  the  dominant  magnetic  field  component  h  as  a  function 

z 

of  z  for  x  *  a/2  and  y  =  b/2  .  The  magnetic  field  h  is  a 

z 

maximum  at  the  center  of  the  resonator  and  has  the  same  kind  of  z -variation 


as  do  the  electric  fields  given  by  (5.10).  The  importance  of 
properly  centering  the  ferrite  sample  in  order  to  achieve  maximum  field 
strength  and  uniformity  is  obvious.  The  arrow  indicates  how  far  the 
edge  of  the  ferrite  extends  when  placed  at  the  resonator  center.  The 
magnetic  field  variation  over  the  ferrite  sample  is  about  1$.  Thus  h^ 
may  still  be  taken  approximately  as  the  average  magnetic  field  strength 
over  the  ferrite  sample  volume.  The  transverse  magnetic  field  at  the 
sample  edge  is  less  than  1$  of  the  maximum  field  strength.  Hence,  as 
far  as  field  uniformity  is  concerned,  the  maximum  deviation  of  the  mag¬ 
netic  field  lines  from  the  z  direction  is  somewhat  less  than  about  0.5°  - 

The  presence  of  the  frequency  tuning  plunger  causes  the  field  to 
distort  somewhat  from  the  configuration  described  above.  In  particular, 
it  causes  an  increase  in  the  lateral  (x-y)  components  of  magnetic  field. 
The  diameter  of  the  frequency  plunger  was  almost  as  large  as  the  resonator 
dimension  a  .  The  position  of  the  plunger  was  symmetrically  located 
parallel  to  the  z-axis  of  the  resonator.  The  distance  of  closest  ap¬ 
proach  for  the  plunger  of  0. 06  in.  is  indicated  in  Fig.  5.4.  The  dis¬ 
tance  corresponding  to  two  successive  turns  of  the  plunger  screw  is  also 
indicated.  The  resonator  was  operated  with  P  >  2  to  minimize  dis¬ 
tortion  of  the  field  while  still  being  close  enough  to  have  an  appreciable 
effect  upon  the  resonator  frequency. 

In  the  exact  solution  of  the  fields  for  the  spherical  dielectric 

2Q 

resonator  operating  in  the  mode,  it  has  been  shown  that  radiation 

of  power  from  the  resonator  takes  place.  This  will  also  occur  in  the 
rectangular  resonator.  As  the  plunger  is  moved  closer  to  the  resonator, 
a  greater  portion  of  the  total  solid  angle  is  shielded  against  radiatidn 
by  the  plunger.  Hence  it  is  expected  that  the  radiation  will  be  reduced 
as  the  plunger  becomes  closer  to  the  resonator.  Such  an  effect  was  ob¬ 
served  experimentally  by  an  increase  in  the  Qq  of  the  resonant  system. 
The  increased  by  about  the  same  percentage  whether  the  ferrite  wa3 

present  or  not. 


Hi 
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£.  MEASUREMENT  OF  THE  UNLOADED  ^ 

In  order  to  accurately  evaluate  the  rf  magnetic  field  strength  at 
the  ferrite  during  the  threshold  experiments  it  is  necessary  to  know 
the  unloaded  Qq  of  the  ferrite -loaded  resonator  and  the  variation  of 
Qq  with  Hdc  .  The  Q  was  measured  at  low  power  levels --far  helow 
that  power  level  at  which  the  instability  thresholds  occurred.  The 
low  power  level  Q's  were  adequate,  however,  since  the  effective  re¬ 
sistance  (or  imaginary  part  of  the  magnetic  susceptibility)  of  the  fer¬ 
rite  does  not  change  until  loss  mechanisms  due  to  growth  of  elastic 
or  spin  wave  instabilities  are  introduced.  Besides,  it  would  be  much 
more  difficult  to  measure  the  precise  value  of  the  Q  at  or  near  the 
high  threshold  power  levels  because  average  pulse  heating  effects  would 
be  significant,  particularly  if  the  half -power  method  were  employed. 

If  the  sweep  frequency  method  were  employed,  a  wide  bandwidth  high  power 
amplifier  would  be  required. 

1.  Basic  Relations 

The  relationship  between  the  loaded  ,  the  unloaded  , 
and  the  external  Q  ,  of  a  resonant  system  is  given  by^° 


1 


The  external  Q>ex^  is  related  to 


«o  by 

% 

P 


(5-17) 


(5.18) 


where  P  Is  the  coupling  coefficient.  Thus  the  unloaded  nap  he 
given  in  terns  of  4.  by 


«0  *  (1  +  P)  ^  .  (5.19) 
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The  coupling  coefficient  0  for  an  overcoupled  system  is  given  by 


0  = 


0 


and  for  an  under coup led  system  by 


0 


(5.20a) 


(5.20b) 


where  r^  is  the  VSWR  on  resonance  and  the  VSWR  is  related  to  the 
reflection  coefficient  T  by  the  following: 


r 


i  +  r 
i  -  r 


(5.21) 


The  external  Qgxt  given  by  (5*18)  is  a  quantity  dependent  upon 
the  coupling  between  the  microwave  power  rource  and  the  resonant  system. 
For  a  metal  cavity  this  coupling  is  usually  through  an  iris  or  a  probe. 
For  the  dielectric  resonator  system  the  resonator  is  placed  inaide  the 
waveguide  and  the  coupling  is  adjusted  by  means  of  a  movable  short. 

2.  Q-Circle  Method 

Initially,  an  attempt  was  made  to  measure  the  unloaded  Q0  by 
the  Q-circle  method  using  instruments  designed  for  a  1  kc  square  wave 
modulated  signal.  However,  the  time  involved  In  measuring  the  VSWR 
and  locations  of  the  VSWR  minima  was  sufficiently  long,  in,  spite  of  the 
low  power  level  (~  1  mW)  of  the  signal  generator,  that  dielectric  heating 
effects  distorted  the  results.  The  dielectric  heating  effects  observed 
were  due  to  the  microwave  power  absorption  of  the  resonator,  as  well 
as  to  the  general  ambient  temperature  increase  between  measurements. 

These  heating  effects  were  compensated  for  to  some  degree  by  taking 
each  set  of  readings  in  an  alternate  sequence  about  the  resonant 
frequency.  Nevertheless,  it  was  decide’  that  other,  more  precise  methods 
of  determing  the  Q  would  be  used  in  order  to  circumvent  these  heating 
effects. 
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E.  MEASUREMENT  OF  THE  UNLOADED  ^ 

In  order  to  accurately  evaluate  the  rf  magnetic  field  strength  at 
the  ferrite  during  the  threshold  experiments  it  is  necessary  to  know 
the  unloaded  of  the  ferrite-loaded  resonator  and  the  variation  of 
Q0  with  Hdc  .  The  Q  was  measured  at  low  power  levels --far  below 
that  power  level  at  which  the  instability  thresholds  occurred.  The 
low  power  level  Q's  were  adequate,  however,  since  the  effective  re¬ 
sistance  (or  imaginary  part  of  the  magnetic  susceptibility)  of  the  fer¬ 
rite  does  not  change  until  loss  mechanisms  due  to  growth  of  elastic 
or  spin  wave  instabilities  are  introduced.  Besides,  it  would  be  much 
more  difficult  to  measure  the  precise  value  of  the  Q  at  or  near  the 
high  threshold  power  levels  because  average  pulse  heating  effects  would 
be  significant,  particularly  if  the  half -power  method  were  employed. 

If  the  sweep  frequency  method  were  employed,  a  wide  bandwidth  high  power 
amplifier  would  be  required. 

1.  Basic  Relations 

The  relationship  between  the  loaded  ^  ,  the  unloaded  , 
and  the  external  of  a  resonant  system  is  given  by^° 


The  external  **  related  to  by 


(5-17) 


(5.18) 


where  0  is  the  coupling  coefficient.  Thus  the  unloaded  may  he 
given  in  terms  of  by 

«0  -  (1  +  P)  ^  .  (5.19) 
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Hie  coupling  coefficient  £  for  an  overcoupled  system  is  given  by 


0  = 


(5.20a) 


and  for  an  under coup led  system  by 


0 


(5.2Cfc) 


where  rQ  is  the  VSWR  on  resonance  and  the  VSWR  is  related  to  the 
reflection  coefficient  r  by  the  following: 


r 


1  + 

[rL 

1  - 

m 

(5.21) 


The  external  Qext  given  by  (5 .  l8 )  is  a  quantity  dependent  upon 
the  coupling  between  the  microwave  power  source  and  the  resonant  system. 
For  a  metal  cavity  this  coupling  is  usually  through  an  iris  or  a  probe. 
For  the  dielectric  resonator  system  the  resonator  is  placed  inside  the 
waveguide  and  the  coupling  is  adjusted  by  means  of  a  movable  short. 

2.  (^-Circle  Method 

Initially,  an  attempt  was  made  to  measure  the  unloaded  Qq  by 
the  Q-circle  method  using  instruments  designed  for  a  1  kc  square  wave 
modulated  signal.  However,  the  time  involved  in  measuring  the  VSWR 
and  locations  of  the  VSWR  minima  was  sufficiently  long,  in,  spite  of  the 
low  power  level  (~  1  mW)  of  the  signal  generator,  that  dielectric  heating 
effects  distorted  the  results.  The  dielectric  heating  effects  observed 
were  due  to  the  microwave  power  absorption  of  the  resonator,  as  well 
as  to  the  general  ambient  temperature  increase  between  measurements. 

These  heating  effects  were  compensated  for  to  some  degree  by  taking 
each  set  of  readings  in  an  alternate  sequence  about  the  resonant 
frequency.  Nevertheless,  it  was  decided  that  other,  more  precise  methods 
of  determing  the  Q  would  be  used  in  order  to  circumvent  these  heating 
effects. 
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3.  Sweep  Frequency  Method 

The  sweep  frequency  method  was  employed  to  measure  the  Q  for 
experiment  I.  A  block  diagram  of  the  experimental  setup  is  shown  in 
Fig.  5*5*  A  mechanical  sweep  generator  was  employed  to  drive  the  fre¬ 
quency  dial  of  the  L-band  signal  generator.  A  sweep  frequency ‘generator 
of  the  wide  band,  backward-wave  oscillator  type  would  have  been  ideal, 
but  none  was  available.  Since  the  frequency  was  swept  over  a  50  Me  range 
the  heating  effects  were  greatly  reduced  and  ambient  temperature  drift 
was  no  longer  a  problem.  The  mechanical  sweep  generator  provided  a 
voltage  output  proportional  to  the  angular  displacement  of  the  frequency 
dial.  The  output  voltage,  which  was  then  linearly  proportional  to  the 
differential  frequency,  was  fed  to  the  horizontal  input  of  the  oscillo¬ 
scope.  The  vertical  input  of  the  oscilloscope  came  from  the  feflec tome ter , 
which  indicated  the  reflection  coefficient  of  the  dielectric  resonator 
system.  A  camera  was  utilized  to  take  photographs  of  the  oscilloscope 
trace  of  reflection  coefficient  vs  frequency. 

Given  that  Q  ^  is  constant,  the  type  of  coupling  (over  or  under¬ 
coupled)  can  be  determined  by  noting  the  dependence  of  Qq  upon  po  or 
rQ  through  (5*18)  and  (5.20). 

Another  means  of  determining  the  type  of  coupling  follows  from 
(5.18)  and  (5.20),  with  the  knowledge  that  the  ferrite  losses  lncrsase 
as  the  dc  magnetic  field  HQ  is  decreased.  Ths  type  of  coupling  is 
determined  by  observing  the  change  in  reflected  power  on  resonance  when 
increased  loss  is  introduced  into  the  system  by  reducing  Hq  .  Thus, 
for  an  overcoupled  system,  increasing  the  loss  causes  the  reflected  power 
to  decrease;  whereas  if  the  system  is  undercoupled,  the  reflected 
power  will  increase. 

Knowing  whether  the  signal  source  was  over  or  under  coupled  to  the 
resonator  system,  it  was  possible  to  determine  the  Q  of  the  resonator 
system  from  one  of  these  photographs.  First,  the  reflection  coefficient 
rQ  on  resonar.ee  was  determined.  Then  the  coupling  coefficient  $ 
was  found  from  (3*21)  and  (5. 20).  Knowing  0  ,  the  VSHR  at  the  half¬ 

power  points  was  determined  from  a  smith  chart,  or  by  the  analytical 
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FIG.  5.5— Arrangement  used  Id  the  measurement'  of  Q_.  of  the  dielectric  resonator  system  for  experiment  I. 


expression 


30 


P(  1  +  P)  +  1  +  (l  +  P)Vl  +  p^1 
P 


(5.22) 


The  reflection  coefficient  corresponding  to  (5*22)  nay  be  obtained, 
and  the  resulting  frequency  difference  df^yg  between  half -power  points 
nay  be  read  from  the  photograph's  plot  of  reflection  coefficient  vs 
frequency.  With  fQ  as  the  resonant 
defined  by 


(5.23) 


frequency,  the  loaded  ^  is 


Finally,  the  unloaded  may  be  obtained  from  (3. 17). 

4.  Half -Power  Method 

The  half -power  method  was  a  method  employed  to  measure  the  Q  for 

experiment  II.  This  method  was  used  because  it  had  greater  accuracy 

and  was  much  quicker  than  the  sweep  frequency  method.  The  experimental 

arrangement  is  shown  in  Fig.  5-6.  The  pulse  length  of  the  signal  was 

less  than  2  nsec,  and  the  repetition  rate  was  20  cps.  The  incident 

power  level  was  about  10  dB  or  so  lower  than  the  maximum  output  of  the 

signal  generator  (~lmW).  Hence,  the  average  power  introduced  into  the 

resonator  system  was  lower  than  the  average  power  required  for  the  sweep 

.5 

frequency  method  by  a  factor  of  at  least  10  .  Consequently,  microwave 

heating  effects  were  completely  negligible.  Alsq,  any  time  delay  between 
taking  the  data  points  was  completely  justified  as  long  as  the  ambient 
temperature  did  not  change  appreciably  during  that  time. 

First,  the  frequency  of  the  dielectric  resonator  system  was  adjusted 
for  resonance  and  the  value  of  the  attenuator  noted  for  an  oscilloscope 
deflection  of  0.1  volt.  Then  the  system  was  tuned  off  resonance  by  the 
tuning  plunger,  and  the  attenuator  reading  taken  for  the  same  oscilloscope 
deflection.  The  attenuator  was  next  set  at  the  half -power  points  and 
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FIG.  5 •  6 — Arrangement  used  in  the  measurement  of  0.  of  the  dielectric  resonator  system  for  experiment  II 


the  frequency  measured  on  each  side  of  resonance.  Thus  the  loaded 
could  he  determined  directly  from  (5.2J).  Knowing  the  reflection 
coefficient  on  resonance,  could  then  he  determined  from  relations 

(5.19)  through  (5.21). 

5.  Results 

The  results  obtained  for  the  unloaded  as  a  function  of  the 
applied  dc  magnetic  field  HQ  for  experiments  I  and  II  are  shown  in 
fig.  5* 7a.  Normalized  curves  showing  the  approximate  dependence  of 
Qq  upon  the  frequency  tuning  plunger  position  are  shown  in  the  ac¬ 
companying  Fig.  5.7b.  Actually,  these  curves  also  dependsllghtly  upon 
the  magnitude  of  Qq  j  i.e. ,  the  lower  the  ,  the  nearer  the  plunger 
must  he  to  the  dielectric  resonator  for  a  given  relative  increase  in  . 
Nevertheless,  the  curves  shown  in  Fig.  5* 7b  are  sufficiently  accurate 
over  the  range  of  frequency  plunger  positions  for  which  threshold  measure¬ 
ments  were  taken.  The  increase  in  as  the  plunger  is  screwed  inward 
is  attributed  to  a  reduction  in  the  power  being  radiated  from  the 
resonator.  The  increased  by  about  the  same  percentage  whether  the 
ferrite  was  present  or  not.  The  deviation  of  at  the  nearest  plunger 
position  from  that  at  the  off -resonant  plunger  position  was  not  more 
than  15^.  The  position  of  the  frequency  tuning  plunger  had  no  noticeable 
effect  upon  the  resonator  coupling  coefficient  or  the  external  . 

The  decrease  of  at  increasingly  low  dc  field  values  shown  in 
Fig.  5 -7a  is  due  to  the  increasing  proximity  of  the  pump  frequency  to 
ferromagnetic  resonance.  At  large  dc  fields  the  ferrite  loss  becomes 
negligible  and  approaches  the  unloaded  Q  of  tb*l  dielectric 
resonator  without  the  ferrite,  since  energy  storage  in  the  ferrite  is 
negligible  compared  to  that  of  the  resonator.  The  dashes  at  the  far 
right  of  Fig.  5*  7a  indicate  the  measured  values  of  Qq  without  the 
ferrite. 

An  improvement  in  the  resonator  system  from  experiment  I 

to  experiment  II  is  apparent  from  Figs.  5- 7»,b.  In  the  two  cases  the 
Qq'  s  were  measured  with  different  plunger  positions,  but  this  accounts 
for  only  about  a  10jt  change  in  .  The  remaining  change  is  attributed 
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to  a  reduction  In  the  Mount  of  Ducco  cement  at  the  contact  between  the 
resonator  and  the  supporting  dielectric  rod,  where  there  was  strong  rf 
electric  field,  and  to  a  thorough  cleaning  of  the  waveguide  system. 

So  appreciable  dependence  of  the  unloaded  Qq  upon  the  position 
of  the  movable  Short  was  noted.  However,  the  range  of  movement  of  the 
movable  short  to  obtain  proper  coupling  was  small  in  both  experiments. 

The  variation  of  with  temperature  in  the  ambient  temperature 
rsnge  was  net  appreciable.  This  is  as  expected,  on  the  basis  of  the 
small  temperature  dependence  of  tan  6  given  in  Appendix  A.  In  any 
case,  the  measurements  of  Qq  were  taken  at  nearly  the  same  ambient 
temperature  as  were  the  threshold  measurements. 

F.  SUMMARY 

To  achieve  the  large  rf  magnetic  field  amplitudes  required  in  order 
t  reach  the  phonon  threshold,  a  dielectric  resonator  was  employed 
rather  than  a  conventional  cavity.  With  the  strontium  titanate  dielectric 
resonator,  the  large  rf  field  can  be  obtained  with  fairly"  low  rf  power, 
breakdown  is  avoided,  and  a  large  ferrite  filling  factor  is  obtained. 

The  microwave  power  source  and,  in  particular,  the  microwave  circuitry 
associated  with  the  modulation  of  the  klystron  beam  current  are  described. 
Explanations  are  given  for  the  limitations  of  the  rf  power  pulse  from 
the  source.  Calibration  techniques  for  the  rf  power  level  and  dc  mag¬ 
netic  field  are  also  discussed.  Various  details  of  the  dielectric 
resonator  system  are  described— including  the  resonant  mode  configuration, 
coupling  of  rf  energy  to  the  resonator,  effect  of  temperature  changes 
upon  the  resonant  frequency,  and  microwave  heating  of  the  dielectric 
resonator  and  ferrite. 

An  expression  for  the  maximum  rf  magnetic  field  amplitude  (at  the 
position  of  the  ferrite)  was  derived  in  terms  of  the  unloaded  Q  and 
power  absorbed  by  the  resonant  system.  Reliable  measurements  of  the 
unloaded  Q  were  obtained  by  using  the  sweep  frequency  method  and 
the  half-power  method;  results  are  given.  The  unloaded  Q  was  found 
to  vary  significantly  with  the  dc  magnetic  field  owing  to  the  presence 
of  the  ferrite. 
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CHAPTER  VI 


THRESHOLD  MEASUREMENTS 

Saturation  effects  have  been  observed  in  two  experiment*  using  a 
disk  of  single  crystal  ZngY  (0.  in.  diam  x  0.0076  in. )  cut  with  th< 

easy  plane  lying  in  the  plane  of  the  disk.  Actual  measurements  of  ab¬ 
sorbed  power  by  the  system  were  taken  at  the  onset  of  saturation  as  a 
function  of  the  applied  dc  magnetic  field  HQ  .  The  onset  of  saturatic  . 
was  determined  by  observing  the  reflected  pump  pulse  and  noting  the  pow< 
level  at  which  the  growth  of  instabilities  first  became  evident.  There 
was  a  critical  power  level  at  which  the  instability  threshold  was  reach* 
and  a  transient  appeared  in  the  reflected  pump  pulse.  Knowing  the  ab- 
,  sorbed  power  at  the  threshold  and  the  unloaded  Q,  M.  *  function  of 

the  applied  dc  magnetic  field  H0  >  the  apparent  critical  magnetic 

fielu  ampliuirte  (h  )  at  the  ferrite  sample  was  determined  from  (5.16). 

ap 

A.  EXPERIMENT  I 

A  photograph  of  a  typical  power  pulse  incident  upon  the  resonator 
system  for  experiment  I  is  shown  in  Fig.  6.1a*  The  incident  power  pulse 
was  120  tisec  long  and  had  a  time  constant  of  about  309  The  decay  f 

the  incident  pulse  amplitude  was  due  to  sag  in  the  klystron  beam  current 
which  discharged  from  a  large  capacitor.  This  time  constant  could  be 
reduced  by  increasing  the  capacitor  voltage;  however,  this  caused  the 
distortion  of  the  pulse  to  increase.  Therefore,  a  pulse  was  sought  whic 
had  a  maximum  time  constant  for  a  minimum  amount  of  distortion.  The  dou  e 
image  In  the  photograph  of  Fig.  6.1a  was  due  to  ac  pickup. 

«  A  typical  reflected  power  pulse  below  the  instability  threshold  is 

indicated  in  Fig.  6.1b.  The  time  constant  of  this  pulse,  in  general,  wa 
the  same  as  that  of  the  incident  pulse.  This  does  not  appear  to  be  quit* 
true  in  comparing  Fig.  6.1b  with  Fig.  6.1a,  because  the  latter  was  taken 
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Time - 

(a)  Typical  incident  power 

(b)  Reflected  power  pulse 

pulse . 

below  threshold. 

(c)  Reflected  power  pulse  Just  (d)  Reflected  power  pulse 

above  threshold.  above  threshold. 


(e)  Reflected  power  pulse  showing  the 

saturation  effect  In  the  low  dc  field 
region. 


FIG.  6.1 — Pulse  waveforms  for  experiment  I.  Time  scale:  is 
20  usec/divisior.. 
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at  a  slightly  lover  capacitor  voltage#  (The  "ear"  at  the  end  of  the 
reflected  pulse  vas  due  to  frequency  modulation  of  the  microwave  signal. 

The  latter  arose  from  the  variable  load  which  the  klystron  presented  to 
the  cv  signal  generator  at  the  leading  and  trailing  edges  of  the  pulse. 

This  effect  vas  somewhat  reduced,  however,  in  experiment  I  because  the 
reflection  coefficient  vas  high  (|p|  *  0.3)  .  The  gradual  decay  of  the 

ear  at  the  trailing  edge  of  the  pulse  was  due  mostly  to  the  fairly  high) 

RC  product  of  the  crystal  detector  circuit,  but  was  also  due  in  part  to 
the  finite  time  constant  of  the  klystron  beam  current  switching  network. 
Although  the  time  constant  could  have  been  reduced  by  decreasing  the 
resistance  in  the  detector  circuit,  this  vas  not  done  because  it  would 
have  reduced  the  already  low  sensitivity  of  the  reflected  power  measurements 
[see  (4.49)  on  the  sensitivity  requirements]# 

Since  the  dielectric  resonator  system  vas  overcoupled,  the  effect  on 

the  reflected  pump  pulse  of  an  increase  in  ferrite  absorption  wee  to  decrease 

the  magnitude  of  the  pulse.  A  typical  example  of  e  reflected  pulse 

exhibiting  the  initial  appearance  of  increased  ferrite  absorption  is  shown 

in  Fig.  6.1c.  This  waveform  is  typical  of  those  from  which  the  apparent 

critical  fields  (h  )  were  determined.  As  discussed  above,  the  threshold 

up 

was  found  from  the  absorbed  power  level  occurring  at  the  point  where 
absorption  vas  first  noticed.  It  is  noted  from  the  figure  that  the  ferrite 
absorption  transient  builds  up  gradually  to  a  peak  and  then  dedays  back  to 
the  unperturbed  reflected: pulse  amplitude.  This  type  of  behavior  is  very 
similar  to  that  predicted  by  the  phonon  absorption  theory  of  Chapter  IV. 

Figure  6. Id  shows  a  typical  reflected  pump  pulse  at  an  absorbed  power 
level  above  the  threshold  level.  In  this  case  It  is  noted  that  the  ferrite 
absorption  is  greater  and  the  pulse  saturation  is  even  more  pronounced  then 
at  the  threshold  level.  This  behavior  also  exhibits  the  typical  absorption 
behavior  described  in  Chapter  IV. 

In  studying  the  precise  nature  of  the<  ferrite  absorption  it  will  be 
helpful  to  normalise  the  reflected  pulse  waveforms  of  Figs.  6.1c,d  to  an 
arbitrary  level  below  threshold  (Fig.  6.1b).  .  The  results  of  doing'  this1 
are  shown  in  Fig.  6.2.  These  waveforms  shoving  ferrite  absorption  are 
similar  to  those  predlctod  in  Chapter  IV  for  phonon  absorption  in  the 
ferrite.  In  particular,  the  normalised  theoretical  waveforms  of  Fig.  4.8 
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t(20  nsec/di vision) 


FIG.  6. 2- -Drawing  of  ref lac tad  power  waveforms  of  Figs.  6. lb ,c,d 
normalized  to  the  level  corresponding  to  6.1b. 


-  (r) 

are  comparable  to  the  normalized  experimental  waveforms  of  Tig.  6.2.'  7 

Figure.  4.8a  is  comparable  to  Fig.  6£b,the  reflected  power  waveform  before 

the  onset  of  saturation;  and  Fig.  4.8b  is  comparable  to  Fig.  6.2'c,  the 

reflected  power  waveform  when  the  phonon  absorption  first  appears.  In 

each  case,  the  absorption  during  the  pulse  builds  up  to  a  peak  and  then 

decays  gradually.  It  is  seen  that  the  waveforms  above  threshold  power 

levels  (Figs.  4.6c  and  6. 2d)  are  also  similar  in  nature*  Again,  the 

absorption  builds  up  to  a  peak  and  decays.  In  both  theory  and  experiment 

the  absorption  peak  Increases  with  increasing  power  levels. 

The  first  exception  to  the  theory  occurs  in  the  magnitude  of  the 

build-up  time  t  The  experimental  value  c>served  from  Fig.  6.2c  is 

ap 


^ap^exp 


32  Weec 


(6.1) 


which  is  much  shorter  than  the  predicted  value  of  t  *  104  u«ec  for 

experiment  I,  baaed  upon  a  phonon  of  1300.  If  the  linear  theory  worked 

out  In  Chapter  IV  is  correct,  then  the  shorter  observed  build-up  time  simply 

implies  that  the  phonon  la  actually  much  lower  then  1300.  In  fact,  it 

was  shown  in  Chapter  IV  for  (t^)^  ■  32  psac  that  a  ^  of  129  was 

obtained  from  the  linear  theory.  It  might  be  sufficient  to  accept  this 

lower  value  of  as  the  actual  value  except  for  whet  happens  to  the 

reflected  pulse  waveform  at  higher  pump  power  levels.  From  Figs.  6.2c,d 

the  build-up  time  is  seen  to  decrease  with  the  increasing  power  Irnvel* 

This  is  contrary  to  the  behavior  predicted  by  (4.6t)>  where  t _ wee 

ap 

expected  to  increase  with  increasing  power  levels,  end  eight  be  explained 
on  the  basis  of  a  variation  of  ^  with  the  4  tude  of  the  unstable 
phonon  as  was  suggested  in  Chapter  IV.  In  any  case,  it  is  clear  that  the 

"tp€rl“nt  d0“  “*  »«"“•  •  r*lubl*  ^ 87  %  • 


'  Although  the  curves  of  Fig.  4.8  were  calculated  for  the  values  of 
experiment  II,  it  is  clear  from  the  equations  of  Chapter  XV  that  the  shape 
of  the  curves  describing  the  phonon  absorption  would  be  the  sees  for 
experiment  1.  Since  V  is  much  larger  for  experiment  I  than  for  experiment  II 
however,  the  effect  of  phonon  absorption  would  be  proportionally  Idee  for 
experiment  I. 
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The  above  discussion  relates  the  observed  saturation  threshold 
behavior  at  the  higher  dc  magnetic  fields  to  the  phonon  saturation  theory 
of  Chapter  IV.  At  the  lover  dc  magnetic  fields,  quive  different  threshold 
observations  vere  made.  A  typical  reflected  pump  pulse  in  the  low  dc  field 
region  is  shovn  in  Fig.  6.1e.  Here  the  threshold  is  defined  by  the  appearance 
of  the  small  ripples  shovn  in  the  figure.  This  type  of  threshold  is 
attributed  to  the  presence  of  spin  wave  Instabilities,  which  apparently  • 
dominate  the  phonon  instabilities. 

Measurements  of  absorbed  power  vere  caken  as  a  function  of  the  external 
magnetic  field  Hq  at  the  threshold  levels  corresponding  to  Figs.  6.1c 
and  6.1e.  Apparent  critical  fields  vere  obtained  from  these  measurements 
and  are  plotted  in  Fig.  6.3*  The  two  typical  experimental  runs  shovn  vere 
taken  with  increasing  HQ  .  These  runs  vere  also  typical  of  runs  taken 
with  decreasing  HQ  .  In  order  to  verify  the  continuous  nature  of  the 
threshold  curve,  the  experimental  points  for  the  second  run  vere  taken  at 
the  smallest  Increments  of  possible  with  the  magnetic  field  current 
supply  used  in  experiment  I.  This  run  also  gives  a  good  indication  of 
the  experimental  error  to  be  expected  in  the  reproducibility  of  such 
measurements. 

The  ambient  temperature  was  slightly  different  for  the  two  runs 
indicated.  This  implies  a  possible  greater  disparity  in  the  temperature 
of  the  ferrite  Itself  for  the  two  runs.  In  accordance  with  results  from 
the  previous  chapter,  the  temperature  of  the  ferrite  depends  chiefly  upon 
the  length  of  time  it  is  exposed  to  a  given  microwave  field.  It  is  evident 
that  the  apparent  effect  of  slight  temperature  deviations  on  the  thresholds 
from  one  run  to  another  was  not  appreciable  and  is  within  experimental 
error.  However,  for  both  runs  it  Is  expected  on  the  basis  of  the  analysis 
of  Chapter  V  that  at  the  higher  power  levels  there  will  be  some  shift 
toward  lover  HQ  values  due  to  the  dependence  of  the  dmsagnetizing  field 
upon  temperature .  In  fact,  using  the  10°C  temperature  change  estimated 
for  ■  10  V  ,  a  downward  shift  of  8  Oe  would  be  expected  at 
HQ  -  365  0e  . 

The  relationship  between  the  apparent  theoretical  threshold  h 

ap 

for  ths  actual  pulse  and  the  threshold  h  _  for  an  ideal  pulse  was 

cr 


•  Bxperine ntal  points  for  first  typical  run  with  T  *  26°C 


*  Kxpsrlaental  points  for  sscflnd  typical  run  vith  /  T  -  2$°C 
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FIG.  6.3-* -Experimental  and  theoretical  values  Of  h 
for  thdi  ZaiY  diskin'  transverse 
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determined  In  Chapter  IV.  This  ratio  h  / h  is  only  approximately 

apr  cr 

constant  over  the  range  of  field  for  vhich  the  experimental  thresholds 

of  fig.  6.3  were  taken.  The  only  significant  change  in  this  ratio  takes 

place  through  variation  in  the  phonon  absorption  peak  (Pp)ap  Recall 

the  relation  (4.6l)  vhich  relates  P  to  the  reflection  coefficient  r 

P 

and  the  change  in  reflected  power  AP  Using  this  relation,  it  was 

found  from  the  experimental  results  that  (P  )  changes  by  a  factor  of 

P  aP 

50:1  ,  increasing  nonotonically  as  HQ  increases.  The  effect  of  this 

change  in  the  observed  phonon  absorption  peak  upon  the  ratio  h  /h 

ap'  cr 

about  a  mean  value  is  less  than  0.2^  for  the  parameter  values  associated 

with  (4.66),  for  the  effective  «  127  and  the  observed  build-up  time. 

Thus,  the  ratio  h ,  /h  is,  indeed,  a  constant  for  all  practiced,  purposes, 

cr 

and  the  above  assumption  is  thereby  Justified. 

The  dominant  theoretical  phonon  threshold  for  an  ideal  pulse,  given 
by  (3*52),  may  be  written 


h  K  x 


.3/2 


per 


Hq  1  +  4.67  x' 


(6.2).  , 


where 


2C  (4xM  ) 
sy 


b: 


and  it  is  recalled  that 


x'  *  x  ♦  d  , 


where 


H  .A 

4*m  4*m 
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Dk 
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The  apparent  threshold  will  then  be  given  by 


<Vap 


rrrari' 


(6.3) 


vhere  the  parameters  are  defined  as  above.  Since  the  ratio  ba^/hc 
Is  very  nearly  Independent  of  Hq  ,  as  explained  above,  Sq.  (6.3)  eas 
fitted  to  Fig.  6.3  by  adjusting  the  paraoeters  ( “* d  •  11 

was  found  that  the  best  fit  was  obtained  vlth 


and 


9.68 


(6.W) 


d  •  0.022 


(6.4b) 


Expression  (6.3)  is  plotted  vlth  these  pareaeter  values  in  Fig.  6.3,  vhere 
it  is  seen  that  experimental  points  fit  the  theoretical  curve  Quite 
closely. 

From  the  value  of  d  given  by  (6.. 4b),  it  is  found  that 

*  56  °*  • 

Knowing  the  lattice  constant  a  ,  and  the  value  of  k  for  the  unstable 
elastic  wave  mode  (determined  from  (3.28b)  and  (3»57)J,  the 
field  in  the  easy  plane  H^y  becomes 

H«x||  *  *2-3  x  lo7°«  ■  (6-5) 
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From  the  above  value  of  (hmp/hcr)  Kp  given  by  (6.4a)  and  the  theoretical 

value  of  b  /h  -  1.003  for  experiment  I  given  by  (4.66),  the  value  of 
ajr  cr 

K  obtained  is  9*65.  This  value  is  somewhat  larger  than  the  value  of 
P 

K  »  7.74  estimated  in  Chapter  III  and  given  in  (3.60).  The  latter  was 

p  6  3 

based  upon  the  estimates  of  *  1500  and  B^  ■  5  X  10  erg/ era'  .  Taking 

instead  the  effective  value  of  ■  127  given  by  (4.66)  and  the  same 

values  of  C. ,  and  4*M  employed  originally,  a  new  estimate  may  be 

a"  '  s . ;  v  ■  ■  ■  r  ’  cr  re 

obtained  for  .  From  the  expression  for  given  in  (6.2),  the 

following  value  for  is  ob "kilned: 


*  1.54  x  10*  erg/cm^  . 


(6.6) 


Thus  it  appears  that  the  magnetoelastic  constants  for  Zn^£  are  larger  0 
than  they  are  for  YIG  (B  *  '5  x  10^  erg/cm? ). 

With  the  parameters  given  in  (6.4)  the  theoretical  threshold 
expression  (6.3)  fits  the  experimental  points  of  Fig.  6.3  veil.  There 
are,  however,  two  minor  considerations  which  should  be  pointed  out.  The 
first  one,  which  has  already  been  described,  is  the  effect  of  microwave 
beating  on  the  saturated  magnetization  and  hence  upon  the  applied  dc 
magnetic  field  required  to  sustain  a  given  internal  field.  This  effect 
explains  the  deviation  of  the  experimental  points  above  the  theoretical 
curb?  at  the  higher  threshold  levels  shown  in  Fig.  6.3*  Since  the  effect 
of  microwave  heating  at  the  higher  thresholds  decreases  the  demagnetizing 
field,  the  theoretical  curve  shown  in  Fig.  6.3  more  closely  approximates 
a  curve  of  experimental  points  taken  at  the  same  ferrite  temperature ,  rather 
than  one  of  points  taken  at  the  same  ambient  temperature.  The  second 
consideration  has  to  do  with  deviation  of  the  theoretical  threshold  from 

that  shown  in  Fig.  6.3  for  the  lover  H  values.  In  deriving  the  trans- 

u  2  2 

verse  phonon  thresholds  the  assumption  was  made  1  thht  o>  «  o>  .As  was 

p  m 

mentioned  In  Chapter  III,  this  condition  is  fulfilled  very  veil  with  dc 

2  <  2 

fields  above  where  oqT  a  0.1  o>m  .  However,  in  the  transition  range, 

**02  <  Ho  <  ho3  1  theoretical  phonon  threshold  is  somewhat  lover  than 
indicated  in  Fig.  6.3.  If  both  of  the  above  effects  are  taken  into 
account  in  the  interpretation  of  the  experimental  points  and  theoretical 
curve  of  Fig.  6. 3>  it  is  expected  that  slightly  different  values  for  the 
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exchange  constant  d  and  the  multiplicative  constant  would  he 

obtained.  However,  the  corrections  which  might  be  asde  would  be  within 

experimental  error  of  the  points  of  Fig.  6.3. 

At  dc  fields  ~  HQ2  in  Fig.  6.3  the  reflected  power  waveforms 

exhibited  a  different  kind  of  saturation  effect.  This  Is  attributed  to 

the  growth  of  spin  wave  instabilities.  The  spin  wave  thresholds  for  an 

ideal  pulse  were  plotted  as  a  function  of  internal  dc  field  in  Figure  2.3. 

These  are  related  to  the  critical  Held  h  by  the  familiar  definition 

cr 


(A  ) 

'  p  er 


cr 

80 


(6.7) 


For  the  actual  pump  pulse  the  apparent  threshold  is 

(V.P  ‘  f~  <V«  •  <6-6> 

ncr 

It  will  now  be  assumed  that  h^^/h^  is  a  constant ,  Independent  of  Hq  , 
and  that  this  constant  is  approximately  unity.  This  assumption  is 
reasonable  in  view  of  the  results  for  the  build-up  of  phonon  instabilities 
(Chapter  IV).  Then  the  apparent  spin  ware  thresholds  are  given  approximately 
by  Fig.  2.3,  where  the  spin  wave  Q  has  been  assumed  to  be  m  193  .  A 
comparison  of  the  theoretical  curves  of  Fig.  2.3  with  the  experimental 
points  of  Fig.  6.3  shows  that  the  experimental  thresholds  art  larger.  It 
is  tentatively  deduced  from  this  that  tbs  spin  wave  Q  Is  lower  then 
193.  By  matching  the  theoretical  on-resonance  second  orddf  threshold 
curve  (2.30)  to  the  level  established  by  the  three  data  points  near 
Hq  -  270  Oe,  an  estimate  Q  *  23.4  ,  is  obtained  for  tbe  effective  spin 
wave  Q  .  Using  this  value  of  Q  ,  theoretical  values  of  the  first  end 
second  order  spin  wave  thresholds  have  been  plotted  in  Fig.  6.3.  It  is 
seen  that  the  lowest  threshold  changes  from  spin  wave  to  phonon  at 
Hq  s  Hq2  The  second  order  off -resonance  spin  wage  threshold  is  seen 
to  Increase  very  steeply  above  Hq  -  HQ2  ,  end  is  the  dominant  theoretical 
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threshold  only  over  a  negligibly  snail  range  of  HQ  *. 

At  further  reduced  magnetic  fields,  ,  the  experimental 

thresholds  are  seen  to  increase  rapidly- -in  a  region  where  the  very  low 
first  order  threshold  is  expected  to  be  dominant.  The  explanation  for 
this  probably  lies  in  the  fact  that  in  this  region  the  internal  dc  magnetic 
field  is  insufficiently  large  to  saturate  the  magnetization  of  the  ferrite. 
Thus  the  magnetization  no  longer  is  uniform  throughout  the  .sample  volume , 
and  additional  magnetic  relaxation  processes  maty  be  operating. 

B.  EXPERIMENT  XI 

This  experiment,  similar  to  experiment  I ,  was  undertaken  to  observe 
in  greater  detail  the  more  exact  nature  of  the  saturation  effects  on  the 
reflected  pump  pulse.  The  magnet  used  for  this  second  experiment  had  a  more 
uniform,  better  regulated,  and  larger  range  of  magnetic  field  than  did  the 
one  employed  in  experiment  I.  The  magnetic  field  was  sufficiently  uniform 
for  it  to  be  measured  directly  with  the  NMR  probe.  The  larger  range  of 
the  magnetic  field  HQ  allowed  measurement  of  thresholds  corresponding 
to  the  maximum  power  output  of  the  signal  generator.  This  made  it  possible 
to  p...  through  th.  -ulmum  fl.ld  at  which  th.  fifth  order  .pin 

wave  instability  would  be  permitted  on  resonance. 

Also,  in  this  experiment  the  level  at  which  the  reflected  power  could 
be  detected  was  20  dB  lower  than  in  experiment  X,  permitting  a  reduction 
in  the  reflection  coefficient  r  of  the  resonant  system;  and  the  unloaded 
Qq  of  the  dielectric  resonator  system  was  higher  than  it  was  for  experi¬ 
ment  X,  with  the  result  that  the  incident  power  level  P,  was  reduced. 

The  effects  of  lower  r  and  P1  on  the  relative  change  in  reflected 
lower  AP due  to  a  given  change  in  ferrite  absorption  are  seen  by  an 
examination  of  (4.49),  This  shows  that  APr  is  increased;  the  saturation 
offsets  can  therefore  be  obser^/ed  In  greater  detail.  In  this  experiment 
much  better  photographs  of  saturation  phenomena  on  the  reflected  pulee 
were,  in  fact,  obtained. 

A  photograph  of  a  typical  klystron  beam  current  pulse,  as  seen  at 
the  output  of  a  current  pulse  trana&rmer,  is  shown  in  the  top  half  of 


Fig.  6.4a.  Some  indication  of  the  sag  in  the  current  discharging  from  the 
high  voltage  capacitor  in  the  klystron  modulator  is  evident  in  the  first 
30  usee.  The  increased  sag  beyond  30  usee  is  due  to  sag  inherent  in  the 
current  transformer  as  veil  as  in  the  beam  current.  The  amplified  rf 
power  pulse  at  the  output  of  the  klystron  detected  at  the  directional 
coupler  is  shown  in  the  bottom  half  of  Fig.  6.4a.  The  pulse  was  90  ueee< 
in  duration  and  had  a  time  constant  of  about  124  usee.  This  pulse  has  a 
shorter  length  and  a  smaller  time  constant  than  the  incident  pulse  of 
experiment  I  because  there  was  an  equipment  failure  and  unavoidable 
alterations  had  to  be  made  in  the  klystron  modulator.  For  each  experiment, 
a  pulse  of  maximum  time  constant  was,  of  course ,  sought. 

A  typical  reflected  power  pulse  below  the  instability  threshold  is 
indicated  in  Fig.  6.4b.  The  ears  on  the  reflected  pulse ,  due  to  the 
frequency  modulation  described  previously,  are  more  pronounced  than  they 
were  for  experiment  I  because  the  reflection  coefficient  i As  much  smaller 
( |r|w  0.08).  The  time  constant  of  the  ears  is  smaller  than  for  experiment 
I  because  the  resistance  In  the  crystal  detector  circuit  was  decreased. 

It  was  possible  to  decrease  the  resistance  for  experiment  II  because  the 
sensitivity  (or  the  relative  change  In  reflected  power)  was  greater,  for 
reasons  mentioned  above. 

In  this  experiment  the  dielectric  resonator  system  was  again  over¬ 
coupled.  Thus  the  effect  of  an  Increase  in  ferrite  absorption  was  to 
decrease  the  magnitude  of  the  reflected  pump  pulse.  A  typical  example  of 
a  reflected  pulse  exhibiting  the  Initial  appearance  of  ferrite  absorption 
is  shown  in  Fig.  6.4c.  This  vavmform  is  typical  of  those  from  which  the 
apparent  critical  fields  (hap)  were  determined.  It  is  apparent  that  the 
ferrite  absorption  gradually  builds  up  to  a  peak  end  then  decays  as  In 
experiment  I.  However,  in  experiment  II  the  initial  appearance  of  the 
absorption  peak  is  more  clearly  seen,  and  the  similarity  to  the  waveform 

r  ■  i 

predicted  by  the  theory  of  Chapter  IV  1»  even  more  narked. 

Figure  6.4d  shows  a  typical  reflected  pump  pulse  at  an  absorbed 
pc  /©r  level  above  threshold.  The  increase  in  ferrite  absorption  is  tv>  . 
apparent.  At  a  still  higher  level  above  threshold,  the  reflected  puap 
pulse  appears  as  shown  in  Fig.  6.4e.  Hire  the  absorption  peak  appears 
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Time  (20  ^sec/division) 


(a)  Klystron  beam  current  and  rf 
output  power  waveforms. 


(c)  Reflected  power  pulse  at 
threshold. 


(b)  Typical  reflected  power 
pulse  below  threshold. 


—  Time  (10  ^sec/division) 

(u)  Reflected  power  pulse  above 
threshold . 


(e)  Reflected  power  pulse  much 
above  threshold. 


FIG.  6.4 — Pulse  waveforms  for  experiment  II. 


slightly  truncated- - indicating  the  presence,  possibly,  of  pulse  saturation 
due  to  the  change  in  r  ,  as  described  in  Chapter  IV.  When  the  changes 
in  reflected  power  become  as  large  as  those  shown  In  Fig.  6.4e,  this  ' 
effect  should  be  appreciable.  If  the  ferrite  absorption  Increases  much 
more  than  that  indicated  in  Fig.  6.4e,  the  reflection  coefficient  will 
"go  through  zero"  at  the  absorption  peak  and  the  dielectric  resonator 
system  will  become  undercoupled. 

In  studying  the  precise  nature  of  the  ferrite  absorption  it  is  again 
helpful  to  normalize  the  reflected  pulse  waveforms  of  Figs.  6.4b-e  to  an 
arbitrary  level  below  saturation.  The  results  of  the  normalization  are 
shown  in  Fig.  6.5.  These  waveforms  showing  ferrite  absorption  are  similar 
to  those  shown  in  Fig.  4.8  for  the  theoretical  phonon  absorption  based 
upon  the  parameters  of  experiment  II.  In  particular,  it  is  noted  that 
the  predicted  waveforms  of  Figs.  4.8a-c  correspond  with  the  obstrved 
waveforms  of  Figs.  6.5b»d.  In  both  theory  and  experiment  the  absorption 
peak  increases  with  increasing  power  levels  above  threshold.  The  shape 
of  the  absorption  is  approximately  the  same  in  each  case.  There  appear, 
however,  to  be  the  same  deviations  from  the  theory  as  were  observed  in 
experiment  1. 

The  first  exception  to  the  theory  occurs  In  the  magnitude  of  the 

build-up  time  t  .  .  The  experimental  value  observed  from  Pig.  6.3c  Is 

ap 


(t  ) 

'  ap  exp 


26  usee 


which  is  shorter  than  the  predicted  value  of  t  •  60  usee  for  experiment 
II,  based  upon  a  phonon  of  1500  .  It  is  noteworthy,  though,  that 

the  experimental  value  of  t  is  closer  to  the  predicted  value  for  the 
parameters  of  experiment  II  than  for  those  of  experiment  I.  This,  as  will 
be  discussed  below,  is  <  attributed  to  the  lover  ferrite  absorption  level 
of  experiment  II.  From  the  linear  theory  of  Chapter  IV  it  may  be  concluded 
that  the  shorter  observed  build-up  time  implies  that  the  phonon  &  Is 
actually  much  lower  than  15 00.  In  fact  it  was  shown  In  Chapter  IV  that 
a  of  294  was  obtained  for  the  linear  theory  for  (t  ■  28  usee. 

As  explained  previously,  this  value  of  ^  may  probably  ba  regarded  as  an 
effective  value  which  can  be  conaidered  as  decreasing  with  increasing 
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FIG.  6. 5 ‘-Drawing  of  raflactad  powar  vava forms  of  Figs.  6.1b,c,d 
normal! 2«d  to  tha  laval  corrasponding  to  6.1b. 
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phonon  amplitude.  Thus  the  phonon  damping  ie  considered  nonlinear.  By 

this  Interpretation  of  the  results ,  it  is  possible  to  explain  the  decrease 

in  the  observed  build-up  time  t  with  increasing  absorption  levels, 

ap 

which  can  be  clearly  seen  in  Fig.  6.5. 

The  above  discussion  correlates  the  observed  instability  thresholds 
at  the  higher  dc  magnetic  fields  with  the  theory  of  Chapter  IV.  At  the 
lower  dc  magnetic  fields,  quite  different  threshold  observations  were  made. 

In  this  experiment  several  different  kinds  of  saturation  effects  were 
observed  for  HQ  <  HQ2  .  A  typical  example  of  the  reflected  power  wave¬ 
form  taken  in  this  region  is  shown  in  Fig.  6.6a.  The  dc  field  of 
Ho  "  ^ 00  for  these  waveforms  Is  within  14  Oe  of  HQ2.  lech  successively 

higher  waveform  in  the  photograph  occurs  for  a  higher  incident  power  level > 
and  hence  a  higher  field  strength  at  the  ferrite.  The  lowest  trace  in 
the  photograph  shows  a  slight  evidence  of  ripples  et  the  lowest  power 
level.  At  e  slightly  higher  power  level  the  ripples  show  their  somewhat 
random  nature  and  their  amplitudes  increase.  The  third  waveform  shows 
superimposed  upon  the  random  ripplee  the  eppearance  of  the  absorption 
peak  associated  with  the  phonon  theory.  At  the  highest  power  #sv«fom  Indicated 
in  the  photograph,  the  ripplee  are  dearly  dominated  by  the  phonon-like 
behavior.  At  da  fields  higher  then  300  Oe  the  amplitude  of  the  email 
random  ripples  became  less  and  less,  until  at  about  HQ  -  360  Oe  they 
vanished  completely.  At  fields  less  then  300  Oe  the  amplitude  of  the 
ripples  became  larger,  but  they  were  still  dominated  by  the  phonon 
absorption  dip  described  shove  until  the  point  Hq  •  HQ2  was  reached. 

For  Hq  <  H  random  ripples,  similar  to  those  described  above,  appeared 
quite  suddenly  at  a  critical  power  level.  This  is  believed  to  define  the 
second  order  spin  wave  threshold.  For  HQ  >  UQ2  the  appeareaoe  of  the 
absorption  dip  is  believed  to  define  the  phonon  threshold,  even  when 
random  ripples  are  superimposed.  The  reason  for  the  superposition  of  the 
two  patterns,  suggesting  concurrent  spin  wave  end  phonon  instabilities , 
is  not  known. 

For  Hq  <  Hq,  another  kind  of  saturation  effect  we*  observed.  This 
is  «hown  in  Fig.  6.6b,  for  which  Hq  -  263  0#  .  There  wee  a  threshold, 
marksd  by  the  suddsn  appearance  of  relaxation  oscillations  in  the  form 
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Time 


(a)  Saturation  effects  for  HQ,  slightly 
greater  than  . 


(b)  Saturation  effects  for  Hq  leuu  than  . 


(c)  Saturation  of fee to  for  HQ  slightly 
leuu  than  . 


FIG.  6. 6- -Reflected  power  pulse  vavoforae  shoving  various  saturation 

effects  m observed  in  experiment  II.  Time  scale  lu  10  |iuec/di vision. 

-  1^4  - 


of  an  attenuated  sinusoid.  One  reason  for  the  attenuation  was,  of  course, 
that  the  incident  power  decreased  during  the  pulse  and  eventually  reached 
a  level  below  the  threshold.  This  waveform  was  reproducible  and  bad  a 
relaxation  frequency  of  about  250  kc.^  This  effect  was  observed  only 
at  dc  field  values  between  Hg  and  HQ1  ,  and  is  attributed  to  the  first 
order  resonant  spin  wave  instability.  For  the  slightly  higher  dc  field 
H0  «  265  Oe  a  similar  behavior  was  observed,  as  is  shown  in  the  top 
waveform  of  Fig.  6.6c.  In  this  case,  however,  as  the  power  level  was 
Increased  the  waveform  shown  in  the  bottom  of  Fig.  6.6c  was  observed. 

Here  the  attenuated  sinusoid  broke  up  into  somewhat  random  oscillations 
as  the  power  was  Increased  further.  This  effect  is  not  understood. 

Msasuremsnts  of  absorbed  power  were  taken  as  a  function  of  the 
applied  magnetic  field  HQ  at  the  threshold  level  for  each  type  of  ferrite 
absorption.  Apparent  critical  fields  were  obtained  from  these  measurements 
and  are  plotted  in  Fig.  6.7.  Proceeding  as  before  with  expressions  (6.2) 
and  (6.3),  a  "best  fit"  of  the  theoretical  phonon  threshold  curve  to  the 
experimental  points  of  Fig.  6.7  is  obtained  with 

h 

K.  -  10.35  (6.9a) 

hcr 

d  -  0.022  .  (6.9b) 

Expression  (6.3)  is  plotted  with  these  parameter  values  in  Jl^>6.7>  where 

it  is  seen  that  the  experimental  points  fit  the  theoretical  curve  closely. 

The  value  of  d  in  (6.9b)  is  the  same  as  was  obtained  from  the  threshold 

ourve  for  experiment  I.  Thus  the  calculated  exchange  field  is  the  same 

as  before,  equation  (6.5)*  r' 

From  the  value  of  (h* ap/bcr)^  6iv*n  by  (6.9a)  and  the  theoretical 

value  of  h  Jh  *  0.961  for  experiment  II,  given  by  (4.79),  the  value 
a y  cr 

“  “"^Such  relaxation  oscillations  have  frequently  been  observed  in 
parallel  pump  experiments.  They  are  not  completely  understood. 
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FIG.  6.7- -Experimental  and  theoretical  values  of  h  vs  HQ  for  the  Zn2Y 
disk  in  transverse  pumping  experiment,  II.  p 
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of  K  obtained  ic 
P 

K  *  10.77  r 
P 

which  is  comparable  to  the  value  of  obtained  for  the  first  experiment. 
Again,  taking  the  above  value  of  ,  the  effective  value  of  294 

given  by  (4.67)  end  the  same  values  of  C.*  and  4*11  used  originally  , 
the  following  estimate  for.  is  obtained: 

B3  -  0.959  x  10T  erg/cm3  .  (6.10) 

This  result  is  comparable  to  the  result  (6.6)  for  experiment  I. 

As  in  experiment  I,  the  thresholds  attributed  to  the  second  order 
spin  wave  Instability  may  be  fitted  to  the  theory  by  assuming  a  spin  wave 
Q  given  by  Q  *  18.0  ,  and  the  theoretical  first  and  second  order 
thresholds  sure  plotted  in  Fig.  6.7,  assuming  this  value  of  Q  .  The 
results  are  comparable  with  those  of  experiment  I. 

C.  SUMMARY  AND  CONCLUSIONS 

Two  experiments  were  performed  using  a  disk  of  single  crystal  Zi^Y 
cut  with  the  easy  plane  lying  in  the  plane  of  the  disk.  In  each  experiment 
the  apparent  critical  magnetic  field  was  found  to  be  a  function  of  the 
applied  dc  magnetic  field.  The  instability  threshold  (or  critical  field) 
was  determined  by  the  absorption  peak  in  the  reflected  pump  pulse* 
Photographs  of  the  incident  and  reflected  power  pulse  waveforms  are  shown 
and  irregularities  in  the  waveforms  due  to  the  nature  of  the  power  source 
are  explained.  Irregularities  on  the  reflected  pulee  waveforms  due  to  the 
growth  of  ferrite  instabilities  are  discussed  in  detail. 

In  the  dc  magnetic  fie]d  region  where  phonon  instabilities  were 
expected  to  be  dominant,  the  shape  of  the  ferrite  absorption  on  the 
reflected  pump  pulse  was  nearly  the  sans  ae  that  predicted  by  the  phonon 
absorption  theory  of  Chapter  IV.  The  observation  of  instability  thresholds 
was  found  to  be  approximately  as  predicted  by  the  theoretical  spin  wave 
and  elastic  wave  threshold  calculations  of  Chapters  II  and  III.  At  applied 
dc  magnetic  fields  above  the  point  where  the  seond  order  spin  wavs  threshold 
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becomss  off -resonant,  the  experimental  threshold  fields  were  found  to 
exhibit  the  sane  functional  dependence  upon  dc  field  as  the  theoretical 
phonon  threshold  curve.  Prom  this  curve  an  estimate  was  obtained  for  the 
exchange  field  H^y  of  Za2Y  in  the  easy  plane.  This  functional 
dependence  strongly  indicates  that  the  magnetoelastic  effect  is  the 
dominant  mechanism  and  that  the  observed  thresholds  were  due  to  growth  of 
unstable  elastic  waves  propagating  and  polarized  in  the  easy  plane.  At 
applied  dc  magnetic  fields  below  the  point  where  the  second  order  spin  wave 
threshold  becomes  off -resonant,  the  experimental  threshold  fields  were 
found  to  exhibit  qualitatively  the  same  dependence  upon  dc  field  as  the 
theoretical  first  and  second  order  spin  wave  threshold  curves,  with  some 
exceptions. 

Explanations  for  these  exceptions  are  suggested  but  more  work  must 
be  done  in  order  to  understand  the  precise  nature  of  these  instabilities, 
and  it  should  probably  be  done  with  ideal  pulses.  In  the  phonon  threshold 
case,  however,  the  use  of  actual  pulses  permitted  the  determination  of 
both  the  magnetoelastic  constant  and  ths  phonon  ;  but  there  Is 
some  question  as  to  the  reliability  bf  these  values. 
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APPENDIX  A 


PHYSICAL  CONSTANTS  FOR  Zr^Y  ,  YIO  ,  ’AND  SrTi  0^ 

This  appendix  presents  a  table  of  those  constants  which  were  useful 
during  this  study,  with  special  Interest  In  the  single  dryutal  properties 
of  each  Material.  The  single  crystal  values  of  the  saturated  Magnetisation 
were  assumed  to  be  the  same  as  the  polycrystalline  values  given,  since  the 
crystal  anisotropy  was  overcome  at  11  kOe. 

The  values  of  saturated  Magnetisation  and  anisotropy  at  2$>°C  were 
used  extensively  throughout  this  study  and  are  given  below  for  reference: 

General  Constants 


4*M#  -  2624  Os 

Ha  -  8460  Oe 

P  •  Ha/4*H#  •  3.22 

Uip/2*  -  1.32  kMe 

q  m  fD^/Zu^  *  0.0898 
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Property* 

Zn£ 

YIO 

Strontium 

Titanete 

Formula 

^n2Ba2^e12°22 

Y3Pe5°12 

SrTi03 

Molecular  Weight  • 

1428 

738 

183 

Density  (gm/cm3) 

5.46 

5*17 

5.1337 

Curie  temperature  (°K) 

403 

560 

Crystal  type 

Hexagonal 

Cubic 

Cubic 

Lattice  constant^*)  of  unit  cell  (X) 

a-5.88 

c-14.5 

a-12.37 

a-3.9038 

initial  permeability  at  29°C 

Polycrystal 

Single  crystal 

■ 

28 

700  eot. 

Saturation  magnetisation  4xM  (Oe)  at  20°C 
Polycr^stal  at  11  KOe 

Powder  at  11  KOe 

Temp,  derating  (Oe/°C)  near  r.t. 

2700.,, 

262831 

-8.4 

1700 

-3*7 

Anisotropy  field  H.(Oe)  at  20°C 

Temp*  derating  (Oe/°c)  near  r.t. 

9000  32,33 
-60  est.>:,JJ 

^35,9 

-0.7  est 

3,9 

• 

Dielectric  constant  at  27°C  and  1.31  kMo 
Temp,  derating  (l/oc^1)  near  r.t. 

Lous  tangent  at  1.3  kMc 

Temp,  derating  (l/°Ca>'L)  near  r.t. 

10  est. 

2771'  28 

-1.093  ° 

2.  2 . 5j'x10-7i(U 
<  3.63x10  ' 

"1  o  *1 

Specific  heat  c(cal  gnf  C  ) 

Slant io  constants  (erg/cm3) 

e^lOrfrlO** 

c^-7.64xlOU 


0.17  est. 


34 


0.2  out. 


41 


c.. -26.9x10 


11 


36 


* 

The  source  for  these  date  lu  reference  9  unleu w  otherwise  Indicated. 
*8ee  Chapter  V. 
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APPENDIX  B 


GLOSSARY  OF  HUKCIPAL  SYMBOLS 


A 

Ai 

V 


^Ap^ap 

(n) 
cr 


(V 


^Vij 


cr 


Generalized  vector  which  may  represent  a  ,  £  ,  etc. 
Component  of  generalized  vector  along  coordinate  direction  i 
Component  of  generalized  vector  in  rotated  coordinate  system 
Pumping  angle  ■  b/HQ 
Apparent  threshold  angle  »  ^ap/^o 

1L 

n  order  spin  wave  threshold 

Phonon  threshold  as  a  function  of  rotation  angle  tQ  ;  i,J 
refer  to  the  direction  of  polarization  and  propagation, 
respectively 

Phonon  threshold  (minimized  with  respect  to  t^) 


0 

b1#b2,b3,b^ 
c&ftnT»'otc ' 


i 


Time -independent  spin  wave  amplitudes 
Magnetoelastic  constants  for  uniaxial  crystal 
Elastic  constants  for  uniaxial  crystal 

Normalized  exchange  field  -  Dk2/4*Mg  ;  ferrite  disk  diameter 

Free  energy  of  crystal 

Anisotropy  energy 

Magnetpcrystalllne  energy 

Elastic  energy 

Magnetostrlctlve  energy 

Phonon  energy --the  portion  of  X  in  the  unstable  mode 

V 

l 

Time  average  stored  phonon  energy  in  sample 
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G 

fi 

HA 

Sa 

»d 

Sd 

^dip 

Sex 

So 


H 


On 


h 

h 


ap 


cr 


m 


mp 


h 

~P 


y 


K  yK 

mp  pm 


Frequency  modulation  parameter 
Total  magnetic  field  «  H(r,t) 

Anisotropy  field  along  "hard"  or  c-axis 
Anisotropy  field 

DC  demagnetizing  field  in  the  plane  of  the  ferrite  disk 
Demagnetizing  field  for  uniform  precession 
Dipolar  magnetic  field  or  demagnetizing  field  for  spin  wave 
Exchange  field 

Exchange  field  constant  in  easy  plane,  normal  to  easy  plane 
Applied  (external)  dc  magnetic  field 

th 

The  applied  dc  field  above  which  the  n  order  spin  wave 
instability  is  no  longer  allowed  on  resonance 

Amplitude  of  applied  (external)  rf  pump  magnetic  field 

Apparent  critical  field 

Critical  field  or  threshold  field 

Value  of  h  at  beginning  of  pulse 

Value  of  h  when  the  apparent  threshold  is  reached 
m 

Maximum  rf  magnetic  field  amplitude  for  dielectric  resonator 
Applied  (external)  rf  pump  magnetic  field 
Nonzero  component  of  h 

Subscripts, each  of  which  takes  on  the  coordinate  values 
1,2,3  or  x,y,z 

Magnetic  and  phonon  coupling  constants  for  spin  wave  and 
elastic  wave  equations,  respectively 

Propagation  vector 

Elastic  wave  propagation  constant  on  resonance 
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n 


M 

N 

N 


n 


P 


P 

a 


^Vap 

P 

r 


P 

rm 


P 


Q 

^ext 

*L 

^0 

s 

<1 

B 

Bk 


Magnetization  vector 

Saturated  magnetization 

Effective  number  of  growing  phonon  modes 

Demagnetizing  factor  in  the  plane  of  the  ferrite  disk  « 

Order  of  spin  wave  instability  ■  2ta/o 

Number  of  turns  of  frequency  tuning  plunger  screw  out  from 
position  of  closest  approach  to  dielectric  resonator 

Peak,  average  power  absorbed  by  resonant  System 

Change  in  ferrite  absorption  due  to  growth  of  instabilities 

Incident  pulse  power  level 

Value  of  P^  at  beginning  of  pulse 

Phonon  power  absorption 

Phonon  absorption  at  apparent  threshold 

Reflected  pulse  power  level 

Change  in  reflected  power  level  due  to  AP^ 

Value  of  Pf  at  beginning  of  pulse 

Normalized  anisotropy  constant  *  H^Asii g 

Spin  wave  Q 

External  Q  of  resonant  system 
Loaded  Q  of  resonant  system 

Unloaded  Q  of  resonant  system;  ferroaMgnetlc  resonance  Q 

Phonon  or  elastic  wave  Q 

Frequency  norwLUzetion  parueter  -  V2^ 

Elastic  displacement  vector 
Elastic  wave  amplitude 
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£ 

£k#£k 


8,8 


k 


U#u 

X 


x,y,z 


x  ,y  ,z 


Spatial  coordinate  vector 

Time -independent  elastic  wave  amplitudes 

Standing  wave  ratio  at  resonance 

Exponential  phonon  growth  parameter 

Build-up  time  of  phonon  instabilities  to  the  apparent 
threshold  level 

Peak,  average  stored  energy  in  resonant  system 
Internal  magnetic  field  parameter  -  H/4xM# 

p 

Effective  Internal  magnetic  field  parameter  *  (H  +  Dk  )/4hM 

fit 

Crystallographic  coordinates  for  uniaxial  Znjf  ;  coordinates 
used  to  describe  dielectric  resonator  fields  d 

Coordinates  in  rotated  system  at  an  angle  with  respect 
to  fixed  system 


2 

rk 

2b 

a 


Unit  magnetization  vector  ■ 

Spin  wave  amplitude 
Uniform  precession  amplitude 

Spin  Vave  component  of  the  unit  magnetization  vector 
Reflection  coefficient  of  resonant  system 
Gyromagnetic  ratio  ■  ge/2mc 
Strain  tensor  component 

Effective  relaxation  frequency  in  the  presence  of  magneto¬ 
elastic  coupling 

Spin  wave  relaxation  frequency  •  o>/2Q  * 

Spin  wave  relaxation  frequency  «  u^/2Q 
Phonon  or  elastic  wave  relaxation  frequency 
Spherical  coordinate  used  to  specify  orientation  of  £ 
Density  of  ferrite  sample 
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♦o 

^Vopt 

a) 


% 

“>0 

°>p 


^0 

u>/Q 


Decay  time  constant  of  Incident  power  pulee 

Spherical  coordinate  used  to  specify  orientation  of  fc 

Angle  between  rotated  and  fixed  coordinate  systems 

Optimum  value  of  for  pumping  of  most  unstable  phonon  mode 

Frequency  of  unstable  spin  wave  or  elaatic  wave  mode 

Aniaotropy  field  frequency  parameter  -  yH^ 

Applied  dc  field  frequency  parameter  -  yH^ 

» 

Saturated  magnetisation  frequency  parameter  *  y4xMg 
Ferromagnetic  reaonance  frequency 

Pump  frequency  (i.e.,  frequency  of  applied  rf  magnetic  field) 
Spin  wave  resonant  frequency 

Spin  wave  resonant  frequency  containing  modulation  component 
Elastic  wave  resonant  frequency 

Elastic  wave  resonant  frequency  in  the  presence  of  mtgneto- 
elastic  coupling  and  pumping  by  the  uniform  precession 

Elastic  wave  resonant  frequency  in  the  presence  of  magneto- 
elastic  coupling  and  xero  pump 

Phenomenological  damping  parameter  for  spin  waves 
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was  varied  from  the  first  to  the  fifth  by  adjusting  the  position  of 
the  spin  wave  manifold  relative  to  the  pump  frequency.  The  nature 
of  the  eaturatlon  was  investigated  by  observing  the  waveform  of  the 
reflected  pump  pulse. 

Theoretical  on-resonance  spin  wave  thresholds  up  to  fourth  order 
and  off -resonance  spin  wave  thresholds  up  to  second  order  were  derived. 
The  first  order  theoretical  spin  wave  threshold  was  the  lowest  threshold 
in  the  range  of  low  dc  magnetic  fields.  At  higher  dc  fields  the  second 
order  theoretical  on-  and  off-resonance  spin  wave  thresholds  dominated 
the  higher  order  thresholds.  Taking  the  magnetoelastic  effect  as  the 
dominant  mechanism,  first  order  phonon  thresholds  were  derived  for 
transverse  and  longitudinal  elastic  waves  propagating  in  the  easy  plane 
and  along  the  c-axli .  The  lowest  theoretical  phonon  threshold  was  for 
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and  at  frequenclee  below  the  bottom  of  the  spin  wave  manifold,  this 
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wave  threshold.  Theory  was  developed  for  the  growth  of  phonon 
instabilities  from  thermal  level  in  order  to  relate  pulse  saturation 
effects  to  the  phonon  absorption. 

Transverse  pumping  thresholds  Observed  experimentally  were 
comparable  with  the  above  theoretical  predictions  and  strongly 
Indicate  that  first  order  phonon  instabilities  dominate  the  higher 
order  spin  wave  instabilities  at  frequencies  below  the  bottom  of  the 
spin  wave  manifold.  These  results  together  with  the  phonon  growth 
theory  made  it  possible  to  Obtain  rough  estimates  for  the  phonon  Q 
one  magnetoelastic  constant/  and  the  exchange  field  of  Zn^Y  in  the 
easy  plane.  The  observed  pulse  saturation  effects  and  the  results 
of  two  experiments  for  two  different  phonon  absorption  levels 
indicated  the  presence  of  nonlinear  elastic  wave  damping. 
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